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John T. Bolljahn, 1918-1960 


John T. Bolljahn was born in Oakland, Calif., on 
June 10, 1918, and passed away at Palo Alto, Calif., on 
June.25, 1960, in his forty-second year. He was the son 
of Herman Bolljahn, a construction engineer who had 
emigrated to California from Germany as a young man, 
and Lillian A. Smith, who was born in Glendale, Calif. 
He married Harriette M. Pedrotti on September 19, 
1942, and fathered three children, a son, Frederick, 
who passed away in 1948 at the age of 2 years, and 
daughters Mary and Claire, born in 1947 and 1951. He 
received his education at the University of California 
in Berkeley, from which he received the B.S. degree in 
May, 1941, and the Ph.D. degree in January, 1950, both 
degrees in electrical engineering. 

Dr. Bolljahn was on the staff of the Naval Research 
Laboratory, Washington, D. C., from August, 1941, 
through December, 1945, as Antenna Group Leader of 
the Aircraft Section, Radio Division. He returned to the 
University of California in January, 1946, to continue 
his graduate studies and to participate, as a staff en- 
gineer, in the antenna research program conducted at 
the University under the auspices of the Office of Naval 
Research. During this interval he served the University 
also as a lecturer in electrical engineering. In October, 
1949, he joined the newly-formed Aircraft Radiation 
Systems Laboratory of the Stanford Research Institute 
in Menlo Park, Calif., as assistant supervisor. He as- 
sumed a series of increasing responsibilities at the In- 
stitute, and was named Assistant Director of Engineer- 
ing Research in 1956. He retained an association with 
the Institute until his death. In 1956, with three others, 
he formed Granger Associates, serving as a director of 
that corporation and, in March, 1960, became its Execu- 
tive Vice President. 

Dr. Bolljahn was a member of Eta Kappa Nu, 
Tau Beta Pi, and Sigma Xi. He was named a Fellow 
of the [RE in March, 1960, “for his contributions to the 
theory and practical design of antennas.” He served at 


various times on scientific advisory committees to sev- 
eral government bodies and industrial groups. He was 
active in the affairs of the IRE in many capacities, and 
served on the Administrative Committee of the Profes- 
sional Group on Antennas and Propagation from 1953 
totlOS 7. 

Jack Bolljahn’s efforts in behalf of his profession 
were numerous, of exceptional quality, and of great sig- 
nificance to the advancement of the radio art. He made 
outstanding contributions to the theory of electrically 
small antennas and to practical methods for their design 
and test. He added important findings to the theory of 
electromagnetic scattering. His work on the properties 
of coupled, distributed elements led to fundamental ad- 
vances in the design of filters and transmission com- 
ponents. His understanding of the special problems of 
antennas for aircraft brought to every worker in that 
field a unifying clarification of the basic phenomena in- 
volved and a wealth of new, intensely practical, design 
approaches. He wrote 29 technical reports on a wide 
range of topics. Ten papers reporting his research were 
published in scientific journals, and he made many 
presentations to technical meetings. In addition, he was 
in the process of completing a number of additional 
papers. One of these “Synthesis of Modulated Corru- 
gated Surface-Wave Structures,” appears in this issue of 
the TRANSACTIONS. It is hoped that others can be pub- 
lished at a later date. 

Any outline of Dr. Bolljahn’s career and his con- 
tributions to the literature cannot begin to convey the 
breadth and depth of the man. His personal and pro- 
fessional influence extended beyond his immediate asso- 
ciates to hundreds of others to whom his wisdom, his 
analytical skills, his keen perception, his selfless dedica- 
tion, his warm and generous spirit, brought new appre- 
ciations of the challenge and the rewards of his profes- 
sion. It is a privilege to have known him and an honor 
to have been his friend. —J. V. N. GRANGER 
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The John T. Bolljahn Memorial Award 


By action of the Administrative Committee of the Professional 
Group on Antennas and Propagation, the annual prize for the best 
paper published in the Transactions of the Professional Group on An- 
tennas and Propagation is to be known as the John T. Bolljahn Mem- 
orial Award. The purpose of the award, in memory of the late Dr. John 
T. Bolljahn, is to recognize outstanding papers in the field of antennas 
and propagation by the award of a cash prize, and thereby to stimulate 
increased interest in the field. 

In order to make the amount of the award substantial enough to 
serve as a stimulus for the generation of meritorious papers, the $200 
presently allotted is to be increased by the proceeds from a fund sup- 
ported by contributions. In this way it is hoped that the honorarium 
could total $500 annually. A fund-raising committee has been ap- 
pointed by Dr. E. C. Jordan, Chairman of the Administrative Com- 
mittee. This Committee, acting on its own behalf and not as agents of 
the Institute of Radio Engineers, is in the process of soliciting con- 
tributions to this fund, which then will be given as a gift to the IRE for 
the purpose of the John T. Bolljahn Memorial Award. Co-chairmen of 
the fund-raising committee are as follows: 


M. Katzin 
5001 College Avenue 
College Park, Md. 


J. V. N. GRANGER 
974 Commercial Street 
Palo Alto, Calif. 
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contributions 


Synthesis of Modulated Corrugated Surface-Wave Structures* 


J. T. BOLLJAHN}, FELLOW, IRE 


Summary—An exact procedure for the design of a surface that 
will support a prescribed group of surface waves simultaneously is 
described. The surfaces resulting from the application of this pro- 
cedure are found to be modulated (i.e., to contain periodic waves in 
the direction of propagation) and to have surface impedance values 
that vary periodically in the direction of propagation. Finite-length 
surfaces of this type are hence reminiscent of Simon’s “cigar an- 
tenna.” Several surface designs are presented, and some general 
properties of the surfaces and the composite waves they support are 
discussed. Certain limitations on the design procedure are described. 


I. INTRODUCTION 


HIS PAPER describes a method for designing a 

modulated corrugated surface-wave structure to 

support a given set of surface-wave modes. Inter- 
est in this problem stems from recent work by Simon 
and Biggi,! Thomas and Zucker,? and Pease.® 


* Received by the PGAP, July 13, 1960. The work described in 
this paper was sponsored by the AF Cambridge Res. Ctr. under Con- 
tract AF 19(604)-3502 with Stanford Res. Inst., Menlo Park, Calif. 

} Stanford Res. Inst., Menlo Park, Calif. Granger Associates, 
Palo Alto, Calif., at the time of his death, June 26, 1960. 

1 J. C. Simon and V, Biggi, “Un nouveau type d’aérien et son ap- 
plication 4 la transmission de télévision 4 grand distance,” L’Onde 
Elec., vol. 332, pp. 1-16; November, 1954. (In French.) 

2 A. S. Thomas and F. J. Zucker, “Radiation from modulated sur- 
face-wave structures—I,” 1957 IRE NATIONAL CONVENTION REC- 
ORD, pt. 1, pp. 153-160. 

3 R. L. Pease, “Radiation from modulated surface-wave struc- 
aap tes 1957 IRE NationaL CONVENTION REcoRD, pt. 1, pp. 
161-165. 


Simon and Biggi’s work! describes the so-called cigar 
antenna, which is a disk-on-rod structure in which the 
disk diameters and/or spacings are varied periodically 
in the axial direction. Thomas and Zucker? consider the 
synthesis of desired radiation patterns using such an- 


tenna configurations, on the theory that the aperture ~ 


distribution along the antenna may be expressed as the 
superposition of a number of discrete modes, each in the 
form of a simple trapped wave. Pease* considers the 
approximate solutions for waves bound to a dielectric 
slab in which the dielectric constant is modulated pe- 
riodically in the direction of propagation. The latter 
two papers?® consider not only end-fire radiators, but 
also apertures which radiate in other than the end-fire 
direction. 

The approach used here is exact in the sense that it 
determines the characteristics (i.e., shape and surface 
impedance) of that surface (or family of surfaces) which 
satisfies the boundary conditions in the presence of any 
prescribed set of trapped waves. The analysis, like those 
in Thomas and Zucker, and in Pease, is applicable to a 
two-dimensional model. A more interesting solution 
from the practical design standpoint would be that for a 
cylindrical end-fire structure. Unfortunately, the tech- 
nique developed is not valid for the cylindrical geom- 
etry. This limitation is discussed in the Appendix. 
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A further limitation is that the analysis is inherently 
applicable only to the corrugated-surface type of struc- 
ture and not to the dielectric slab form of trapped-wave 
surface. This limitation arises since the analysis postu- 
lates that there is no energy flow across the interface 
between the guiding structure and the free-space region 
above this structure. In a dielectric slab radiator sup- 
_ porting more than one mode, energy can travel in the 

_ direction of propagation both within the slab and in 
the medium above the slab, and a periodic interchange 
of energy can occur between the two regions across the 
interface. The restriction that the interface be a surface 
across which no energy passes is therefore not applicable 
in this case. In a corrugated-surface structure, there is 
no energy flow below the interface, and hence the re- 
striction is a valid one. 


IJ. ANALYsIs 


We shall confine attention to the two-dimensional case 
(Fig. 1) for which we may represent the simple trapped 


E 
y DIRECTION OF 
> PROPAGATION 


x 


Fig. 1—Coordinate system. 


wave by the field components, 


1 ce e7 ey iBx 
we 
Beg 
Es = — — gy ibe 
we Qa 
H, = — eB 
a 
= B? — 2, (1) 


We wish now to postulate a group of such waves having 
different a’s and #’s, but the same frequency, and to 
find the flow lines represented by the directions of the 
‘real part of the complex Poynting vector. These flow 
lines or surfaces are the surfaces across which there is no 
energy flow, and are hence the surfaces which, if de- 
signed to present the proper surface reactance, can sup- 
port the set of waves postulated. The imaginary part of 
the complex Poynting vector can be used, as described 
later, to define the surface reactance required along the 
surfaces defined above. 
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Q37 
For a group of m waves, we may write 
| OOS w Ga Vee Des euler AY 3 
1 n n 
= d, — Je Hs e7any—Bnx . ey Bn —any+IBnx 
WE n An n An 


. B, 
— dy ae » Bye 7 onda * De — e-anytiBnx (2) 


WE xn n An 
Separating real and imaginary parts,‘ 


we Re (E X H*) 


n B 2 
= | (22) bare 
p=1 \OAp 


sD 
+o 


p=1 q=1 OpAq 


Lae oP) 

- pPaq 
+ ay Dds 
p=1 q=1 ApAg 


we Im (E X H*) 
© Ct BB, 


=4.5 


p=1 q=1 ApAq 


*B 
+a] ¥ — e-Papy 


p=1 Ap 


(8, + B,)e-‘ert#0¥ cos (Bp — a) 


(aq — ape “ren U sin (65 18,)a. (3) 


(Bq ao: Spescrtcan sin (Bp = Bq)% 


B,Ba 


‘2353 


p=1 q=1 pq 


(ap + ape otaev cos (By — ae. (4) 


Specializing to the case n=3, and taking (6:—B2) 
= (8:—83) =4, the expressions become 


we Re (E X H*) 


Bx Be Bs 
— Tike |: ae ea + B2 ae e7 2e2u +. B;? Soe —2a3sy 


Qa, Q2 Q3 


+ cos 6x | BB. (Bi + Bs) 


e7 (e1ten)y 
a1a2 


+ BoB; (B2 + Bs) cesta | 


A2as 


+ cos 26x Ez (8s + 63) ciateon | 


aA1A3 


(a2 = a1) e7 (aitae)y 


aia? 


+ dy {sin bx Ez 


+ BoB; (a3 — a2) eter | 


A203 


(as — a) eevtenn | (5) 


A1AsZ 


+ sin 26% | BBs 


4 A more general formulation may be made in which the relative 
phases of the various modes at x =0 are arbitrary. The cases consid- 
ered in this report, however, are amply covered by the more limited 
(and compact) expressions given. 
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we lm (E x H*) 


Bi — B) 
=—4, fsin 6x Ez ia al en (arto2)y 


a1a2 


+ By B; (B2 — Bs) eietens| 


A2Q3 


+ sin 26x [ain eintaon | 


A1a3 


BY B? B;? 
— dy SE ob —— @ 7 2a2y + —— ¢@ 2a 


eal a2 as 


a1 + ae) 
+ cos 6x | BB: ame en (artas)y 


a1a2 


+ eo costes | 


aA2As 


+ cos 26x Ez (a + 3) cintenn lh e (6) 


a1a3 


To calculate the surfaces across which no energy 
flows, we note that since there are no sources of energy 
in the region under consideration, 


V-Re (E X H*) = 0. 


Assuming that the surfaces sought are of the form shown 
in Fig. 2 (there will actually be a family of such surfaces 


Fig. 2—Application of Gauss’s theorem to surface. 


corresponding to different values of yo), we may apply 
Gauss’s theorem to the region shown cross-hatched and 
obtain 


dx 


y=2Xo 


Yo 
_ f Re (E X H*), 
0 


dy -{ Re (Bex); 
z=0 0 


+ { Rex H, 
0 


L=29 


Substituting from (3) and (4) and carrying out the 
integrations, we find that 
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ee 
. 2a;* 


= cos 6x |= (===) e(artaa)y 


Q1Q2 \Q1 ai a2 
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+ cos 26x BBs (2 *) en (airtaa)y 
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= B, eek. e2e1v0 Be Bs e240 B;? Bs 
2 3 


2013 ae a3 


+ BBs (Ao ge (artaz)y0 


12 \ ay + a2 


= BBs (=) ge (estas) yo 


ara; \ a2 + a3 


+ BiB; ==") e7 (aires) yo | 


13 \ Ay fe Q3 


Bo ! Bs 
—2oiy + Be? rere + Bs? ae 7 2as 
Q2 a3 


e7 22340 


(8) 


For a given value of yo, the right-hand member of this 
equation becomes a constant. The desired surface for 
that value of yo is determined by substituting values of 
x into the equation and finding the corresponding values 
of y which lead to satisfaction of the equation. 

Up to this point, the surfaces discussed have been 
defined only as mathematical surfaces with fields both 
above and below them. If we wish to place a physical 
boundary on one of these surfaces which will retain just 
the fields above the surface, it is clear that this bound- 
ary must present a surface reactance which varies from 


point to point in accordance with the ratio of the tan- — 


gential E field to the tangential H field along the bound- 
ary. Thus 


E, Im(E X H*) X Re (E X H*) 


X=— = 
Hy | H|?| Re (EX H*)| 


(9) 


Surface shapes and reactance functions have been 
calculated for a number of different wave groups, and 
these are shown in Fig. 3(a)—3(d). It is noted that the 
first of these is designed for a group of three waves while 


the others are all for two-wave groups. The basic objec- — 


tive of an end-fire antenna employing such a guiding 
surface is to permit the launching of a wave having v/c 
slightly less than unity on a surface having a more tightly 
trapped-wave system than is obtained on a single-mode 
surface for the same value of v/c. It is noted in these 
curves that in most cases the average value of surface 
impedance is considerably higher than that of the sur- 
face required to support the fastest wave alone. Since 
the attenuation of the fields normal to the surface is, to 
a good approximation, proportional to the local surface 


450 Vr SAV 
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(b) 
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(d) 


Fig. 3—Modulated corrugated-surface design. 


reactance, it follows that the objective of achieving a 
relatively tightly bound wave system is achievable with 


this type of structure. The surface reactance reaches a 


maximum in most cases in the troughs of the modulated 
surfaces. It is seen, therefore, that in these cases the 
feed system for such an antenna should be placed in a 
trough in order to couple most effectively to the wave 


" system to be excited. This is not a general result for all 


surfaces, however, as is shown by the results in Fig. 3(d). 


III. Some PROPERTIES OF SURFACES AND FIELDS 


A question of interest in connection with the design 
of such surfaces is: Can a plane-corrugated surface (i.e., 
one without undulations) be designed to support a finite 


group of trapped waves by simply varying the surface 
impedance in the direction of propagation? It is readily 
shown with the formulation used here that the answer 
is no. Referring to (3), the y component of Re (EX H*) 
for a group of m waves at a fixed value of y is seen to be 
of the form 


Re (E x H*), eee 


Pree Bee. ; 
= 3 pS (aq — ap) sin (6, — By)x, (10) 
WE p=l g=l1 pq 
where 
B, = Bie", 
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and 


i eis. — 
Be = Boert2!. 


If n» is finite, this expression is a finite series in 
sin (8,—8,)x. Since there is no combination of coeffi- 
cients which can reduce this expression to zero for all x, 
it follows that there is no finite group of waves for which 
the required surface is planar. 

Another point of interest in fields of the type con- 
sidered here is the occurrence of regions in which the x 
component of Re (EX H"*) reverses direction. The fact 
that such regions occur is easily demonstrated in the 
case where only two waves exist by considering the 
phase relationships between the total fields EZ, aud is 
first at x=0 and then at x=7/6. At x=0, the two E, 
and the two H, components are in phase. At x=7/6, 
because of the different phase velocities, the two Ey, 
components are 180° out of phase, as are the two H, 
components. Since the ratio E,/H, is different for the 
two waves, there will always be a range of relative 
amplitudes of these two waves for which the resultant 
value of E, at x=7/6 is 180° out of phase with the re- 
sultant value of H, at the same point. This means that 
the x component of Re (EX H%*) is negative at x=7/6 
and for values of y for which the relative wave ampli- 
tudes are within the proper range. Since the relative 
wave amplitudes will go through all possible values as a 
function of y, with the wave having the larger a becom- 
ing negligible at large positive values of y and vice versa, 
it is seen that regions of negative Re (EX H*), will oc- 
cur with any pair of waves. 

The characteristics of the flow lines of Re (EX H*) in 
the vicinity of this reversal of the x component has 
been investigated in detail for the two-wave case in 
which By=B.=1, 1:/c=0.975, v2/c=0.852. The flow 
lines for this case are shown in Fig. 4. Of particular inter- 
est is the fact that some of the lines close on themselves, 
indicating that stored energy in the form of circulating 
waves exists in these regions. 


IV. CONCLUSIONS 


It has been shown that an exact procedure exists for 
the design of a modulated corrugated surface to support 
a specified group of surface waves. In the illustrative 
examples of the application of this procedure, it has 
been shown that the surface-impedance level on such a 
surface attains values well in excess of the correspond- 
ing values for the single-mode surfaces designed for the 
fastest wave in the group. This supports the point of 
- view that a modulated surface can contain a wave with 
slightly less than the velocity of light as a part of a rela- 
tively tightly bound wave system. 

The analysis used is limited in application to two- 
dimensional corrugated surface radiators. It is further 
limited in that it provides no information on the band- 
width characteristics of the structures designed. 
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Fig. 4—Energy-flow pattern in vicinity of 
region of energy reversal. 
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Fig. 5—Cylindrical coordinates system. 


APPENDIX 
CYLINDRICAL RADIATORS 


The end-fire type of surface wave on a cylindrical 
surface is described by a combination of the first-order 
TE- and the first-order TM-cylindrical wave functions,® 


oy We 
E, = — jhay, oe eae paris 
or 
h oy 
Eo = — — ayy + jopb: — 
Yr or 
E, = (R* — W)awn 
2 
re) 
in ee 
wuUr r 
—jk? oy h 
Hy = he 
ayn or Yr 
H, = (Rk? = h*)bw, (11) 


55), As Sena tton, "hs Theory,” McGraw-Hill Book 
Con ines, New York, N. Y.; 1941 Z ak eee 


1961 


where 
Yi = cos 0H, (./k? — h? r)e-ihe, 


To calculate the wave constants appropriate to a 
given cylindrical corrugated surface (it is assumed that 
there are many corrugations per surface wavelength), 
one first chooses a relationship between a; and }; such 
that He=0 at r=a, and then selects the value of h 
which causes F,/H, to reduce to 7X;, where X, is the 
surface reactance. 

In principle, it is possible to employ two or more 
cylindrical surface-wave modes having different a and 6 
constants and different values of # in the same way that 
the planar modes were employed in this paper. Thus, 
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using the analogies L,-E,, E,—-E,, and H,—-H,, it 
should be possible to define surfaces analogous to those 
developed here. When this procedure is carried out, 
however, it is found that the boundary condition E,=0 
on the surface is not satisfied in general by the wave 
groups used to generate the surface function, and hence, 
the synthesis procedure is not applicable in this case. 
This difficulty would not arise in the axially-symmetric 
TM solutions, and an analogous formulation could be 
developed to provide the characteristics of multimode 
cylindrical structures for such modes. The axially- 
symmetric TM modes are of no interest for end-fire 
antenna applications, however, since they produce radi- 
ation-pattern nulls in the axial direction. 


Magneto-Ionic Faraday Rotation of the Radio Signals 
on 40 Mc from Satellite 1957. (Sputnik I)* 


E. V. SORENSEN, ASSOCIATE MEMBER, IRE 


Summary—Some of the fadings of the 40-Mc signals from Sput- 
nik I during morning transits are assumed to originate in Faraday 
rotation in the ionosphere. This has been investigated for the transit 
of October 12, 1957, 5:30 GMT, for which the orbit was known with 
reasonable accuracy. 

The fading records show an increasing fading rate during the last 
part of the transit. This is contrary to the simple theory based on the 
assumptions of a flat earth and a homogeneous ionosphere, from 
which a constant fading rate should be expected. The main reason 
why simple theory did not apply in this case, is that the satellite 
travels through the twilight zone from night towards day, i.e. from 
low to high ionization. 

The Faraday rotation is calculated by means of data for the orbit, 
magnetic field maps and predictions of the ionosphere conditions. 
The result is in good agreement with the observation when a high 
ionization between the altitude of maximum ionization and the satel- 
lite is assumed. The effect of refraction is estimated by means of the 
perturbation theory. It increases the Faraday rotation from 0 to 19 
per cent during the transit. The general conclusion is that, even 
though the fading behavior of a transit near sunrise is explainable to 
some extent, it is still difficult to deduce detailed information on the 
state of the ionosphere from it, especially when the frequency is so 
low that refraction becomes significant. 


* Received by the PGAP, May 19, 1960; revised manuscript re- 


ived, September 20, 1960. ; 
at aetwave Lab. of the Danish Academy of Tech. Sciences, 


Copenhagen. 


I. INTRODUCTION 


T has been shown! that a rather simple inter- 
| pretation of the fading rate due to magneto-ionic 

Faraday rotation of plane-polarized radio waves 
emitted from a satellite is possible when the following 
approximations apply simultaneously. 


1) The curvature of earth can be neglected in the 
geometry. 

2) The satellite travels at a constant height with a 
constant velocity. 

3) The height-profile of the electron density is inde- 
pendent of the geographic position between re- 
ceiver and satellite. 

4) The magnetic field is homogeneous. 

5) The signal frequency is so high that refraction is 
negligible. 

6) The quasi-longitudinal condition is fulfilled, z.e. the 
angle between the magnetic field and transmission 
path differs from 90 degrees by at least a few de- 
grees. 


1S, A. Bowhill, “The Faraday-rotation rate of a satellite radio sig- 
nal,” J. Atmos. and Terrest. Phys., vol. 13, pp. 175-176; December, 
1958. 
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|Fig. 1—Fading record for the transit October 12, 1957, 5:30 GMT. Trace A, Horizontal dipole. Trace B, Vertical monopole. 
Bandwidth of receivers: 145 cps. Amplitude response: nearly logarithmic. 


The fading rate will then be constant in time, inde- 
pendent of receiving site and proportional to: 

1) the total amount of electrons in a vertical column 

below the satellite, 

2) the magnitude of the magnetic field component 

along the satellite orbit, 

3) the ratio between velocity and height of the satel- 

lite. 

For some satellite transits the fades are satisfactorily 
explained by this simple picture. In other cases, where 
not all the approximations apply simultaneously, the 
fading behaviour is more complicated, although the 
fades in all probability originate in Faraday rotation. 

Some of the recorded transits of the first Russian 
satellite 1957a (Sputnik I) launched October 5, 1957 
belong to the last category and have given rise to some 
discussion.? Such a transit will be analyzed in the fol- 
lowing. The Faraday rotation will be calculated point by 
point on the basis of available data for the orbit, the 
magnetic field and the ionosphere. 

The only starting approximation is the assumption of 


quasi-longitudinal conditions. This will be justified — 


later. The calculated values will be corrected approxi- 
mately for the effect of refraction. The results agree 
fairly well with the observations for large values of the 
electron density between the height of maximum ion- 
ization and the satellite. Especially the approximations 
(3) and (5) have proved to be serious over-simplifica- 
tions in the present case. 


2 “Radio observations of the Russian satellites, a discussion meet- 
ing,” Proc. IEE, vol. 105, pp. 81-115 (particularly pp. 101-104 and 
113); March, 1958. 
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Fig. 2—Number of fades per minute vs time. 


Blackband, et al.’ have demonstrated how absolute 
information on horizontal gradients in the ionosphere 
may be obtained from such transits. This requires that 
the point corresponding to quasi-transversal propaga- 
tion can be located on the fading records, a condition 
which is not always fulfilled. 


II. OBSERVATIONS 


An investigation of the fading characteristics of the 
40-Mc radio signals from 1957a@ (Sputnik I) as received 
in Copenhagen revealed the following: 

1) During the early transits in the evening and the 
first part of the night, when the satellite passed in alti- 


’V. T. Blackband, B. Burgess, I. L. Jones, and G. J. Lawson, 
“Deduction of ionospheric electron content from the Faraday fading 
of signals from artificial earth satellites,” Nature, vol. 183, pp. 1172- 
1174; April, 1959, 
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tudes of about 200 kilometers and in directions from 
south-west to north-east, the signals faded with the 
constant rates of 6.8 and 13.6 fades per minute inde- 
pendent of the satellite position. These fades are as- 
sumed to originate in the satellite spin, and accordingly 
they are dismissed as being unimportant from the view- 
point of propagation effects. 

2) The late transits in the last part of the night and in 
the morning, when the satellite passed in altitudes be- 


tween 400 and 500 kilometers and in directions from . 


north-west to south-east, showed fades of gradually in- 
_ creasing rate. In general, the fading started about the 
_ time of nearest approach with a rate of some 5 fades per 
minute and lasted until the signals disappeared, at 
which time the rate was about 50 fades per minute. 
Simultaneous recording of the signals in a vertical mono- 
pole and a horizontal dipole revealed that the fades 
occurred nearly in time quadrature on these antennas, 
as was to be expected for the Faraday rotation type of 
fading, when the transit takes place near the horizon 
such as here. 

Fig. 1 shows such a record of the latter half part of the 
transit October 12, 1957, 5:30 GMT (the lines marked 
“tune” indicate that the common local oscillator for the 
two receivers had been adjusted manually to compen- 
sate for the Doppler-frequency shift). In Fig. 2 the 
Faraday rotation fading rate, as found by an examina- 
tion of the record, is plotted against time. Sudden jumps 
amounting to one or two times the RPM of the satellite 
spin may occur in such a curve. In the present case 
however, the fading rate varies in a fairly continuous 
way, indicating that the fades almost certainly originate 
in Faraday rotation. 

Fig. 3 shows the corresponding part of the orbit as 
established from scattered information. The different 
points with the corresponding indications of time have 
been used in the calculations of the Faraday rotation. 

If the assumption that these fades originate in Fara- 
day rotation is correct, the variation with time of the 
phase difference between the ordinary and the extraor- 
dinary ray may be obtained from graphical integration 
of the fading rate g(t) as indicated in Fig. 2. 


o(t) — o(to) = an f g(i)dt_ rad. (1) 


If we choose fo as the time when the ray path is per- 
pendicular to the magnetic field, #(to) will be very close 


_ to zero. The corresponding satellite position is somewhat 


to the north of point A in Fig. 3. This time cannot, how- 
ever, be found from the fading record as it corresponds 
to an infinitely long fading period. Accordingly we have 
‘chosen the time (5:30.19 GMT) from which the fade 
can be traced with reasonable certainty. This corre- 
sponds to point A in Fig. 3, for which the angle between 
the ray path and the magnetic field is 67 degrees. The 
result of the integration is plotted in Fig. 4 together 
with the calculated curve. 
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Fig. 3—Assumed orbit for 1957 during transit 
October 12, 1957, 5:30 GMT. 
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Fig. 4—Variation with time of the phase difference ¢ between or- 
dinary and extraordinary rays. Curve A, obtained from ob- 
served fading rate. Curve B and point C, calculated by means of 
predicted data for the magnetic field and for the ionosphere. 
Different N(h) profiles. Curve D, calculated curve corrected for 
the effect of refraction. 
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Ill. THEorY 


As shown by Appleton in 19324 the refractive index 
of the ionosphere, taking the magnetic field of earth into 
consideration is given by: 


m2 = 1— - : ; (2) 
nie en A VE See 
a1 (a — 1)? 
where: 
mo” wH cos 6 wH sin 6 
ee ane ie or 


H is the local magnetic field in oersteds 
6 is the local angle between the ray path and the 
magnetic field 
N is the local ionization in electrons per cm? 
e is the charge of an electron: —4.8-107!° esu 
m is the mass of an electron: 9.11-10~*8 g 
c is the velocity of light: 310!° cm/sec 


If the quasi-longitudinal condition holds, (2) may be 
reduced to 


(4) 


where mp is independent to the magnetic field. 

In the present analysis, the “worst” value of 8 occurs 
in point A Fig. 3 (0=67 degrees). The magnetic field is 
here 0.4 oersted. If an ionization of N=4-10° electrons 
per cm® is assumed, and the frequency is 40 Mc, we 


find: 


Calculated from (2) | Calculated from (4) 


1—n 1.024-10 1.018-10 
a 1.002-102 0.996-10 
An =n —N2 Ag MYON OR Dslevoules 


which indicates that (4) is a reasonable approximation 
in the following analysis. 

Under the same condition the two magneto-ionic 
components of the radio wave from the satellite will be 
nearly circularly polarized TEM fields rotating in op- 
posite directions. When they arrive at the receiver, the 
phase difference between them will be: 


a R R 
as lf nds; — i madss| rad. 
G s Ss 


As the satellite moves, ¢ will vary with time. Each 
alteration amounting to 27 will turn the plane of polar- 
ization of the resulting received signal by 7 radians and 
produce one fade. The function @(¢) is determined from 
calculation of ¢ for a suitable number of points of the 


a (5) 


4E. V. Appleton, “Wireless studies of the ionosphere,” J. JEE, 
vol. 71, pp. 642-650; October, 1932. 
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orbit (points A-G, Fig. 3). From this the quantity p(t) 
—(to) may be compared to the values found from the 
graphical integration outlined in Section II. 

A rigorous calculation of the integrals in (5) is out of 
the question as it requires the tracing of the two rays 
through a medium with an inhomogeneous magnetic 
field and electron distribution. An approximate solution 
may assume the following form, where (5) is rewritten 
toe 


w R wo CR 
o~ PA f= m)ds= F=f nds. (0) 
CvYs5 Cvs 


Here the integral is taken along the optical path be- 
tween S and R, and the refraction is taken into con- 
sideration through the approximate correction factor F 
(Section V and Appendix). 

Substitution of (3) and (4) in (6) gives 


e? 


ac?-m*:f? 


o~F f _N) -H(s)-cos6(s)ds. (7) 


The integral in (7) is convenient for graphical solution 
when the variations of NV, H and cos @ along the optical 
path are known. The establishment of the electron 
density data, for which there is considerable uncer- 
tainty, will be outlined in the following. The mapping 
of the two other quantities rests upon charts of the geo- 
magnetic field® and the following information on the 
orbit: 


1) The inclination of the orbit is about 65 degrees. 

2) The transit in question is an overhead transit for 
the Royal Aircraft Establishment radio station in 
Lasham, England.® Incidentally this point is al- 
most coincident with transit point as seen from 
Copenhagen. 

3) The vertical motion of the satellite was derived 
from a curve giving height and latitude vs the 
time since the preceding northward equatorial 
transit.” 


IV. THE IONOSPHERE CONDITION 


The first assumption in establishing the N(s) data is 
that only the F, layer is present. This is very probable 
in view of the hour of the day. 

The maximum electron density Nm and the corre- 
sponding altitude h» are mapped by means of charts of 
the predicted fo/:—and MUF (4000 km)—frequencies 
for October 1957, (Slough Bulletin A—No. 128), in con- 


5 E. H. Vestine, L. Laporte, C. Cooper, I. Lange, and W. C. 
Hendrix, “The Description of the Earth’s Main Magnetic Field and 
Its Secular Change, 1905-1945,” Carnegie Inst., Washington, D. C., 
No. 578; 1947. 

_ ® Radio observations of the Russian satellites, a discussion meet- 
ing, op. cit., p. 94. 

TH, Uyeda, T. Ishida, H. Shibata, and M. Mambo, “On the re- 
ception of radio waves from Russian Earth Satellite I,” J. RRL 
(Japan), vol. 20, pp. 135-141; April, 1958. 
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Fig. 5—Dependence of the MUF factor Xq on height k,, of maximum 
electron density and semi-thickness Y,, of a parabolic F, layer 
(Appleton and Beynon). 


nection with Appleton-Beynon’s theory on MUF cal- 
culation.® NV, is given by foFs: 


Nm = 1.24-104(foF2 Mc)? e/cm? (8) 
hm is found by means of the MUF-factor defined by 
MUF(d) 
tq = ————— - (9) 
SoF 2 


Fig. 5 (taken from Appleton and Beynon’) shows that 
this quantity is a sensitive function of hn, whereas it is 
fairly independent of ym/ho where ym is the semi-thick- 
ness of the parabolic V(t) profile assumed in the theory, 
* and ho is the height of the lower boundary above the 
- ground. The values thus found of NV, and hm below the 
satellite and above Copenhagen during the fading- 
period of the transit are shown in Fig. 6. The motion 
of the satellite through the twilight zone from night 
towards day explains the increasing value of Nm, and 
this again is the principal reason for the increasing fad- 
ing rate characterizing the morning transits of this satel- 
lite as observed from Copenhagen. 

The knowledge of NV, and h» is not sufficient to ob- 
tain N(s). We still need to form an idea of the remain- 
ing part of the N(h) profile, and here the greatest 
amount of uncertainty is met. According to Appleton 


8E. V. Appleton and W. J. G. Beynon, “The application of 
ionospheric data to radio communication problems: Part I,” Proc. 
Phys. Soc,, vol. 52, pp. 518-533; July, 1940; “The application of 
ionospheric data to radio communication problems: arte wel roc. 
Phys. Soc., vol. 59, pp. 58-76; January, 1947. 
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Fig. 6—Variation with time of N,, and h,, below satellite, as derived 
from ionosphere predictions for October 1957. The values over 
anny aoa during the transit were: N,,~4.3-4.6 105 e/cm, 

We m. 


and Beynon’ a fairly good approximation to the M(h) 
profile below h» is a parabola with yn/ho~0.5. 

We know as yet very little about the electron distribu- 
tion above hm. An extension of the parabola above hn» 
gives too few electrons to account for the magnitude of 
Faraday rotation observed. This is also the case if we 
make use of a Chapman distribution!® with usual values 
of the scale height.!! This is consistent with the few 
rocket measurements of the upper ionosphere reported 
so far! and with recent theories of whistler propaga- 
tion!® and the generation of low frequency signals in 
the exosphere,!* which all indicate that the electron 
density above hm decreases more slowly with the height 
than previously assumed. In view of the relatively low 
satellite height we therefore make the convenient as- 
sumption that the electron density is constant and equal 
to NV, between h, and the satellite. As shown in Fig. 4 
this gives a fairly good fit to the observed Faraday rota- 
tion. Blackband, et al.* have used the same approxima- 
tion in their work, and their results proved to be con- 
sistent with vertical incidence soundings. Fig. 7 shows 


® Appleton and Beynon, op. cit., vol. 59, p. 65. 

10S. K. Mitra, “The Upper Atmosphere,” The Asiatic Soc. of 
Calcutta, India; 1952. 

1 The result of [$(¢)—¢#(to)]X1/F where ¢ corresponds to point 
G. Fig. 3 is shown in Fig. 4 for the case of a Chapman distribution 
with scale height H=50 km. The distribution is given by 


Ny = Nmei!?2-2 — e®) 
where 
h — ln 
a : 


2 V. I. Krassovsky, “Exploration of the upper atmosphere with 
the help of the third Soviet Sputnik,” Proc. IRE, vol. 47, pp. 289- 
296; February, 1959. 

13 R. A. Helliwell and M. G. Morgan, “Atmospheric whistlers,” 
Proc. IRE, vol. 47, pp. 200-208; February, 1959. 

“4 R. M. Gallet, “The very-low-frequency emissions generated in 
the earth’s exosphere,” Proc. IRE, vol. 47, pp. 211-231; February, 
1959. 
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Fig. 7—Examples of the functions N(s) and H(s)-cos 6(s). 
—_——_—_—— WN(s) according to N(h)—profile I (Fig. 4) 
------ N.(s) according to N(h)—profile II (Fig. 4) © 
—— ~— —— Approximation to N,(s) for calculation of refraction. 


the variation of NV and H cos 6 along some of the integra- 
tion paths. The last quantity is remarkably constant 
along each path, which means that the Faraday rota- 
tion is a measure of the amount of electrons between 
receiver and satellite. 


V. EFFECT OF REFRACTION 


A perturbation calculation of the refraction factor F 
in (6) is given in the Appendix. The underlying assump- 
tions are: 


1) The refraction is so small that the values of NV, 1 
and cos @ along the optical path can be used. 

2) The product NA cos @ is constant along the iono- 
spheric part of the optical path. The validity of 
this approximation may be judged from Fig. 7 
path R-G. 

3) The tilt of the ionosphere boundary relative to 
earth is negligible around the penetration-point 
of the ray. 


The refractive index of the ionosphere along one of the 
transmission paths under consideration will now be a 
constant depending only on the satellite position. The 
result is 


a 


F~wi+ (1 — no) cot? u (10) 


a 

where mp is defined by (4); a and 0 are the lower and 
upper parts respectively of the optical path as divided 
by the ionosphere boundary, and wu is the inclination of 
the optical path at the penetration point (see Fig. 8). 
F was found to increase from approximately 1 to 1.19 
during the transit. The corrected curve for (¢(t) —¢(to)) 
is shown in Fig. 4. 
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Fig. 8—Geometry involved in calculation of 
refraction correction factor. 


VI. COMMENTS 


As may be seen from Fig. 4, the calculated total num- 
ber of fades in the period under consideration agree 
very closely with the observed number. The slope of the 
curves, based upon calculation or observation, reflecting 
the instantaneous fading rate differ somewhat more, 
especially towards the end of the period. This may be 
due to local deviations from the predicted ionosphere 
conditions. 

The results indicate a relatively high ionization above 
the height of maximum ionization. 

The effect of refraction is far from negligible even at 
the relatively high frequency of 40 Mc. In our case it 
increases the total number of fades by 19 per cent. 

The instantaneous fading rate cannot be taken as a 
measure of the electron content in a vertical column 
below the satellite (as predicted by simply theory) when 
the ionization varies along the orbit. If the refraction 
can be neglected, it is still a measure of the number of 
electrons between the satellite and the observer. How- 
ever, this condition is not sufficiently fulfilled at 40 Mc. 

In view of this examination we conclude that detailed 
information on the ionosphere can hardly be expected 
from fading records corresponding to transits near sun- 
rise, unless higher frequencies are used or more refined 
and sensitive methods of analysis are developed. 

Finally, reference should be made to additional theo- 
retical and observational work on this subject. Aitchison 
and Weekes’ have gained ionosphere data from Fara- 
day rotation rate measurements at 20 and 40 Mc for 
1957a2. Their analysis includes a vertical component 
of the satellite motion, and the effect of refraction is 


_ _% G, J. Aitchison and K. Weekes, “Some deductions of ionospheric 
information from the observations of emissions from satellite 1957 
rata J. Atmos and Terrest. Phys., vol. 14, pts. 3 and 4, pp. 236-243; 
1959, 
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estimated. They use a horizontally stratified ionosphere 
model and exclude transits near sunrise. 

Garriott!® has used two techniques to obtain electron 
density data from radio observations of Sputnik III 
- during a period of eight months. One method involves a 
_ determination of the total rotation angle at the transit 
F point. The other is based on rotation rate measurements. 
A complete account of all the approximations is given 
and the relevant corrections are made. Horizontal 
gradients appear to be included in their computation 
program. 

APPENDIX 
CALCULATION OF THE EFFECT OF REFRACTION 


The geometry of the problems is shown in Fig. 8. The 
ionosphere-approximation applied in this model is ex- 
plained in Section V. 


General Method 

First the magnetic field is neglected and the electrical 
distance WV in radians between satellite and receiver is 
expanded in terms of the variable: x=1—mo: 

The magnetic field then is introduced as a perturba- 
tion, and the phase difference between the two resulting 
rays is calculated as 

o= AY = — ——— An. (11) 
The electrical length of the refracted ray path A—B 
-in Fig. 8 is 


= = (4 BU — x). (12) 


The ray segments A and B may be expanded in terms 


16 Q, K, Garriott, “The determination of ionospheric electron con- 
tent and distribution from satellite observations,” J. Geophys. Res., 
vol. 65, pp. 1139-1158; April, 1960. 
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of the infinitesimal angle A@9 by means of simple geo- 
metrical relations: 
Avene AG ear Ag? 
B=6+ LAd+ NAP+:::. 


(13) 
(14) 


I 


A@ in turn may be expanded in terms of x =1—m by 
means of the law of refraction: 


Ag = Px + Ox? + + 
substitution of (12)—(15) in (11) gives 


LP.=(M + N)P? 
b= [142 ; % 


(15) 


ps |= ban, (16) 
Cc 


The correction factor F in (6) is now recognized as 
the term in the brackets. Note that the first order ap- 
proximation to ¢ requires both the linear and the quad- 
ratic terms of (13) and (14), but only the linear term of 
(15) as Q is absent in the coefficient to x. This is a con- 
siderable simplification in the detailed calculations. The 
result of the somewhat comprehensive procedure is 
(referring to Fig. 8): 


o~ E + (1 — mo) a Cots u| oan. (17) 


Note that the curvature of earth is absent in this ap- 
proximation. 
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On the Guided Propagation of Electromagnetic Wave Beams* 


G. GOUBAUT, FELLOW, IRE, AND F. SCHWERINGT{, MEMBER, IRE 


Summary—Any field in a half-space can be described by a con- 
tinuous spectrum of cylindrical waves. If this spectrum comprises 
substantially only waves whose propagation constant is very close to 
the plane wave propagation constant, the field can be resolved into a 
set of elementary wave beams which are characterized by Laguerre 
polynomials. They satisfy orthogonality relations like the wave 
modes in a waveguide. The elementary beams or “beam modes” can 
be reiterated and guided by reconstituting the cross-sectional phase 
distribution at certain intervals. Reiterative beams are utilized in the 
beam waveguide. The finite size of the phase resetting devices effects 
a modification of the reiterative beam modes and causes diffraction 
losses. These losses decrease very rapidly with increasing diameter 
of the phasing devices. 


INTRODUCTION 
Pilea wave beam emitted by a highly directional 


antenna, such as a parabolic antenna, has within 

the Fresnel region a substantially uniform diam- 
eter. However, the cross-sectional amplitude and phase 
distributions vary substantially along the path of the 
beam.! These variations depend on the field distribution 
at the antenna. The question may be asked whether 
there are beams in which the cross-sectional amplitude 
distribution is repeated at a certain distance from the 
origin of the beam. If there are any such beams, the 
original phase distribution can be reconstituted by suita- 
ble phase shifting means so that a new Fresnel region is 
formed which has the same field configuration as the 
original one. Repeating this iteration process by passing 
the beam along a structure of uniformly spaced phasing 
devices (phase transformers), the beam is guided with- 
out expansion of energy. 

The main objective of this paper is to show that there 
are indeed reiterative wave beams. A further objective is 
the calculation of the diffraction loss caused by the 
finite dimensions of the phase transformers. 


REITERATIVE WAVE BEAMS 


We consider wave beams which propagate along the 
z-axis of a cylindrical coordinate system, Fig. 1. The 
general solution of Maxwell’s equations for a homogene- 
ous half-space free of sources can be constructed of par- 
tial solutions with the components? 


Jip ee. cos (md + am) 

By = ee sin (md + am) 
(4) 

E, = Em cos (md + am) 


* Received by the PGAP, June 16, 1960. 

t+ USASRDL, Fort Monmouth, N. J. 

1S, Silver, “Microwave Antenna Theory and Design,” M.I.T. 
Rad. Lab. Ser., McGraw-Hill Book Co., Inc., New York, N. Y., 
p. 173; 1941. 

2J. A. Stratton, “Electromagnetic Theory,” McGraw-Hill Book 
Co., Inc., New York, N. Y., pp. 360-371; 1951. 


5 ae fe ‘om Sin (mp + am) 
Hs = + \/~ saa cos (mp + ain). 
vy 
ede Va om sin (md + am), (1a) 
with 
fo =f. $2 {(k £ h) Imaal) 
+ (k = h)Im—1(yp) § eo ™*ydy 
joa = =f. fn (V){ (B + h)Im+i(ve) 
— (k F h)Im—s(yp)} e~*ydy 
Em = ¥j f fn (Y)Imlyp)ey dy (1b) 
0 


Fig. 1. 


fn (y), fm (y) are arbitrary functions of the real 
variable y, subject to the condition that the integrals 
in (1b) and the integrals over the first derivatives of the 
integrants with regard to p and g exist. fm(y), Fin (Y) 
represent the amplitude spectra of the elementary 
cylindrical waves which form the beam. The propaga- 
tion constant h is determined by the relation R? =y+h’, 
where k= 27/2 is the plane wave propagation constant 
in the half space. Jm1, Jm and Jmn41 are the Bessel func- 
tions of the order m—1, m and m+1, respectively. m is 
any positive integer including zero. 

For symmetry reasons, a field representation in terms 
of hybrid waves has been chosen in preference to the 
more conventional one in E-waves and H-waves. There 
are two alternatives for combining E- and H-waves of 
the same order to hybrid waves symmetrical in E and 
H. They are distinguished in (1a) and (1b) by the super- 
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oS to) 


scripts (+) or (—). Pure E- or H-wave fields are ob- 
tained by superposition of two hybrid fields with f,, (y) 
= Ff (7). 

The elementary waves are propagating waves in the 


range 0<y<k and evanescent waves in the range 
Y<k< om. 


Assume the beam contains only propagating waves 


and fn? (y) and fm(y) are real functions. Then, as 
_ evidenced by (1a) and (1b), the transverse field com- 
ponents at z=+29 and z= —2 are conjugate complex. 


Assume, furthermore, that there is a phase transformer 
which, when inserted into the beam at z=20, transforms 


the phase distribution at this plane into a phase dis- 


tribution identical with that at z= —2zo. Then, the orig- 
inal beam is reconstituted; this means, the field dis- 
tribution in the space range —20<z<-+29 is repeated 
in the range +29<2<+3zo. If several of these phase 


transformers are placed at regular intervals of 220, a 
_ reiterative beam is obtained, which is “guided” by this 


phase shifting structure. 


The question as to the physical realizability of the 
phase transformers depends on the structure of the 
beam. Fortunately, there are beams where the required 


_ phase transformation is obtainable with very simple 


means, at least within a finite cross-sectional area. These 
beams are characterized by a y-spectrum substantially 


- limited to values of y?k?. In other words, the functions 


fm (y) have appreciable values only within a range 
0<y<7 where y.2«k?. In this case, the propagation 


| constant h in the phase terms of (1b) can be approxi- 


mated by h~k—vy?/2k, provided (y.?/2k)z0 remains 
within the order of 27. This imposes a limitation of the 
go-range in which the approximation is valid. In the 


amplitude terms of the transverse field components, 


(1b), # can be replaced by k. Thus, (1b) simplifies to 


(+) —jke [°~ (+) +5 (7?/2k)2 
BO = Be abe [fe Imire yay 
0 
he en a 43eP[one Y" 
Emm = + jke fin (VY) Im(vp)e ate (2) 
0 


Since the field is uniquely determined by the trans- 
verse field components, only these components are con- 
sidered in the following. 

Due to the restriction on fn‘*)(y), the transverse elec- 
tric as well as the magnetic field components have the 
same cross-sectional distribution. The phase trans- 
formation which is required for reconstituting or iterat- 
ing such a beam is indeed feasible, at least theoretically, 
as shown in Appendix I. 

Forthcoming deductions prove it very advantageous 
to represent fn (y) in (2) by the following series of 

£2) = Dionmfma (Y) = 


orthogonal functions: 
No 
An,m\~_ 
n=0. n=0 Yo 


2 
yi 2). 2 
“Ln ma(75) ee ees AS) 


Yo 


(+) 
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where Yo is an arbitrary constant. 


cena L | 


AN ds Se 


(—x)* 


4! 


are the generalized Laguerre polynomials,’ and 


n} ” ) YP ae 
=2—__ [4 () 
NE eee 


2 
sy bee (=) e271) (=) 
Yo” Yo 


are the coefficients of expansion. The Laguerre poly- 
nomials L,,1(x) form for each / a complete system of 
orthogonal functions for the range 0<x< with the 
weight function «'e~*. 

The fields derived from fo (y) and fo (7) are linearly 
dependent, since the field components, ((1a); (1b) 1, 
with the superscripts (+) and (—), transform into each 
other if a is replaced by t—a@p. Therefore, only one of 
the functions, namely fo, is considered in the follow- 
ing. 

Now, (2) can be evaluated, and we obtain the field 
distribution functions* 


cE) 


An ,m 


(4) 


1 


(+) 


=e Hen 


+ gi2n aretg Leo"/k) 1 ym™t1T may 


«ea 1/2 (079?) | a (yok) 2) | 


The beam can thus be considered as the superposition 
of elementary beams having the field components 


(+) gkz 


(+) jke 
| i a 


zi Fo sm ne 


(+) ymett 


AE u 
m,n Ina ( ) 
(1 + gp?) 12@nE1)) 1+ 2? 


{ 1 u? 
ale —_—— 
i 2 1+ 


=A 


+j|Ontm£1+4, 


1 uv | (5a) 
GBC Re ame : a 
. 21+? 
(+) Bee) es € (4) 
oe = + EL one =-+ 4/ Toa 
M 


(Sb) 


l| 
| 
ae 


3 W. Magnus and F. Oberhettinger, “Functions of Mathematical 
Physics,” Chelsea Publishing Co., New York, N. Y., pp. 84-86; 
1955. 

4A. Erdely, “Tables of Integral Transforms,” McGraw-Hill Book, 
Co., Inc., New York, N. Y., vol. 2, p. 43; 1954. 
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with 


Y Yo" yo\? 
w= e=(%) (ip), v= = (7) ce), 


(+) 
WW = 


ie ee 
aL) eyo: 
For m=0, only the upper sign is valid. 
The elementary beams or “beam modes” satisfy the 
same orthogonality relations as the modes in an ordi- 


nary loss-less waveguide 


f (ee SHO \edF 
F 


fee * n+m + 1)! 
E \/~ A eels e \ ) as Own! 
Ml Yo n! 
ib CEs, x Hy esd = 0, (6a) 
F 
* * wr (n+tm dt i)! 
eee 
F Yo" n! 
Gye) ’ 1ln=n' 
Enka war = 0 with Orn = (6b) 
F Onn’, 


where the integration is performed over any plane z= 
constant; e, is the unit vector of the z-direction; Em,» 
and H,,,,“*) denote the transverse field vectors of the 
beam modes and E,,,,*, Hmn»‘t)* the conjugate com- 
plex vectors. The proof of (6) is given in Appendix II. 

Eq. (6a) shows that the total power of a beam com- 
prising several beam modes is the sum of the powers of 
the individual beam modes. The power of a mode is 
given by 


ee eee 
mM 


2 


Yo n! 
for m ~ 0 
fee: +r (n+1)! 
po — sf Ave ase get uae 
2 mn Yo n!} 
form=0. (7) 


The vector functions Ep. and Ea form a com- 
plete system of orthogonal functions. With reference to 
(1a), we note that each vector function represents two 
modes characterized by an=0 and an=m/2, respec- 
tively. Any given field distribution in a plane z-constant 
can mathematically be developed into the mode func- 
tions. However, this development has physical meaning 
only if the corresponding amplitude distribution func- 
tions fm‘ (y) satisfy the condition introduced at the be- 
ginning of this paragraph; viz, that fn*)(y) has appreci- 
able values only within a range 0<y<y, where VWeKRk. 
Similarly, any individual beam mode has only physical 
reality if its amplitude distribution function 


(+) Woven 2 
tee y) = (=) Ins (2 e1/2y7/y0") (8) 
Yo 70" 
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decays sufficiently rapidly in y. Otherwise, the field 
associated with the mode would not be consistent with 
Maxwell’s equations. Thus, there is an upper bound 
for Yo/k, which depends on the order n, m of the 
Laguerre polynomial. 

Each beam mode, as defined in (5), has a plane phase 
surface at z=0. The phase distribution in a plane 240 
is given by the factor 


uy 


1 
exp {i | (an +m+1+ 1)arctgv — fe + =|} 


= exp {j[Jn.0(2) + ¥(0,2]}, 


where 
Drs) Yo" 
Vmin (2) = (20 + mt 1 + pare (= :) (9a) 


is independent of p and 


Q 1 
V(p, Z) See 


2 yo \? 

a 
is independent of the order n, m of the Laguerre poly- 
nomials. If a phase transformation is performed at the 
plane z= +20 by a factor e~#¥ 0) the field distribution 
at = —2Z of each beam mode is repeated at z= +2, 
disregarding the constant phase shift 2m,‘ (go). Thus, 
any elementary beam can be restored at intervals 229 


by phase correcting means which effect a radially de- 
pendent phase shift 


¥(p) = — 2p, 0) + vo 


(9b) 


as 9 & yop” - Yo, 
14 C 3) 
k 


where wo is an arbitrary constant. 

Although every single beam mode is iterated by the 
same phase transformation, (10), wave beams composed 
of several modes are not iterated because the independ- 
ent phase shift (9a) is different for every mode. A sys- 
tem of equally spaced identical phase transformers guide 
such a beam, 7.e., prevent the energy from spreading 
out, but the field distribution is not repeated from inter- 
val to interval. The conditions are similar to those in an 
ordinary waveguide where a field consisting of several 
modes is guided with continuously varying cross-sec- 
tional field distribution. 

The concentration of energy in a beam mode can be 
characterized by the fraction 7m, of the total energy 
which passes through a plane z=constant within a circle 
of radius p=p,. If a certain ymn is considered, p, varies 
with z. The p,(z) which characterizes the envelope of 


(10) 


May i 


| 
| 
| 
| 
| 
di 
| 


‘ 


1961 


constant energy flux can be derived from (5). Substitut- 
ing x=u?/1+v?, one obtains 


(4) n! fee £1 7 ( y}e xq 11 
UT Tas x nim x “ 
( 1)! A ymt1 é x ( ) 


0" Pn 


70 : 
ae ee :) 
k 
This equation shows that the dependency of m,n‘ on 
p, and z is only contained in xo. The relation between 


p, and 2 for 7m,n‘* =constant, 7. e., x»=constant, can be 
normalized to 


xo = 


(12) 


with 


p,/Roas a function of z/z9 with vp as parameter is plotted 
in Fig. 2. If the distance 22) between two successive 
phase transformations and the parameter yo character- 
izing the set of beam modes are given, all the curves 
p,(z) are similar and are represented by one of the 
curves of Fig. 2. In other words, assuming the same 2p, 
each curve describes the energy flux distribution for the 
entire mode system with the parameter Yo. 

The question can be asked: What value of yo yields 
the highest energy concentration, 7.e., the smallest “field 


Fig. 2—Envelope curves of constant energy flux 
in normalized coordinates. 
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diameter” at the locations of the phase transformers, 
if their mutual spacing is given? The desired yo value is 
obtained with (12) by introducing the requirement that 
p,/Ro, as a function of vo, is a minimum for =z. The 
requirement is satisfied, if vo = (yo?/k)20=1. 

This condition for maximum energy concentration is 
independent of ” and m and, therefore, is the same for 
each beam mode. The corresponding phase transforma- 
tion is 

k 
¥(o) = + — p* + Yo. (13) 
2Z0 

The case of maximum energy concentration is illus- 
trated in Fig. 3 which shows various curves m1,9~ 
=constant for the beam mode m=1, n=0, (—). Asa 
numerical example, we consider a beam at the frequency 
of 35,000 Mc, whose phase is adjusted at intervals 
of 10 m (g.=5m). The diameter of the beam for 90 per 
cent energy flux varies from 25 cm (g=0) to 35.4 cm 
(g=20), for 99 per cent from 35.6 cm to 50.4 cm, and for 
99.9 per cent from 43.6 cm to 61.6 cm. 


JUSTIFICATION OF THE SIMPLIFICATION UNDERLYING 
THE THEORY 


The simplified field equations which led to the beam 
mode system had been obtained with the assumption 
that the amplitude spectrum fn‘*)(y) has considerable 
values only for y<k. Since the amplitude function for a 
beam mode as introduced in the preceding section is 


m+1 2 
fan) = (=) J Bp caraash (=) e127 170") | (14) 


Yo Yo 
40 
- 
VK p 
@ 
ie es 3.6) te 
pS a 99 i 32 = oI) 
28 092 
| 325 ls 
oo 
22> 
24 FQ % A 
oe 
20 
ie cs 16 1b | 
Ah A 12 
ie IIe 
04 if Hea 
-10 -08 -06 -04 -02 0 02 04 06 08 10 
ee. 
Zo 


Fig. 3—Envelope curves of constant energy flux for the beam mode 
n=0, m=1, (—) in the case of maximum energy concentration 
at z= +20. 
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where Ln m+i(y?/Yo2) is a polynomial in y?/yo’, this as- 
sumption requires that 7» is sufficiently small. The field 
distribution of the mode in the plane z=0, where the field 
has the highest concentration, is, according to (5), 


(+) (+) 
es = gm ,n 


= (—1)"by0?(vop)™*1Ln ms1(yo2p2)e-*/200"#, (15) 


This equation shows that the condition of small yo im- 
plies large beam diameters. 

In order to get an estimate as to the meaning of “large 
beam diameters,” consider a beam mode with the param- 
eters m=1, n=0, (—), where 


—1/2 (72/702) © —1/2(yop) 


= (>) 
fio) ee E,i,o = — Egi,o = Rye 
atz=0. (16) 


If 99 per cent of the energy pass the z=0 plane within 

a circle of radius p=5A, the corresponding value of 

Y0/k, as obtained with (11), is yo/R=6.8X10-?. Thus, 

(-) —110(7"/k”) 

Tian e * 

At y/k=0.25, fi,o(y) is already as low as 107%. If 

fi,o (y) is assumed to be zero for y/k>0.25, the rela- 
tive error of the energy integrals , 


(-) ;2 (-) ;2 
ef | ei l'ar =a f | Hee lar 
F F 
is 0.23 per cent. The phase error introduced into the 
phase term when / is approximated by k—vy?/2k is 
below 0.12 radian at z=. These errors decrease rapidly 


if the beam radius is increased. Higher modes require 
larger radii for the same magnitude in errors. 


EFFECT OF Cross-SECTIONAL LIMITATIONS 


In the preceding discussion, the assumption had been 
made that the phase transformation is performed over 
the entire plane z= 2 and corresponding planes at inter- 
vals of 2zo. 

Actually, the phase transformation can only be 
achieved over a finite area. A practical solution for the 
phase transformation is the use of dielectric lenses, as 
will be explained in another paper.® 

Although the field of the beam modes, particularly the 
lowest one, decreases rapidly in the radial direction, 
some distortion will be introduced if the phase is cor- 
rected only within a finite area. The energy outside 
this area can be considered lost, since the uncorrected 
field diverges and is radiated into space. The beam mode 
system which has been derived with the assumption of 
unlimited phase correction will be modified in two re- 
spects: 1) The field functions will depend on the dimen- 
sions of the phase correcting devices, and 2) there will 
be an attenuation factor which is caused by the loss 
of energy bypassing the phase transformers. This at- 


5 J. R. Christian and G. Goubau, “Experimental studies on a beam 
waveguide for millimeter waves,” this issue, pp. 256-263. 


(17) 


tenuation factor does not include any absorption or 
reflection losses which may occur in actual phase cor- 
recting structures. 

In the following, the assumption is made that the 
phase correction at the planes z=(+21+1)z0, where 
1=0, 1, 2,-+-, is performed within a circle of radius 
R, while the area outside this circle is covered with an 
ideal absorber which prevents the field from being re- 
flected or transmitted into the succeeding space section. 
The modes of this system are characterized.by the con- 
dition that the field distribution in successive phase cor- 
rection planes z=(+2i—1)z) and z=(+2i+1)z is the 
same, except for a complex constant whose absolute 
value determines the loss caused by the finite area of 
phase correction. 

The field at z= —2o, after being subjected to the phase 
transformation, is according to (2) 


(+) GE) (+) 
Fam = + Egm = Em (p, — £0) 


—j (7 /2k)z9 


fee ees 
ee } fo aT Gee vay; (13) 
0 


where fn‘+)(y) may be complex. E,,%(p, —zgo) can be 
considered as the Hankel transform of fm‘) (y)ei— ("2% 
of the order m+1.° Therefore, if Em‘*)(p, —zo) is given 
fm‘*(y) is obtained by the inverse transformation 


2 R 
(+) jkz9(1—y / 2k) (4) 
ee i) Em (p; —20) Jmsx(yp)pdp. (19) 
0 


ee 
Fm ewes 


The integration is only extended from zero to R instead 
of infinity, because EZ,,‘*)(p, —zo) is limited to this range 
by the absorbing screen. 

The field at z=+29 before being subjected to the 
phase correction is obtained by inserting fm‘*)(y) of (19) 
sm 


-} R 

(+) —j2kz (+) 

Be tae if f Eo Goan 
y=0 n=0 


‘Imer(yp)eXY!*mdnydy. (20) 


The integration over y can be performed 
if Ime ilyp) J ma(y) ef 1k) 0 dry 
0 


k k 
Se peed fe Sy Co en 3 (k 1420) (179) | a4. 
J] 220 (= m) ve 


Thus, 


(4) Pub ae Br aes k 
220 0 220 


e771 (k1 429) (1+ adn. (22) 


° W. Magnus and F. Oberhettinger, of. cit., pp. 136-137. 


ee a ae 
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This field, after being subjected to the phase correc- 
tion, becomes zero in the range p> R, which is covered 
by the absorbing screen. The field is somewhat modified 
within the range 0<p<R because of the boundary con- 
ditions. However, since we consider beams of large diam- 
eters, the modification is limited to a narrow zone along 
the rim p=R and may be neglected. Hence, the field 
behind the phase transformer is obtained by applying 
the phase transformation to the center portion 0<p<R 
of the incident field (22). 

Obviously the effects of the limitation of the beam 
will be minimized if the phase correction is adapted to 
maximum field concentration at z=, as given in (13). 
Then, the field at =z, after the phase has been trans- 
formed, is 


(+) 5 ((k/22z9) p40) 
Jars (p, +zo)e 


This field shall be the same as the field E,,+)(0, —20), 
except for a constant complex factor ». Thus, we can 
write 


Oecip <k (23) 


(+) 7 ((&/220) p+) (+) 


Em (p, +20)e = pEm (p, —20) - 
forO<p<R. (24) 
Inserted into (20) we obtain 
"yee (p, ee ee 
k 4 (2ke R 6 Byes 
oy easy 4 (2k sal eo, MY i (k/420) 9 
220 0 


k 
*Imt1 (— no) ndn. (25) 
220 


We introduce the following abbreviations: 


(t) —i(k/Azoe® Ct) JL —yekeotv0) 
Em (p, —20)é lim (X59 ae =4 


ae Rk NRE 
— = —n, a= — R. 
2Z0 aie Mes 220 
Then (25) transforms into the following normalized ho- 
mogeneous integral equation of the second kind? 


+= (26) 


F(a) = 4 f "PP @Inalabede, (27) 


whose eigenfunctions and eigenvalues are real and only 
depend on the quantity a.° 

If the phase corrections at = (+ 27+1)z» are thought 
to be performed in two equal steps, then, Fn‘ (x) repre- 
sents the field distribution in between these two steps. 


7 This equation can also be derived with Fresnel-Kircheoff’s dif- 
aa eae This approach has been used by G. Goubau and 
R. Christian, “A New Waveguide for Millimeter Waves,” presented 
at the URSI-IRE Fall Meeting, San Diego, Calif. ; October 19, 1959. 

®R. Courant and D. Hilbert, “Methoden die Mathematik 
Physik,” Springer-Verlag, Berlin, Ger., vol. 1, pp. 96-138; 1937. 
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Since F,,)(«) is real, the field in the middle planes of 
the phase transformers has uniform phase. 

The eigenfunctions F,,~(«) and Fm‘(x) form each 
a complete orthogonal system. Together they determine 
an orthogonal system for the transverse field vector. 
They characterize the modes of the system with cross- 
sectional beam limitation. The eigenvalues g are greater 
than one. 1/| q| = | 2| represents the attenuation of the 
field between two successive phase corrections. Any 
beam of arbitrary cross-sectional field distribution can 
be considered as the superposition of these modes. If. 
such a beam is passed through a large number of phase 
transformers, the “higher modes” which correspond to 
the higher eigenvalues are gradually damped out until 
only the lowest mode remains. 

If a approaches infinity, the modes transform into the 
unlimited beam modes. The eigenvalues become +1, 
since there is no beam limitation and, therefore, no 
diffraction loss. 

The following discussion is restricted to the mode 
group having the parameters m=1, (—). These modes 
are polarized in one direction. They are the only modes 
whose transverse field is not zero at the axis. 

In order to solve the integral (27) one can proceed 
in the following manner: 

The kernel Jo(x€) (for m=1, (—)) can be expanded 
into a series of the mode functions of the unrestricted 
system 


Tox) = 230 (1) La(22)-LalEDeUEY, (28) 


n=0 
where 


Ibs — Loa. 0 


The series converges for any x and &, except zero and 
infinity. At x =£=0, the series alternates between zero 
and two. 

An approximate solution of the problem can be ob- 
tained if Jo(xt) is replaced by the first NV terms of the 
series. The resulting approximated integral equation can 
be solved rigorously. Expanding Fi (x) = F(x) also into 
the mode functions 


N 
F(x) = Qo dnLn(a*e-O", 


(29) 
n=0 
the approximated integral equation yields 
N 
a On Ln (x2) e~ (1/2) 
n=0 
N ON ; 
= 4 ey > (—1)"BurdpL,(a2)e- 1/22", (30) 
r=0 n=0 
where 


Rael fi Lala?) L,(2)eF'nd, 
0 
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Fig. 4—The first eigenvalue of integral (23). Curves I, II and III 
are obtained by expanding the kernel Jo(xé) into the system of 
mode functions using 1, 3 and 5 terms, respectively. Curve IV 
has been computed with an IBM computer applying power series 
development. 


(30) resolves into a system of linear homogeneous equa- 
tions for the coefficients a,, since the equation must be 
satisfied for any x 


N 
De (De = 1) dbnn) on = 0 0 < r = N, 
n=0 


1lr=n 
eee. (31) 


nr 


The determinant of this system must be zero. This con- 
dition leads to an algebraic equation of Nth degree. The 
roots are the desired approximations for the eigenvalues. 
The corresponding coefficients a, determine the ap- 
proximated eigenfunctions. 

Fig. 4 shows curves for the first eigenvalue as a func- 
tion of a=~W/(k/2%)R calculated in this manner for 
N=1, 3 and 5 (curves I to III). Curve IV has been ob- 
tained with an IBM computer using for Jo(x&) a power 
series expansion. The ordinate is given in a db scale and 
refers to p? which determines the diffraction loss caused 
by the finite apperture. Fig. 5 is a plot of eigenfunctions 
for a=1.8, 2.0, 2.2 and 2.4, as obtained with the com- 
puter. 

In order to obtain an estimate for the attenuation of 
the higher modes, the eigenvalues for n=2 and 3 have 
been calculated with N=3. The results are given in 
Table I. Presumably, these values are not very accurate, 
since higher modes require larger N for the same degree 
of accuracy. 


Fig. 5—The first eigenfunction of integral (23) 
for various values of a. 


TABLE I 

a p(n=1) b(n =2) p(n =3) 
1.4 0.797 —0.040 0.102 
1.6 0.911 —0.142 O.113 
1.8 0.968 —0.324 0.158 
Ze0 0.990 —0.538 0.236 
De 0.997 —0.727 9.348 
2.4 0.999 —0.859 0.487 


0 2m 4m 6m 8m 10m 12m 4m 


Fig. 6—Calculated diffraction loss per km of a beam waveguide at 
35,000 Me for various beam diameters 2R as a function of the 
transformer spacing D = 22». 


The diffraction loss of a beam at a frequency of 35,000 
Mc, guided by an array of phase transformers over a 


- distance of 1 km is shown in Fig. 6. Each curve refers 


to a certain diameter of the beam measured at the loca- 
tions of phase transformation. The loss is very small as 
long as a2 2.4 and increases rapidly with decreasing a, 
z.e., with increasing z9. The spacing increases for the 
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same amount of loss more than with the square of the 
beam diameter. F 


APPENDIX I 
Ideal Phase Transformer 


Assume at = 29 a thin layer (thickness 6) of dielectric 
material whose € and p are functions of p and ¢, but not 
of z. A wave beam of the kind described by (12) shall 
pass through the layer without producing a reflected 
field and, thereby, change its phase distribution in a 
predescribed manner. 

The transverse field vectors E; and H; at both sides 
of the layer are 


i 
Ey = E(p, ) /~H = ¢, X. E(p, >) 
€ 
for z = 20, 
E; = E(p, pein), \/~ H; = ez X E(p, d)e* 
(= 


forz = 2+ 6. (32) 


a(p, @) is the desired phase shift which is a function of 
p and ¢. e, is the unit vector in the z-direction. 

Since the transverse field components are continuous 
at the interface of any two media, E; and H; are the 
same inside the layer. The longitudinal field components 
of the incident field are reduced by the factors €/e and 
u/ui, respectively, after the field penetrates the inter- 
face. If e, and mw: are sufficiently large, the longitudinal 
field components within the layer are negligible. Thus, 
in spite of the variation of the field with p and ¢, the 
field inside the layer can be treated as an ordinary plane 
wave. The condition for no reflection at the interfaces 
is €/u=6/m. The required phase shift inside the layer 
is a(p, 6) = —w/m1e5+2rn. The distribution functions 
of « and sz, which are required to effect this phase 
shift are, therefore, given by 


1 oe 
i ait = =4/ ease ey - G3) 
w a 


1 
ui(p, 6) = —A/ — (2xn — a(p, $)). (34) 
63) (2 
Since the phase shift depends on the products «4-6 
and y4-6, the thickness 6 can be made arbitrarily small 
if « and p; are sufficiently large. Instead of varying 4 
and ju, the thickness 6 may depend on p and ¢. The 
integer n may have different values for different ranges 


of p. 
APPENDIX II 
Benos of the Orthogonality of Beam Modes 
- We first normalize the beam modes by setting 


(4) n! a ae : 
ie ee ae 12, (35) 
Ms ee m + If 
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First (6b) shall be proven. The integrals in these equa- 
tions are extended over a plane z-constant. Thus, 


f BORO oar 
F 


0 Qa 
even Che (+) (+)* 
= f f { | DRAG DES, yn! + PF pete ‘ een wn? } dodpdp 
p=0% ¢=0 


(¥)* 


(+) 
f a4 . En wah 
F 
Pi (END aan ae 58 (2) 
= f if {Ee ello mn! =e Emin Em’ n’} dppdp. 
p=0 o=0 


Eymn® and Egmn™ are replaced by the right hand side 
of (5a), substituting 


u? Yo Pp 
x= = 
1+ v 


we obtain with the normalization (35) 


(+) (+)* 
f Enn' En wF =e 
ed 


j{2(n—n')+ (m—m')} arctg v 


n! n' | \ 1/2 
; ee is 1! +m! a 1)! 
f gilemtm ENT naa (%) Ln? m'1(%) € 7dx 
2=0 


1 Qa 
_— cos {(m — m')o + (Qn — 


21 J 9 
(+) 
it Ena: 


n! n'! ee 
ie =e 1)! (n! +m’ ¥ 1)! 


Om’) } do (36) 


(F)* j{2(n-+n’)+ (m4+m’)+242} arctg v 


Em’ nF ne 


2 
, {, ah l2(mtm!) Tm (%) Ln’ sm’ F1(K) eax 
z=0 


—{- cos {(m + m’')o + (a + Om’) | dd. (37) 


Eq. (36) applies to all values of m and m’. Eq. (37) 
does not exist for m=m'=0, since it applies to modes 
with different superscripts (+) and (—), while in the 
case m=O, there are only modes with one superscript 
which has been chosen to be (+). 
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The integrations with regard to ¢ in (36) and (37) 
are readily performed 


1 Qa 
—{ cos {(m — m')¢ + (am — Om) } db = 8mm (38) 
217 Jo 


Qn 


1 
— cos { (m + m')o + (am + am’) | dd == 0) 
2a 0 


(39) 
The orthogonality relations of the Laguerre poly- 
nomials yield 


n,n! . 


(n+ k)! ah 
n!} 


i * Dn n(X) Ln (ae "dx = 
0 


Eq. (39) inserted into (37), and (38) and (40) inserted 
into (36) prove the orthogonality relations (6b). 
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The proof of (6a) is now readily obtained, since with 
(5b) the relations 


Be Grr G5. 
Go. x Heo = (/~ En’ En’ yn! 
(ey Se) 
[E2x HP = 4/4 eee CAD 
€ 


can be easily verified. Thus, (6a) is ms themed 
identical with (6b). 
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Experimental Studies on a Beam Waveguide for Millimeter Waves* 


J. R. CHRISTIAN}, mMemBER, IRE, AND G. GOUBAU{, FELLOW, IRE 


Summary—The beam waveguide utilizes reiterative wave beams 
which are guided by resetting the cross-sectional phase distributions 
at periodic intervals. The major purpose of the measurements was to 
obtain data concerning the inherent losses of the guide, including the 
diffraction loss, caused by energy by-passing the phase transformers, 
and the dielectric and reflection losses of the transformers. Reso- 
nance measurements with an “open cavity,” consisting of a section of 
a beam waveguide terminated by plane reflector walls, yielded the 
over-all losses of the phase transformers. In order to isolate the dif- 
fraction loss, pulse measurements were made between two reflectors, 
slightly curved to reset the phase distribution in the reflected beam. 
These measurements also gave information on the build-up of the 
reiterative wave beam. The experimental data support the theory 
which predicts that the diffraction losses of reiterative beams can be 
made extremely small. 


I. INTRODUCTION 


HE waveguide discussed in this paper differs from 
“ae waveguides in that the energy is not trans- 

mitted in the form of a wave mode. A wave mode 
is characterized by a certain propagation constant and 
a certain field pattern which is everywhere the same 
along the guide. The field in a beam waveguide consists 
of a wave beam, that is, a bundle of waves characterized 
by a spectrum of propagation constants. The field dis- 
tribution varies along the guide, but is periodically re- 
peated at certain intervals. The guidance is accom- 
plished by passing the beam through “phase transform- 


* Received by the PGAP, June 16, 1960. 
1+ U.S. Army Signal Res. and Dev. Lab., Fort Monmouth, N. J. 


rs,” that is, phase shifting devices which reset the cross- 
sectional phase distribution in the beam at predeter- 
mined intervals. 

Beams which can be guided in this manner are called 
reiterative wave beams. The theory of reiterative beams 
is contained in a preceding paper.! There are sets of 
reiterative beams called beam “modes” which satisfy 
the same orthogonality relation as the modes in an 
ordinary waveguide. Beam modes of the same set re- 
quire the same phase transformation for their guidance. 
The attenuation differs for the various modes. The low- 
est mode whose radial intensity distribution follows es- 
sentially an e-@" curve, is of particular interest since 
it has the lowest attenuation for a given beam diameter; 
it is this mode which is utilized in the beam waveguide. 

Fig. 1 shows a schematic drawing of a beam wave- 
guide. It comprises uniformly spaced phase transform- 
ers which represent the actual guiding structure and 
terminations for launching and receiving the desired 
beam mode. The phase transformers are designed to ad- 
vance the phase in the outer portion of the beam relative 
to the center according to the relation 


Rp? 


1G. Goubau and F. Schwering, “On the guided propagation of 
electromagnetic wave beams,” this issue, pp. 248-256. 
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where @¢ is the phase advance at the radius p, D, the 
spacing of the phase transformers, and k=27/\. The 
maximum phase shift @max=RR?/D (R is the radius of 
the phase transformer) is directly related to the diffrac- 
tion loss Lp of a phase transformer, that is, the loss 
caused by the cross-sectional limitation of the beam. The 
relationship between ¢@max and Lp is shown in Fig. 2 
which is essentially curve IV of Fig. 4 in the above 
quoted paper.! The loss drops below 1/100 db if @maz ap- 
proaches 27. Thus the maximum phase shift for all 
practical cases is approximately 27. 

Phase shifts which obey a quadratic relation in p can 
be accomplished by use of dielectric lenses. However, it 
shall be emphasized that the lenses in this application do 
not function like optical lenses. Geometric optical con- 
ditions would require much larger values of max, 1.€., 
much higher operational frequencies. The focal length 
of the lenses which follows from (1) is D/2. Under geo- 


PHASE TRANSFORMERS 


HORN 


LAUNCHER 


Fig. 1—Schematic of a beam waveguide. 


DIFFRACTION LOSS—» db 


Fig. 2—Calculated diffraction loss as a function of 
maximum phase shift of a phase transformer. 
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metric optical conditions each lens would image the 
preceding lens into the succeeding lens, 


Il. MopaL BEAM WAVEGUIDE 


In order to demonstrate the existence of reiterative 
wave beams and to obtain some information on the 
performance of a beam waveguide, a simple model was 
built for frequencies around 23 kMc. A photograph of 
this model is shown in Fig. 3. The phase transformers 
were dielectric lenses set in wooden frames. Their diam- 
eter was 20 cm and their mutual spacing approximately 
one meter. Lauching and receiving were accomplished 
by means of horns in combination with lenses. The aper- 
ture of the horns was circular and had a diameter of 10 
cm. Standard K-band rectangular waveguide was used 
to feed the horns. The wave beam emitted from the 
horns was substantially linearly polarized. 

Two types of phase transformers were used in the 


Fig. 3—Photograph of a model beam waveguide. 


model, both employing different methods for minimiz- 
ing reflection, since reflection causes transmission losses 
and undesired standing waves. 

One type (Fig. 4) employed a single lens of polysty- 
rene foam with a dielectric constant of only 1.2. Because 
of the low dielectric constant the reflection loss is less 
than 0.02 db. The dimensions are indicated in Fig. 4. 

The other type of phase transformers comprised a pair 
of identical lenses mounted back to back, as shown in 
Fig. 5. The design of the lenses was based on the well- 
known fact that the phase of a plane wave can be 
shifted without reflection by passing the wave through 
a pair of appropriately spaced parallel dielectric plates 


of equal thickness. Each pair of adjacent elements, indi- 


cated in Fig. 5 by 1 and 2, has been treated like a pair 
of parallel dielectric plates effecting the desired phase 
transformation at the location of these elements. Al- 
though such a design procedure is only justified for 
phase transformers of extremely large diameters, it is 
apparently applicable to diameters of only 16 wave- 
lengths, since the phase transformer exhibited negligible 
reflection at their design frequency of 23 kMc. 
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Fig. 5—Double-lens phase transformer. 
Radii of curvature, Ri = 75 cm, R2 = 216 cm, ¢ = 2.55. 
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Fig. 6—Schematic of field distribution measurement setup. 
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A number of field measurements was made to explore 
the cross-sectional field distribution within the space 
between the phase transformers. Fig. 6 shows a sche- 
matic of the measuring setup. An open-ended flanged 
waveguide was used as the probe for most measure- 
ments. A rectangular horn, 2 cm across its aperture, was 
used to investigate the lower level field amplitudes 
(more than 25 db below Emax). The indicator consisted 
of a calibrated crystal detector connected to a milli- 
voltmeter. Sa 

The field distribution in the first sections of the beam 
waveguide depended largely on the adjustment of the 
launcher, that is, on the mutual distances between the 
horn and the lenses. This adjustment affects the field 
distribution of the beam developed by the launcher. If 
this distribution deviates from that of the lowest beam 
mode which is propagated by the lens structure, higher 
beam modes are excited simultaneously. In spite of their 
greater attenuation, the presence of these modes is still 
noticeable within the first sections of the beam wave- 
guide. It was found that even after the beam had passed 
eight transformers, the field distribution was still some- 
what dependent on the adjustment of the launcher. 

In order to obtain representative field measurements 
of the reiterative wave beam with a small number of 
available phase transformers (8 and 10 respectively), 
the launcher was carefully adjusted to minimize the 
development of higher modes. A criterion for this 
adjustment was that the field distribution remained sub- 
stantially unchanged after the beam had passed seven 
phase transformers. 

Fig. 7(a)—7(d) shows measured cross-sectional dis- 
tributions at four different distances (D/4, D/2, 3D/4, 
D) behind the eighth phase transformer. The foam 
lenses were used in these measurements but the double 
lens transformers yielded similar results. The abscissa 
indicates the distance from the beam axis, and the or- 
dinate indicates the ratio [E(p, z)/E(0, D)]? where 
E(p, 2) is the field at the distance p from the axis and 
E(0, D) the field in the center of the beam at the dis- 
tance D. The broken line refers to the E plane and the 
dotted curve to the H plane. According to the theory, 
both curves should be identical if the beam diameter is 
sufficiently large. The deviation between the E- and 
H-plane curves may have been caused by the fact that 
the beam diameter of 16 wavelengths was not large 
enough and/or that the launching horn which radiates 
differently within the E and H planes causes the excita- 
tion of an azimuthally periodic beam mode of also rela- 
tively low attenuation. 

Fig. 7(d) which refers to the field distribution at the 
distance z=D, that is, the location of a phase trans- 
former, has been measured for a wider cross-sectional 
range, thus showing the distribution of the field which 
bypasses the phase transformer. The solid curve is the 
computed field distribution for Vk/D-R=2.2. This 
curve has been taken from Fig. 5 of Goubau and 
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— Joke Ke 
(c) 


(d) 
Fig. 7—(a) Cross-sectional field distribution at a distance z=D /4 behind the eighth phase transformer. (b)—Cross-sectional field distribution 


behind the eighth phase transformer. (c)—Cross-sectional field distribution at a distance z=3D/4 behind the eighth 


phase transformer. (d)—Cross-sectional field distribution at a distance =D behind the eighth phase transformer. 


Schwering.! The actual value of /k/D-R for the phase 
transformers used in the model was somewhat smaller, 
namely 2.15. 

In order to obtain information on the efficiency of the 
beam waveguide, the over-all transmission loss was 
measured between the input terminal of the launching 
horn and the output terminal of the receiving horn. The 
measured loss was 2.5 db using 8 foam lens transformers 
and 1.9 db using 10 double lens transformers. In both 
cases the spacing D was 1.1 meter. The loss remained un- 
changed when the wooden frames were covered by 
aluminum. This observation is of particular interest 
with regard to the theory which ignores specific bound- 
ary conditions at the rim of the phase transformers." 


The transmission loss consists of two components, the 
“iteration loss” which is caused by the phase transform- 
ers and the “launching loss” caused by the terminations. 
Since the launching loss is the same in both cases, the 
iteration loss of the foam lenses is apparently larger 
than that of the double lens configurations. As will be 
shown later, the iteration loss of a foam lens was 0.11 db 
and that of a double lens transformer approximately 
0.05 db. Therefore the launching loss was approximately 
1.5 db. Since the largest part of the measured transmis- 
sion loss was launching loss, it was not possible to deter- 
mine the iteration loss by adding or removing a phase 
transformer. The adjustment of the terminating horn- 
lens combinations was quite critical and the variation 
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of the launching loss was of the same magnitude or even 
greater than the loss of an individual phase trans- 
former. Less critical was the alignment of the phase 
transformers, although transverse displacements of the 
transformers divert the beam from its original path. 
No noticeable change in loss was produced by a slight 
tilt of the phase transformers, since a tilt causes only a 
second-order effect on the phase distribution. It was 
found that the spacing of the phase transformers could 
be varied within 10 per cent without affecting the total 
transmission loss. 


III. ITERATION Loss 


The iteration loss of a phase transformer consists of 
the diffraction loss and the reflection and absorption 
losses. The theoretical diffraction loss can be obtained 
from Fig. 2 as a function of (k/D) R®. The reflection loss 
pertains mainly to the phase transformers, which are 
uncompensated for reflection. The power reflection co- 
efficient I of one surface of a dielectric lens can be cal- 
culated approximately. For values of ¢>1.2 one obtains 

(4./e — € — 1) 


eae 

VYet1 4A evfe (e — 1)? 
giginitetener) 
(RD)? 1—e* 


where a= /k/D-R, and e¢ is the relative dielectric con- 
stant. The first term represents the reflection from a 
plane dielectric interface; the second term takes into 
account the curvature of the interface. The reflection 
loss of a single lens phase transformer, assuring no inter- 
action between the surfaces, is 


Lr = —2[10 logy (1 — I')|db. (3) 


The absorption loss within the dielectric depends not 
only on the dissipation factor of the dielectric material 
but also on the field distribution. One obtains for the ab- 
sorption coefficient the expression 


e,(a@? <1) 
Ve—1 


where e,=dissipation factor. The absorption loss is then 
given by 


, (4) 


n~ 


IDA = —10 logio (1 re! n) db. (5) 


The reflection and absorption losses are plotted indi- 
vidually in Fig. 8 as a function of the dielectric constant, 
assuming a 2.16 value for a. Two curves are shown for 
the absorption loss for values of ¢-=0.0001 and 0.0002. 
Curves representing the total dielectric loss of the phase 
transformer, the summation of reflection and absorp- 
tion loss, are also shown. 
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Fig. 8—Phase-transformer dielectric losses as a function of ~ 
dielectric constant for an a value of 2.16. 


In order to eliminate the launching variable and to 
obtain measured data as to the individual phase-trans- 
former loss, a parallel plate resonator was built using 
one section of the beam waveguide. This resonator con- 
sists of two plane reflector plates, one placed D/2 before 
and the other placed D/2 after the phase transformer 
(Fig. 9). This resonator behaves like a transmission-line 
section short-circuited at both ends. The resonant fre- 
quencies are approximately given by f=uc/2D, where 
n is an integer. By measuring the Q of the resonator, the 
loss of the phase transformer can be obtained? from the 
expression 


aD 
a One 8.68 db, (6) 


0 


where Q=f/Af, and Af is obtained from the half-power 
points of the resonance curve. The above expression 
does not consider a loss in the reflector plates which is 
indeed negligible. 

Since Q’s in the order of 20,000 or more were expected, 
a very stable oscillator was required. An atomicron 
stabilized Klystron oscillator which had a possible fre- 
quency variation of +3 Mc was employed to obtain the 
Q measurement. The coupling to the resonator was per- 
formed at the center of each plate through small slots. 
Measurements were made with the foam phase trans- 
formers. The double-lens transformer unfortunately 
could not be measured in the experimental set up be- 
cause its reflection compensation occurred at a fre- 


? See, for instance, E. L. Ginzton, “Microwave Measurements,” 
McGraw-Hill Book Co., Inc., New York, N. Y., pp. 469-470; 1957. 
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Fig. 9—Schematic of parallel-plate resonator. 


quency which was outside the range of the measuring 
equipment. In order to measure Q as a function of a 
the radius R of the phase transformer was varied by 
inserting aluminum plates with circular apertures. The 
points in Fig. 10 show the loss LZ calculated from the 
measured Q value as a function of a. Also shown as 
solid curves are theoretical iteration losses assuming 
dissipation factors of e¢,=0.0002, 0.0003, 0.0004. The 
measured values are in good agreement with the curve 
for ¢,=0.0003.% 


IV. Buttp-UP or REITERATIVE WAVE BEAMS 


The experiments conducted on the model beam wave- 
guide indicated that the wave beam had to pass through 
several phase transformers before the repetitive beam 
mode was established. In order to investigate the build- 
up of the reiterative beams, pulse-type measurements 
were conducted on a specially constructed beam wave- 
guide resonator (Fig. 11). This resonator differed from 
the one used in the preceding experiment, in that the 
reflector plates were, in effect, located within the center 
- planes of two adjacent transformers. The advantage of 
this arrangement is that the required phase transforma- 
tion.can be achieved by appropriately curving the re- 
flector plates rather than by using dielectric lenses, thus 
eliminating the dielectric losses. In other words, the line 
section simply consisted of two curved reflector plates 
facing each other as shown. 

The radii of curvature of the reflector plates coincided 
with their spacing which was approximately 30 meters. 
The diameter of each plate was 0.91 meters, correspond- 
ing to a value of $max=1.67. The reflector plates were 
machined from lucite sheet stock, the final surfaces be- 
ing plated with silver. The plates were mounted on steel 
frames, with a provision made on the back of these 
frames to allow a small deformation of the reflector 
plates. In this manner it was possible to adjust the aver- 


3 The manufacturer of the foam material, Emerson and Cuming, 
Inc., Canton, Mass., quotes an average dissipation factor of 0.0002 
for a frequency range of 10° to 3X10" cps. 
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Fig. 10—Measured and calculated iteration loss of_a foam 
dielectric phase transformer as a function of a=/(k/D)R. 
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Fig. 11—Experimental setup for determining diffraction loss. 


age curvature of the reflectors which varied slightly 
with temperature. 

A millimeter wave communication set was used as a 
signal generator, supplying 35 kMc pulses of 70 mysec. 
duration. The pulses were introduced into the cavity 
through a small coupling hole in the center of one of the 
reflector plates. The coupling hole, a square aperture of 
one cm on each side, was fed from standard rectangular 
waveguide. The same coupling scheme was used on the 
receiving side. The receiving system employed a crystal 
detector followed by a series of wide-band chain ampli- 
fiers, the output of which was displayed on an oscillo- 
scope. In order to view any portion of the series of re- 
ceived pulses, a variable delay network was inserted 
between the oscilloscope and the transmitter trigger. 

Fig. 12, a photograph of the oscilloscope display, 
shows the received pulse train. The display does not 


262 


have a flat base line because of the poor low frequency 
response of the receiving system. The pulse train en- 
velope has in the first portion a more or less irregular 
structure which eventually transforms into an exponen- 
tially decaying curve. The exponential portion indicates 
that the reiterative wave beam has been established. In 
Fig. 13(a) and 13(b) the exponential portion has been 
expanded to illustrate the identical shape of all succes- 
sive pulses. The original transmitted pulse had an odd 
characteristic shape which is precisely repeated from 
iteration to iteration, indicating the absence of any 
delay distortions. For Fig. 13(a) the reflector plates 
were spaced 30 meters apart so that the distance be- 
tween recurring pulses is equal to 60 meters. Therefore, 
in this photograph, the pulse train corresponds to a total 
path length of 4.5 km. The attenuation of the pulse 
train indicates a transmission loss of 1.75 db/km. This 
is equivalent to a loss of 0.052 db per iteration. The 
theoretical iteration loss, from Fig. 2 for (k/D) R*=1.6 7, 
is 0.02 db. Including a reflection loss of about 0.004 db 
caused by the finite conductivity of the reflector plate, 
(neglecting coupling loss at the coupling holes) the total 
calculated iteration loss is 6 0.024 db or about one half 
the measured value. This discrepancy, however, is not 
surprising considering the fact that the radii of curva- 
ture of the reflectors were not uniform, and that the 
surface condition because of difficulty in silver plating 
was poor. Also, the reflector surfaces could not be pol- 
ished because of poor adherence of the plated surfaces. 

A similar experiment was conducted with the reflector 
spaced 24.5 meters apart, the curvature of the surface 
wave substantially deformed to achieve an average 
radius of curvature of 24.5 meters. The pulse train for 
this measurement is shown in Fig. 13(b). The measured 


Fig. 12—Envelope of pulses reflected between the reflectors 
of Fig. 11 for a spacing D =30 meters. 
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Fig. 13—(a) Expanded pulse train of Fig. 12 with D=30 meters; 
(b) D=24.5 meters. 
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attenuation in this case corresponds to a transmission 
loss of 1.2 db/km, or an iteration loss of 0.029 db. The 
(k/D)R? value for this case is about 27. No theoretical 
diffraction loss figures have been as yet calculated for 
such high values of @max which could be used for quan- 
titative comparison. But this experiment does illustrate 
that the iteration losses can indeed be made very small. 
It must be said that the alignment of the reflector plates 
was not very critical and that no provisions were re- 
quired for fine adjustment. ~. 

The build-up portion of the pulse train depends very 
much on the curvature of the reflectors. According to 
theoretical expectations the simultaneous excitation of 
several beam modes should result in a gradual amplitude 
increase in the pulse envelope to a maximum, then a 
gradual transition into the purely exponential decay 
which occurs when the higher modes have been sub- 
stantially damped out. Fig. 12 is closest to the theoreti- 
cal expectations. If the curvature is not uniform or not 
of the correct radius of curvature the “effective” re- 
flecting surfaces for the various modes do not coincide, 
that is, the effective path lengths of the various modes 
are slightly different. This in turn results in different 
mutual phase shifts for each mode at each reflection, 
and causes the formation of maxima and minima of the 
pulse envelope as shown in Fig. 14(a). These maxima 
and minima are particularly pronounced in Fig. 14(b), 
where the spacing of the reflector plates was reduced to 
21 meters and the reflectors mechanically deformed. 


(b) 


Fig. 14—Build-up of reiterative wave beam from the setup of 
Fig. 11 for the case of nonuniform reflector curvature. 
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V. CONCLUSION 


The experiments support the theory on reiterative 
wave beams, and demonstrate the feasibility of guiding 
- electromagnetic wave beams with very little loss by 
- resetting the cross-sectional phase distribution at peri- 
_ odic intervals. They also indicate that the diffraction 
_ loss can be made small compared to the material losses 
of the phase transforming devices. The tolerance and 
alignment of the beam waveguide appear not to be very 
critical, although more detailed investigations in this 

respect are warranted. 
Since the beams can be deflected by means of reflec- 
tors or dielectric wedges, bends in a beam waveguide 
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can readily be performed without introducing any dis- 
tortions or additional losses. Because of the very rapid 
decay in the radial direction, the beam is almost com- 
pletely confined to a cylindrical space of the diameter of 
the phase-transforming devices and is not affected by 
obstacles outside this diameter. 
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Parasitic Excitation of Circular Antenna Arrays* 


T. L. SIMPSON}, MEMBER, IRE, AND J. D. TILLMAN{, MEMBER, IRE 


' Summary—A method of calculating the radiation pattern of a 
circular antenna array of parallel monopoles or dipoles using sym- 


- metrical components is presented. The results obtained when only 


one element of the array is excited are not particularly good. How- 
ever, if an additional element is added at the center of the array a 
fairly good directional pattern results. The sidelobe level is depend- 
ent upon the tuning of this central element, and a graphical analysis 
is given which optimizes this tuning. Calculated and measured pat- 
terns are shown and compared. Sidelobe levels of 15 db, and a beam- 
width of 78° are typical. The principal advantage of this antenna lies 
in the fact that the pointing of the beam can be changed by simply 
changing the element that is excited. 


INTRODUCTION 


been devoted to the development of antennas capa- 

ble of radiating directional beams which can be 
pointed in any desired direction without physical motion 
in the structure. In such an antenna the beam-steering 
is accomplished entirely by phase and amplitude control 
of the excitation currents. Although it is possible to 
achieve some control over the direction of pointing of 
a linear array, the angular sector which can be so 
covered is severely limited. If it is desired to point the 
beam to any azimuth angle whatever, some configura- 
tion involving circular symmetry is indicated. 


[: recent years, considerable engineering effort has 
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The circular geometry considered in this regard con- 
sists of a number of identical radiators uniformly dis- 
tributed along the circumference of a circle. If a satis- 
factory pattern can be obtained at all, then steering the 
beam could simply be a matter of rotating the excita- 
tion accordingly. The array characteristic would be in- 
dependent of the direction of pointing. With an array 
having a finite number of elements, an angular rotation 
of the beam would necessarily be accomplished in dis- 
crete steps depending on the number of elements in use. 

This paper contains an analysis of circular antenna 
arrays in which only one element of the array is driven. 
The configuration which is considered consists of m 
identical parallel dipoles or monopoles which are spaced 
uniformly around the circumference of a circle. An addi- 
tional element, not necessarily having the same half 
length, is located at the center. When an element in the 
ring is driven, a pattern with directivity in the plane of 
the ring results. Considerable control of the sidelobe 
level can be obtained by tuning the center element. A 
half-power beamwidth of 80° and a maximum sidelobe 
level of 14 db results from a typical design. 


DETERMINATION OF CURRENTS AND PATTERNS 


The configuration to be analyzed is shown in Fig. 1. 
Consider first the situation where all of the elements in 
the ring are driven by voltages which are symmetrical 
about the diameter through ¢=0. It can be shown that 
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such a symmetrical set of m excitation voltages can be 
resolved into (m+2)/2 or (m+1)/2 (depending on 
whether m is even or odd) sequence sets. For simplicity, 
m will be taken as an even number. The sequence volt- 
ages are 


Qa 
V,™ = V™) cos — ki, (1) 
m 
where 
k=sequence number, k=0, 1, 2-- + m/2, 
i=element number, 1=0, 1,2- +--+ m/2, 
V®) = kth sequence amplitude. 

a 
Ss 


Se 


<— Driven element 


Fig. 1. 


We wish to determine V) so that 


m/2 m/2 


24 
a7 Vs = 2m V) ere — ki. (2) 


Since there are (m/2)+1 values of 7, (2) represents 
(m/2)-+1 simultaneous equations which may be solved 
for the (m/2)+1 values of V™. The result is 


m—1 


1 
Wk) ee ee oS V get 2m lm)hi 


mM 1=0 


(3) 


A similar resolution into sequence components can be 
carried out for the currents, 

The linear equations describing the circuit properties 
of the array are written in terms of Z;;, the mutual 
impedance between elements 7 and j, and the self im- 
pedance. They are 


V; — sy LL oi. 


j=0 


(4) 


IRE TRANSACTIONS ON ANTENNAS AND PROPAGATION 


value of k, simple algebraic manipulation leads to 


V(*) = [@MZ®, (5) 
where 
ar 
FR Te >a Zoi COS — kt. _ (6) 
t=] m 


Z) is called the sequence impedance. 

It has been shown! that the far-field radiation pattern 
radiated by the kth sequence, with phase referred to the 
array center, is given by 


Fe = mI J insu) cos ko 


eo 


+ DS [jo PImp-(2) cos (mp — k)d 


p=1 


+ jm?t* I npr k(Z) Cos (mp + Delt 3 (7) 
where z= (27p/d) sin 6. The leading term has the form 
cos k¢ and is proportional to J;.(z), the kth order Bessel 
function of the first kind. This term is called the prin- 
cipal term, and the remaining terms are called residuals. 
The largest residual is proportional to Jm_x(z), and since, 
if the order is greater than the argument, the Bessel 
function is small, the residuals may often be neglected. 
Neglecting residuals, the pattern is then 


m/2 V &) 


) fPIx(2) cos ko, 


(8) 
which is valid for any excitation. If only the element 


4=0 is driven with voltage V,, the sequence voltages 
are, from (3), 


1 
V®)=—Y, t = 0, m/2 
m 
io 
=—V,: all otherz. (9) 
m 


Thus, driving one element only excites all the possible 
sequences. Substituting this into (8) gives 


= oi ae ec ae oer TAS) 
Haren) Zo Tg pa Z)(z) 


cos k¢. (10) 


This includes the first residual for k=m/2 which is 
equal to the principal term. Eq. (10) is general for any 
number of elements and any radius. 

On the basis of (10), extensive calculations were made 
for arrays of six elements with 27p/\ ranging from 1.5 
to 3.0, and with quarter-wave monopoles as elements. 
Results of spot calculations with a different number of 
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“The use of a ring array asa skip range antenna,” Proc, IRE, vol. 
43, pp. 1655-1660; November, 1955. 
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elements showed no significant differences. Throughout 
this range, the parasitic excitation resulted in a moder- 
ately directional pattern (half-power beamwidths aver- 
aged about 80° in azimuth), but was found to produce 
objectionably high sidelobes, particularly to the rear. 
Since no way was found to reduce the sidelobe level 
by varying the available parameters, a modification was 
sought that would achieve this result without sacrific- 
ing the circular symmetry. 

The patterns can be improved by including a central 
element which is not driven. Designating the self- 
impedance of this element by Z,., the mutual impedance 
to a ring element as Z,,, and the center current by J., the 
circuit equations for the array are 


Vel fort LiL, 


j=0 


(11) 


O= Zm >, Ty + TLee- 


g=0 


(12) 


Now, when (11) is substituted into (3), it is found that 
the addition of center element affects only the zero 
sequence. The new zero-sequence impedance is 


dee 


GOs FO 


(13) 
and the other sequence impedances are unchanged. 
The modified azimuthal pattern factor will be as given 
in (10) but with Z replaced by Z’ and with the con- 
tribution from the center element. The entire pattern is 


1 Vs 
BVP Naas (10. 
] Ji. (z) 
(k) 


+ 2 
fain OL 


cos Kal Ps) 


In (14), F, is the pattern factor of an isolated element of 
the array. 

The self impedance Z.- of the center element may be 
regarded as the sum of the actual self impedance and an 
arbitrary tuning reactance, X ,. If the center element is 
identical to the ring elements, then 


Lee = Zoo am ee (15) 


A method is now needed to choose X, to minimize 
sidelobe level. This will be illustrated for the case of 
m=6. For this case, Zm=Zo1. Defining the total omni- 
directional part of (25) as FP, the azimuthal pattern 
may be written 


3 
PF = fOr > Fe, 


(16) 
k=l 
in which, making use of (14) and (16), 
1 Zor 
(0) = VoF,| ——\| Jol -—"_)|, (17) 
4 ue iS ( @ Zo t+ jx 
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3 g 72 
Chih sis 
uu F VoF. E Pt J\(z) cos@ + 2 Hes J2(z) cos 2 
4° 
+2 oes J3(2) cos 30. (18) 


For convenience in plotting, the expressions of (17) and 
(18) are multiplied by 10% and normalized to Vo. 

For a particular choice of circumference the omni- 
directional term is a function only of Xz and the direc- 
tional term is a function only of ¢. Symbolically, (17) 
and (18) become 


no FO (Xz) and iS F“)() 
SG By ——= ; 
V : V : 


0 0 k=1 


respectively. Thus the omnidirectional component can 
be controlled independently from the directional part 
by tuning X;. 

To illustrate the optimization technique, consider a 
six-element ring of vertical \/4 monopoles 2.40) in cir- 
cumference. Assume a length-to-radius ratio h/a=75. 
Fig. 2 presents a complex plane plot of (17) and (18) 
for these parameters. The omnidirectional term is plotted 
as the locus of values obtained for — © <X,< ~ and the 
directional sum the locus obtained for 0°<@<180°; 
since the excitation is symmetrical the pattern will 
repeat for 180°<#<360°. As the omnidirectional term 
is independent of @, it adds directly to the directional 
term at each angle ¢. Thus, in Fig. 2 a particular 
choice of Xz, (and, hence, of F®) will result in a rigid 
translation of the directional locus. Construction lines 
are drawn in the plot to indicate the direction of trans- 
lation which would maximize the on-axis intensity and 
to define, also, that translation which would minimize 
the secondary radiation. The center of the “secondary 
radiation rectangle,” c, would be translated to the origin 
by the addition of its negative, c’. Evidently this would 
minimize the secondary radiation level. In this fortunate 
case, c is located very nearly opposite the origin from 
the ¢=0° point so that translating the directional locus 
to minimize the secondary level also nearly maximizes 
the on-axis level. By another fortunate circumstance the 
F® locus passes near the indicated optimal value for 
X, in the neighborhood of —75 ohms. 

In practice, the optimal value for Xz is determined by 
testing near-optimal values for secondary radiation 
levels. Fig. 3 shows the azimuthal patterns resulting 
from three ‘near-optimal values of X,0, —75, and —115 
ohms, Evidently —75 ohms produces the minimum side- 
lobe level. 

In order to determine pattern variation over a range 
of frequencies, complex-plane plots of pattern data sim- 
ilar to Fig. 2 were prepared for the arrays whose param- 
eters are tabulated in Table I. A frequency of 300 Mc 
was chosen for the case shown in Fig. 2. Patterns at 
330 Mc and 360 Mc were calculated using the best 
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value of X,, and beamwidths and sidelobe levels are 
included. 

The usable frequency band extends from slightly be- 
low to considerably above 300 Mc. Very poor patterns 
result at 240 and 270 Mc, with little or no improvement 
possible by tuning. Above 300 Mc, the pattern deterio- 
rates much less rapidly. 


o 
a 


NORMALIZED FIELD INTENSITY 


Tchebyscheff Tz pattern 
for-I5db side lobes 


60 0 100 120 
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Fig, 3. 
TABLE I 

Frequency | Circumference | X, Tale Rowen Side-Lobe 

Beamwidth Level 

300 Mc 2.4% —62 iN 78° ‘ nO 
330 Mc 2.64 d —163 g4° 6 as 
360 Mc 2.88 —180 96° = 8.3! db 
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MobDEL STUDY 


Since the results obtained in the preceding section 
for the 2.40 \ ring with center tuning were considered 
promising, this design was tested at a frequency of 300 
Mc. The elements were made of quarter-inch bronze rod 
which resulted in an h/a ratio near 75, the value used 
in calculations. Both the driven and center elements 
were constructed by extending the center conductors 
of 50-ohm coaxial line upward from below the ground 
plane. Center element tuning was accomplished using a 
calibrated coaxial tuning stub. * 

Fig. 4 compares measured with calculated azimuthal 
patterns of the 2.40 ring. Excellent agreement exists 
in the main beam, although some minor discrepancies 
appear at the lower levels—a calculated sidelobe of 
— 13.6 db at $=130° was measured as —12.9 db. Fig. 5 
compares vertical patterns, measured and calculated, for 
the array with good agreement. 

Experimental investigation at the frequencies of Ta- 
ble I verified the indications of calculated performance. 
Below about 290 Mc, the pattern deteriorated rapidly 


regardless of center element tuning. Above the design 
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frequency of 300 Mc, however, the pattern was much 
less sensitive to frequency. Figs. 6 and 7 show azimuthal 
patterns at 300 Mc and 360 Mc, respectively. With in- 
creasing frequency the beamwidth becomes wider and 
sidelobe levels increase. 

Impedance measurements were made at the base of 
the driven ring element. Fig. 8 compares the measured 
input impedance to the calculated value. While the 
pointwise correlation is only fair, the order of magnitude 
and the variation with tuning are in good agreement. 
Evidently it would be possible to tune out the reactance 
component of the input impedance by tuning the center 
element. However, this would disturb the pattern since, 
at this frequency, a negative tuning reactance is re- 
quired for optimal sidelobe reduction. Therefore, the in- 
put impedance must be tuned separately at the input. 
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CONCLUSIONS 


It was found possible to achieve a satisfactory azi- 
muthal pattern with a circular array having six ring 
elements if a tuned central element were introduced. 
A graphical analysis was developed whereby an optimal 
value of tuning reactance could be predicted. 

Parasitic excitation imposes definite limitations on 
the types of patterns obtainable. After the driven ele- 
ment is excited, the currents in the remaining elements 
are determined by the inherent properties of the array 
(except for the tuning control on the center element). 
Such a design is much less powerful, theoretically, than 
an approach involving a synthesis technique for specify- 
ing each elemental current to produce an arbitrary pat- 
tern. However, the simplicity of operation of the para- 
sitic array is so attractive that in many cases it is en- 
visioned that pattern requirements might well be relaxed 
slightly to allow equal competition between driven and 
parasitic circular arrays. 

If single-frequency operation is desired, a very simple 
switching system is possible to steer the beam. A six- 
position coaxial switch of the type which short circuits 
the undesired terminals can be connected to the base 
of the six elements in the ring through half-wavelength 
lines. The driven element is then selected by the position 
of the switch, and the remaining elements are effectively 
shorted at their base. 
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Angular Accuracy of a Phased Array Radar* 


L. E. BRENNAN}, MEMBER, IRE 


Summary—One type of phased-array radar of current interest 
employs an array of separate receiving elements, each followed by an 
individual amplifier. These individual signals are combined co- 
herently to form one or more receiving beams for searching, tracking, 
or performing both functions simultaneously. This paper presents an 
approach to the theory of angle measurement with a phased array of 
this type. 

In the one-dimensional problem considered here, the receiving 
antenna consists of a linear array of individual antenna-amplifier ele- 
ments. The receiver-noise-limited case is considered, in which ac- 
curacy is limited by the additive normally distributed noise present 
in each channel. An expression is derived for the limiting accuracy of 
angular measurement when a single set of samples is available. This 
set of samples is obtained simultaneously, one sample from each 
channel. Next, two methods of implementing the angular measure- 
ments are discussed. These are amplitude comparison monopulse 
and a coherent or phase comparison technique. For large signal-to- 
noise ratios and for either a square law or a linear envelope detector, 
the accuracy of amplitude comparison monopulse approaches the 
theoretical limit. The same accuracy can be achieved with the co- 
herent technique by proper weighting of the individual signals. 


SYMBOLS 
A;, Bi=amplitudes of the two quadrature compo- 
nents of Fi. 
A>, Bz=amplitudes of the two quadrature compo- 
nents of FE». 


N 
ai= >) cos (6+k5+#é). 
k=1 


d,= >, cos (6+k6—ké). 


k=1 


N 
b= Di sin (6+k5+8E). 


k=l 


N 
bo= >) sin (6+h6—2#Eé). 
k=0 
¢=an arbitrary constant. 
D;=difference signal; phase comparison, equal 
weights case. 
D,= difference signal; phase comparison, optimum 
weights case. 
d=spacing between elements of the array. 
i, E,=signals in the two squinted beams formed for 
amplitude comparison monopulse. 
ex = signal (voltage or current) in the kth channel, 
includes noise component. 
F(W) =a function defined in (37). 
k=number of the channel (an index). 


* Received by the PGAP, August 19, 1960. 
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Ei Vig oe yn| 5, ¢) =likelihood function for the 
sample of xz, ye, given 6 and ¢. 
N=number of channels. 
n, = noise component in kth channel. _ 
S=the sum signal (47), used to determine sense 
of the error signal in phase comparison mono- 
pulse. 


Pribar dea ee power ratio in one channel. 


t=time. 


N 
Ui= >> uz cos RE. 


k=1 


N 
U2= Ss Uk sin RE. 


kal 


Ux, Ve = amplitudes of the two quadrature compo- 
nents of 7. 6 


N 
Vi= >> uw cos RE. 


Kr=1 
N 

Vo= >> » sin RE. 
k=l 


w;, =set of weights used in forming difference sig- 
nal for phase comparison monopulse. 
Xk, Ve =amplitudes of the two quadrature compo- 
nents of eé. 
yy =reference phase at center of array. 


A,=difference signal; amplitude comparison 
monopulse, square law detector. 
A,=difference signal; amplitude comparison 


monopulse, linear detector. 
A;=difference signal; coherent signal processing, 
constant weights. . 
A,=difference signal; coherent signal processing, 
optimum weights. 
6=incremental phase shift between adjacent ele- 
ments of array due to target displacement 
from crossover axis. 
*=the estimated value of 6. 
6=angle between axis of array and line of sight 
’ to target. 
6’ =angle between crossover axis and line of sight 
to target. 
\= wavelength. 
£=incremental phase shift between adjacent ele- 
ments of array to provide monopulse squint 
angles. 
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Pi, P2=envelopes of signals in the two channels cor- 
responding to squinted beams for amplitude 
comparison monopulse. 

o=rms noise (voltage or current) in the indi- 
vidual channels. 
0,1=rms error in A; due to noise. 
os2=rms error in A, due to noise. 
0a3=rms error in A; due to noise. 
O,4=rms error in A, due to noise. 
o31=rms error in 6; amplitude comparison mono- 
pulse, square law detector. 
0s2=rms error in 6; amplitude comparison mono- 
pulse, linear detector. 
os3=rms error in 6; coherent signal processing, 
constant weights. 
os4=rms error in 6; coherent signal processing, 
optimum weights. 
=lower bound for the rms error in 6. 
o=reference phase of signal. 
y=NE. 


w=27-signal frequency. 


‘ I]. INTRODUCTION 


a HERE is increasing interest today in phased array 
Bl scans which use a set of separate antenna ele- 

ments, each followed by an individual amplifier, 
in place of a more conventional receiving antenna. These 
individual signals can be combined coherently to form 
several receiving beams simultaneously. The theory of 


angle measurements with this type of array is discussed 


here. For simplicity, only the one-dimensional problem 
is considered, and the antenna elements are assumed to 
be equally spaced in a linear array. 

In practice, the angular accuracy of this type of radar 
can be limited by any of several effects, including re- 
ceiver noise, external noise, component phase or ampli- 
tude errors, and atmospheric refraction. The following 
analysis considers the errors due to receiver noise, 1.€., 
noise originating in the individual antenna-amplifier 
channels. 


¢ 
Il. STATEMENT OF THE PROBLEM 


* 

A linear array of equally spaced antenna elements is 
assumed, each followed by an individual amplifier, as 
illustrated in Fig. 1. When a plane wave is incident at 
an angle @ from the axis of the linear array, the indi- 
vidual amplifier outputs (voltage or current) can be ex- 
pressed i in the form 


er = cos (wt + & + kb) + m, (1) 


where 
e, =signal in kth channel, 
w=2r-signal frequency, 
¢ =reference phase of signal, 
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6=incremental phase shift between channels, 
n; = noise in kth channel. 


The incremental phase shift, 5, is related to the angle of 
arrival by 


O-== —— C08 0, Z 
5 (2) 


where 


\= wave length, 
d=spacing between elements, 
6=angle of arrival measured from axis of array. 


For convenience, the following analysis considers errors 
in the incremental phase shift 6. These can be converted 
readily to errors in angle of arrival 0, using (2). 

It is assumed that the reference phase of the signal 
¢ is unknown, as is generally the case in radar systems. 
The effects of equipment errors, either random or sys- 
tematic, are not considered here, nor are errors in nor 
malizing the signals in the individual channels. 


Direction of arrival 


HA ike eh — <«—— Receiving antennas 


freee sist ‘ Bunonados a 
<+— Amplifiers 
e e 


Fig. 1—Linear phased array geometry. 


The noise components in the individual channels are 
assumed to be independent and normally distributed. 
It will be convenient to express them in the form 


Ne = Uz COS wl + d Sin wt, (3) 


where 


Ux, Ye= Quadrature components of noise in the kth 
channel, 


and 


Uy = Ve? = Ne 4 G. (4) 


The bar above a quantity is used to denote average 
value. The signal-to-noise power ratio in each individual 
channel is then 

S 1 


= ——- 5 
N 20° ) 
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In terms of the above definitions, the specific problem 
considered here is the following: Given a set of signals 
e,, obtained by sampling the outputs of the NV channels 


al 
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a log y' 
dg 


May 


so the expression for a joint estimate must be used. 
From Cramér,! the mean square error for all regular 
unbiased estimates of 5 has the lower bound 


ost? 


of a linear array simultaneously, how accurately can 
the incremental phase shift 6 be estimated? Since a 
single set of e, samples is assumed, this corresponds to 
the single-pulse accuracy of a pulsed radar. 

In Section III, a lower bound is obtained for the rms 
error in estimating 6. Sections IV and V consider two 
specific methods of estimating 6. These are amplitude 
comparison monopulse and a coherent phase comparison 
technique. 


II]. Lower Bounp ON ERROR IN 6 


The accuracy with which the incremental phase shift 
5 can be estimated from one sample of the set of signals 
e. is limited by noise. In this section, a theorem of sta- 
tistics is used to obtain a lower bound on the rms error 
in the estimate of 6. This theorem applies specifically 
to all regular unbiased estimates. Denoting the esti- 
mated value of 6 by 6*, an unbiased estimate is one for 
which the expectation of 6* is equal to 6, 1e., 


E(6*) = 6. (6) 
The conditions for regularity are more subtle and are 
discussed by Cramér.' 

It is convenient to express the signals, e,, in the form 


(7) 


Ck = Xp COS wt — yx SiN wi 


where x; and y, are the two quadrature components of 
one sample from the &th channel. The likelihood func- 
tion for the sample {x:, yx} is then 


N 


I 


Eon MNis 6, Teak rene 
XN yn | ?) (200%)® 


L(x, Seloy 2 


This expression follows directly from the assumption 
that the noise components in the individual channels are 
independent and normally distributed. 

Although we are concerned only with estimating 6, 
there are two unknown parameters in this case, 6 and ¢, 


» — cos (25 + $)]? + [yn — sin (k6 + 4)? 
[xn — cos (ké + #)]? + [yn — sin ( oN 


2c? 


1 Harold Cramér, “Mathematical Methods of Statistics, Prince- 
ton University Press, Princeton, N. J., ch. 32; 1946. 


B{(° log “yi Ai(? log L 
ag 06 


2 dlog L dlog L) 7? eB 
Wa erecra 
a¢ a5 
From (8) 
mee — +S ff, = cos (48 + 6)] sin (40 + @) 
— [y, — sin (kb + ¢)] cos (k6 + ¢)} (10) 

and 

a — + Saf [me — cos (25 + 6)] sin (48 + 9) 


(11) 


Since the noise components in the individual channels 
are independent of each other, as are the two quadrature 
components in each channel, we have 


E| [xm — cos (#8 + ¢)][y: — sin (6 + $)]} = 0 all &,/ (12) 


— [y, — sin (k6 + ¢)] cos (ké + ¢)}. 


and 
E{ [x — cos (R6 + $)] [x1 — cos (15 + ¢) ]} 
= E{ [y, — sin (k6 + ¢)][y. — sin (6+ ¢)]} =0 RA! 
=o? k=l]. 
(13) 


Combining (9)—(13) gives the following expression for 
the lower bound on mean square error in 6: 


120” 
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IV. AMPLITUDE COMPARISON MONOPULSE 


In amplitude comparison monopulse, two beams are 
formed, as illustrated in Fig. 2. Each beam is formed 
by summing all V outputs with the appropriate relative 
phases. The two resulting signals are envelope-detected 
and the difference of the two envelopes is calibrated in 
terms of angle from cross-over 6’, as shown in Fig. 2, 
The approximate incremental phase shift (or angle of 
incidence) must be known a priori to permit formation 
of two beams which are pointed roughly in the direction 
of incidence. A square law envelope detector is assumed 
in Sections IV-A through IV-C, and the linear detector 
case is discussed in Section IV-D. In both cases a large 
signal-to-noise ratio is assumed so that noise cross- 
modulation terms can be neglected. 
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Fig. 2—Amplitude comparison monopulse. 


In terms of the individual channel voltages, ex, the 
sum voltages for the two beams are 


1 
N 

= >> [cos (wt + o + #5 + RE) + m cos (wt + RE) 
1 


+ yy, sin (wt + ké)], 
; N 
Eo = Di ¢/—kE 
1 
N 
= > [cos (wt + ¢ + kb — hE) + u, cos (wt — RE) 
1 
3 + »% sin (wt — Ré)], (15) 
where 


5=incremental phase shift due to displacement of 
cross-over from incidence angle (6’ in Fig. 2) 
£=incremental phase shift for beam squinting. 


Since the noise voltages in the two channels are cor- 
related and will partially cancel in the difference signal 
A, it is important to retain information on the relative 


‘noise phases, 
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The sum signals can also be expressed in the following 
form: 


= A,coswi + By, sin wt 
fi, = Ascoswt + Bz sin wt (16) 
where 
Ay=a,+ U,+ V2 
By = —b,— U2+ Vi 
Az = 2+ Ui — V2 
B.= —bb+U2+NV1 (17) 
and 
N 
a, = >, cos (p + kd + é) 
1 
N 
a2 = >, cos (¢ + kd — RE) 
1 
N 
bi = Do sin (@ + hd + RE) 
1 
N 
by = >. sin (@ + kd — RE) 
1 
N 
Ui = >, um cos ké 
1 
N 
U2, = >> um sin ké 
1 
N 
Vi = > % cos RE (18) 
1 


N 
Vo = >, » sin RE. 
1 


In terms of these quantities, the outputs of the two 
square-law envelope detectors are 


A+ BY 
A? + B,’. 


Il 


(19) 


Finally, the error signal, which is denoted by A; in this 
case, is 
Ai = px” — pe’. (20) 
To obtain the error in the measurement of incremen- 
tal phase shift 6, two quantities must be calculated: the 
rms value of A; due to noise oa,, and the error slope or 
derivative of A, with respect to 6. The rms error in 
estimating 6 is then 


(21) 


TA; 
=) 
08 


Cos 
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A. The Error Slope—Square Law Detector 


First, we will obtain an expression for the error slope 
in the absence of noise. In this case, the difference signal 
is 

Ai = pi? — po? = (a1? + 512) — (a2? + bo®), (22) 
Substituting the sums of (18) for a1, a2, bi, and by, into 
this expression gives 

N N 
sin? — (6+ &) sin? oh (Ol aes) 


2 


we) (GP 


The following identities were used in summing the series 
of (18): 


Ai 


a 1 
cos — — cos (2 + 5)e 
p) 2 


N 
> sin ka —— p) 
1 |e 
2 sin — 
z 
: 1 ee: 
= sin (w+ Ja ~ sin 
> costa = (24) 
1 a 
2 sin — 
2 


Differentiating the expression for error signal A; of 
(23) gives the following value for error slope at cross- 
over: 


OA, 
06 | 5=0 


; Rey ee ee Mees 
N sin (VE) sin? (=) — sin (£) sin? (+) 


ee ers 


B. Derivation of o4,—Square Law Detector 


Next we obtain an equation for the mean square error 
in A; due to noise, 7.e., 


OA,” a (Ay a A,)?. 
From (17), (19), and (20) 


Jo Ay = 2[Uis(ar =a a2) aa Vo(ay = a») a. Uo(by of be) 
+ Vi(b, — b:)| + 4(UiV2 — UeV:). (27) 


(26) 


For the large signal-to-noise ratio case considered here, 
the second term of (27), which consists of noise cross- 
modulation products, can be neglected. The equation 
for 04,” is then 


oa,” = 4[(a1 — a2)?Uy? + (a1 + @2)?V2? + (b1 + b2)?U2? 
+ (be aa by)?V2 aE 2(a1 aaee 2) (by i bo) Ui1U 2 


=} 2(a1 = a2) (bs am bi) ViV 9]. (28) 
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Terms of the form U;V; are zero since the two quadra- 


ture components of noise in the individual channels are © 


independent. 


Next, consider the individual quantities occurring in — 


(28). The quantity U2 can be expressed as follows: 


ee Ee Ss) Fee ae 
UP = ( > um cos Ht) ( > ui cos i) (28) 
1 1 
> Ns ae ee 
= >) um? cos? ké 
1 
N 
= g? > cos? ké 
1 
a i 2N 1 
=e | ev tier act. (29) 
4 sin & 


Since the noise components in the different channels are 
independent, terms of the form a; in (29) were re- 
placed with 


RAl 


U,Ut = 


Un? = 0 all k. (30) 


The corresponding expressions for the other terms in 
(28) involving U; or V; are 


ee : in(2NV + 1 
[i ve=<| an + 1) | (31) 
4 sin & 
V2 =i? (32) 
and 
PATE Ae ape — cos (2N + -*) (33) 
ks ee 4 sin & 


Terms of the form (a@;+a2) or (6:1+62) in (28) are 
functions of 6. For simplicity, only the 6=0 case was 
considered here, which would be closely approximated 
in many cases with closed-loop tracking systems. Sub- 
stituting the expressions of (18) for a1, 01, a2, and by into 
(28) gives 


N 2 N 2 
oa,” = 16 | 7 ( » sin it) + Ta ( > cos 2) 
1 1 
nae Fey av N 
_ Us ( >> sin Ke) ( >> cos ke) | (34) 
1 1 


which after substitution of (29), (31), (32), (33) and fur- 
ther simplification becomes 


NE 
4N Une as + sin 2NVE — 2 sin NE 


P35) 


g 
sin & sin? — 
2 
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C. Error in Estumating 6—Square Law Detector 


Expressions have now been obtained for the error 
slope [see (25)| and the mean square error in A; due 
to noise [see (35) |. These can be combined using (21) 
to obtain the rms error in measuring incremental phase 
shift 6. The incremental phase shift &, which determines 
the monopulse squint angle has been retained as a 
parameter, and it is interesting to see how the error in 
6 varies with this quantity. Combining (21), (25), and 
(35), the mean square error in 6 can be expressed as 
follows when J is large: 


i oa, s 52 
oF ey oa FW) (36) 
06 / 50 
where 
y= NE 
and 
( ee ay 
4y sin? = + sin 2y — 2 sin v) 
FW) =¥3 - - (37) 


(v sin y — 4 sin? “) 


The function F() is plotted in Fig. 3. For Y=z, the 
squint angle is 1/2 beamwidth, 2.e., the centers of the 
two monopulse beams are separated by one beamwidth. 

_ Note that the function F(W) is constant for small y; 
i.e., the accuracy of angular measurement is independ- 
ent of squint angle in this region. This is explained by 
noting that as the squint angle decreases, the error slope 
decreases, while the correlation between the noise com- 
ponents in the two beams increases. The noise therefore 
cancels more completely in the difference signal. The 
error slope and noise in the difference signal A decrease 
proportionately so that the resulting error in 6 remains 
constant. 

For squint angles less than 1/2 beamwidth (1.e., 
Y=7), the value of F(y) is 24, and the corresponding 
rms error in measurement of 6 is 

Ce oes : (38) 

Nee 

This equation is valid for large N, and in this case the 
rms error in 6 approaches the theoretical limit discussed 
in Section III. When N is small, the mean square error 
in 6 approaches the limit given in (14) for small squint 
angles (i.e., a term of the form N*—WN occurs in the 
denominator). 


DD, Linear Envelope Detector 


The corresponding analysis for the case of a linear 
envelope detector will be outlined here. Again a large 
signal-to-noise ratio is assumed and the rms error in the 
difference signal is computed only for the case of 6=0. 
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30 


W~ monopulse squint angle 


F(W)~ mean square error in measurement 
of incremental phase shift, 5 


(e} 1/2 7 
W(Né) in radians 
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Fig. 3—The function F(y). 


With a linear envelope detector, the difference signal, 
which is denoted by A, in this case, is 


A> = pi — p2 (39) 
where pi and pp are defined as before [(15)—(19) ]. 
In the absence of noise, the difference signal is 
_ NG +s) elt 28) 
sin ————— sin ————— 
2 2 
As = (40) 
pol ae eee 
sin ; sin 


Differentiating this expression to obtain the error slope 
gives 


: Né Peet; E 
N sin — cos —— — sin —— cos — 
OAs 2 2 2 2 
e . (41) 
06 |s=0 : 
sin? — 


In addition to the error slope, one must also evaluate 
the rms error in A, due to noise. Using (39), (15), (19), 
and noting that for the large signal-to-noise ratio case 
cross-modulation products of noise can be neglected, 
gives 


Ne == 4/01 407 2a,(U 1 + V2) + 25:(U2 — V;) 
— Jag? + by? + 2ax(U1 — V2) — 2b2(U2 + Vi). 


(42) 

For the case of 5=0 and large signal-to-noise ratio this 

equation reduces to 

As — Ae 

Cia a2) =r Vi(b2— bi) + U2(bi+b2) ar V 2(a1+ a2) 
VJ/ar+b; 


- (43) 


QT 4 


As before, this expression is square and averaged to ob- 
Pano Ag: 

Again the rms error in 6, denoted by oi in this case, 
is the ratio of o,, and the error slope. The same end 
result for rms error in 4 is obtained in this case, namely, 
for the case of large signal-to-noise ratio and a linear 
detector 


o 


2N8 


Fy) (44) 


EV ears 


where F(y) is defined in (37). 


V. CoHERENT METHODS OF ANGLE MEASUREMENT 


Two coherent or linear signal-processing techniques 
are discussed in this section. In both cases, a difference 
signal proportional to the incremental phase shift 6, is 
formed by combining the individual signals e, coher- 
ently (i.e., at RF or IF where phase is preserved). In the 
first case, the signals from each half of the array are 
summed individually with equal weights. The difference 
between these two partial sums is proportional to the 
incremental phase shift 6. This is essentially a phase 
comparison monopulse technique, the difference signal 
being proportional to the difference in phase of the two 
partial sums. The second technique, discussed in Sec- 
tion V-B, is similar, but weighs the different signals e:, 
unequally. 


A. Constant Weights 


First consider the case where the signals are weighted 
equally and the difference signal has the form 


N N/2 
Da foe ee ities 
1 


N/24+1 


(45) 


For convenience, we assume the array contains an even 
number of elements, V. Again, the approximate angle of 
incidence must be known to steer the beams formed by 
the two halves of the array. From (1) and (45), the 
signal component of the difference signal is given by 


N6 
sin* ——=| : 
D; = ; 4 | sin (wi + y) = Assin (wt + y), (46) 
sin 6/2 
where 
N+1 
y=¢ 6 = reference phase at center of array. 


The amplitude of the difference signal A; is a function 
of the magnitude of the incremental phase shift 6, while 
the sense of the error signal is determined by the RF 
phase of D;. In this case, a sum signal must also be 
formed to determine the sense of the error signal 
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N 
S= >> & = K cos (wt + y). (47) 
1 


For small pointing errors, 7.e., small values of 6, the 
quantity K is approximately equal to the number of 
channels, NV. Of course, the sum signal S also contains a 
noise term, so that at very low signal-to-noise ratios, 
errors in sign of the error signal may occur frequently. 
Again we will consider only the large signal-to-noise 
ratio case, and assume that the amplitude of the 
sin (wt-+y) component of D3 is obtained coherently. 

In this case, the rms noise in the difference signal A; 


is 
oss = V/N co. (48) 
The error slope in this case is 
N6 OS SeGe Oe ING 
sin? — cos 6/2 — — sin — sin —— 
0A; 4 2 2 2 
(—)- : (49) 
06 sin? 6/2 
For small values of 6, this reduces to 
(=) = (50) 
8 faze 20 


The rms error in 6 is again obtained by dividing the rms 
error in the difference signal by the error slope, 1.e., 


TAs 4o 


0A3 S N3!2 
Coe 


The limiting accuracy of (14) is not achieved in this 
case when N is greater than 2. For large N, the rms error 
obtained with the coherent technique described here is 
13 per cent greater than the error for amplitude com- 
parison monopulse. 


(51) 


an 


B. Variable Weights? 


In the preceding subsection, it was shown that when 
the signals are all weighted equally, the limiting accu- 
racy (see [14]) is not achieved. From the following 
heuristic argument, one would expect an improvement 
in accuracy when a set of unequal weights is used. The 
error signal obtained by subtracting the two end signals 
from a multi-element array (ey —e1) is greater for a given 
6 than the error signal from the two middle elements of 
the array. However, the rms noise in these two difference 
signals is the same. Hence, a set of variable weights 
which favor the end elements of an array should improve 
the ratio of error-signal to noise. In the following para- 
graphs it will be shown that this is indeed the case, and 
that with an optimum set of weights the limiting accu- 
racy is achieved. 


2 The possibility of improving accuracy with a set of variable 
weights was suggested by I). L. Margerum, of Systems Labs. Corp. 


| 
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In the general case of variable weights, the difference 
signal Dy, is 


N/2 
Des). 


1 


fue N+1 
( » 2w, sin (: — ———\s) sin (wt + y) (52) 
1] 


a 


We len—K+1 = ex.| 


Ag sin (wt + y). 


Again, N is assumed to be an even number, and the 
large signal-to-noise ratio case is considered. As before, 
the sum signal is used to determine the sense of the 
error signal and to extract coherently only one quadra- 
ture component of signal plus noise. 

The mean square noise in the error signal A, is then 


N /2 


oa? = 20? Dd, wr’. (53) 
1 


Again considering the case of small 6, the error slope 
obtained by differentiating A, of (52) is given by 


) 
; : 


Next, one would like to find the set of weights wx which 
minimizes the ratio of os to the error slope, 1.e., the set 
for which 


N/2 


peels 


OAg 
06 


(54) 


OA, 


0 2 SSS 

a (=) = 0 Fa ee 

OW; —— 
06 


This optimum set of weights is proportional to distance 
from the center of the array, 1.¢., 


1h ie la 
Wr = (e- =) 


where c is a constant. 
From (53), (54), and (56), the corresponding mean 


square error in 6 is 


(56) 
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N3 = N 


(57) 


(en) = 
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This is the same mean square error as was obtained in 
Section III for the limiting accuracy and again in Sec- 
tion IV for amplitude comparison monopulse. 


VI. CONCLUSIONS 


A lower bound was obtained for the rms error in angle 
measurement with a linear phased array for the single 
sample, receiver-noise-limited case. It was shown that 
this limiting accuracy is achieved with amplitude com- 
parison monopulse for the more tractable case of large 
signal-to-noise ratios. It was also shown that with 
amplitude comparison monopulse, angular accuracy 1s 
independent of squint angle over a wide range of squint 
angles—from zero to approximately 1 /2 beamwidth. 

When the individual signals are processed coherently, 
as discussed in Section V, and weighted equally, the 
angular error exceeds the theoretical limit by 13 per 
cent. However, when an optimum set of weights is used, 
the limiting accuracy is again achieved for large signal- 
to-noise ratios. It was shown that optimum weighting 
in this case implies weights proportional to distance 
from the center of the array. 

Some of the results obtained here are also valid at 
low signal-to-noise ratios. The expression for limiting 
accuracy derived in Section III is valid at all signal-to- 
noise ratios. However, the results for amplitude com- 
parison monopulse are invalid at low signal-to-noise 
ratio, since the noise cross-modulation terms in the en- 
velope detector output were neglected. Noise cross- 
modulation would predominate at low signal-to-noise 
ratios, resulting in rapid degradation of angular accu- 
racy. The analysis of Section V for the coherent signal 
processing case is valid at low signal-to-noise ratios, 
provided an accurate reference or sum signal is avail- 
able. However, at low signal-to-noise ratios the reference 
signal would also contain a large noise component, re- 
sulting in frequent errors in the sign of the error signal. 

When closed-loop angle-tracking is used, the accuracy 
would exceed that estimated here for the single sample 
case. In effect, with closed-loop tracking many samples 
are integrated to obtain a current estimate of angle. Al- 
though the closed-loop case was not considered explic- 
itly, the results obtained here provide an essential input 
for the analysis of linear arrays employing this tech- 
nique. 
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A Spiral-Doublet Scanning Array* 
JOHN R. DONNELLAN{ 


Summary—An array of eight spiral doublets, each doublet con- 
sisting of a pair of equally-excited two-wire spirals in the same plane 
and wound in opposite senses, has been used to scan a beam of hori- 
zontal polarization over a +40° range simply by rotations of the spiral 
doublets. The combination of a microstrip ring network and twin lead 
provided a balanced feed for each spiral in the array. Sidelobe levels 
were those predicted (25 db) for the broadside pattern and—at the 
test frequency of 1430 Mc—were, in general, below 19 db over the 
80° scan. A novel microstrip feed harness provided isolation between 
the doublets and served to absorb most circulating currents. This 
printed circuit feedboard consisted of judiciously-grouped microstrip 
ring networks and provided the desired feeding coefficients to within 
+0.1 db. An application of the array for tracking purposes is also 
described. 


INTRODUCTION 


T has been shown!” that two equally excited spirals 
] of opposite sense, placed side by side in the same 

plane—a spiral doublet—radiate a combined field 
which is everywhere linearly polarized. The direction of 
polarization is controlled by rotating one spiral of the 
doublet relative to the other. The phase of the radiated 
field is changed, while holding direction of polarization 
constant, by rotating both elements of the doublet equal 
amounts but in opposite directions. 

A beam scan of linear polarization from an eight 
doublet array of spirals has previously been obtained.” 
The direction and width of each main beam was in very 
close agreement with theoretical predictions, but side- 
lobes were excessive. Moreover, the beam scanned suc- 
cessfully 25° more to one side than to the other. Each 
spiral in this array had been excited from a coaxial 
line, z.e., an unbalanced feed. The spiral on the other 
hand is a symmetrical device which indicates the need 
for a balanced input. Further reduction of sidelobe 
levels and over-all improvement of performance of the 
array appeared to depend to a large extent upon im- 
provement of individual spiral performance. Also, an 
improved method for obtaining correct phasing between 
spiral doublets was deemed necessary. With these needs 
in mind, it was decided that a new and improved array 
would be constructed. 


* Received by the PGAP, July 1, 1960; revised manuscript re- 

ceived, October 3, 1960. 
_ + Microwave Antennas and Components Branch, Electronics 

Div., U.S. Naval Res. Lab., Washington, D. C. 

1]. A. Kaiser, “Scanning Arrays Using the Flat Spiral Antenna,” 
Naval Res. Lab., Washington, D. C., Rept. No. 5103; March, 1958. 

2 J. A. Kaiser, “The Archimedean two-wire spiral antenna,” IRE 
TRANS. ON ANTENNAS AND PropaGaTION, vol. AP-8, pp. 312-323; 
May, 1960. 

3J. A. Kaiser, “Spiral Antennas Applied to Scanning Arrays,” 
presented at Electronic Scanning Symp., AF Cambridge Res. Ctr., 
Cambridge, Mass.; April, 1958. 


SPIRALS AND ASSOCIATED COMPONENTS 


Fig. 1 depicts the geometry of the Archimedean spiral 
used in the experimental work described in this report. 
A spiral doublet, therefore, consists of two-spirals in the 
same plane, one wound in the sense shown in Fig. 1, 
the other wound in the opposite sense. An array of such 


doublets leads to an antenna of an arbitrary linear polar- 


ization which can be made to scan by changing the rela- 
tive phase between the spiral doublets. 

A front view of the eight-spiral doublet array used in 
the experimental work is shown in Fig. 2. The spacing 
of the spirals above the ground plane is 0.25 \ and the 
spacing between doublets (and also between two spirals 


Fig. 1—Archimedean spiral antenna, 


Fig. 2—Eight-spiral doublet array. 
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of each doublet) is 0.5 (10.5 cm) at 1430 Mc, which was 
the test frequency. All spirals along the top row are of 
one sense while those along the bottom row are of the 
opposite sense. This array was intended to scan only in 
the horizontal plane with an arbitrary linear polariza- 
tion; horizontal polarization was chosen for the initial 
experimental work. 

Each spiral was fed with 150-ohm twin lead, one end 
of the twin lead being connected to the input terminals 
at the center of the spiral while the other end was con- 
nected to a microstrip ring network (Fig. 3). Attaching 
a two-wire line at points B and C assures that the input 
terminals of the spiral are fed 180° out of phase, a 
condition necessary for optimum spiral performance. 
Thus the rat-race provided a transition (balun) from 
an unbalanced line to a balanced line, and the combina- 
tion of ring network and twin lead provided a balanced 
feed for the spiral antennas. This balanced feed elimi- 
nated the “squint” which had been associated with the 
coaxially fed spirals. 

Aquadag was applied to the areas between the fila- 
ments of the last quarter-turns of each spiral in order to 
absorb any currents which reached the outside edges of 
the spirals. Since the Aquadag was located beyond the 
area from which radiation occurred, there was no effect 
on the currents reaching the radiating region from the 
input terminals. The Aquadag served to attenuate cur- 
rents which, if allowed to be reflected from the ends of 
the spiral, would radiate in a sense opposite to that of 
the currents traveling from the input terminals. These 
two opposite-sense currents would result in an ellipti- 
cally polarized radiated field from the spiral, the degree 
of ellipticity being a function of the relative magnitudes 
of the two currents. Thus Aquadag was employed to im- 
prove the axial ratio of the field radiated from each 
spiral. An axial ratio of 1.5 db with no Aquadag painted 
on the spirals was reduced to 0.5 db when Aquadag was 
applied. 


FEED STRUCTURE 


The novel feed harness used to excite the array is 
diagrammed in Fig. 4. It was composed of judiciously 
grouped microstrip rat-races of the same general con- 
struction as the balun of Fig. 3. The board was carefully 
designed to give the proper Tchebycheff illumination— 
into matched loads—for 25-db sidelobe levels. An 
analysis of design considerations follows. 

Power is fed into the center rat-race at point 1 (Fig. 
4) and divides equally at points 2 and 4, the currents at 
these points being 180° out of phase. The initial hybrid 
network is thus used to divide the input power into two 
equal amounts, one half going to one side of the feed- 
board, the other half going to the opposite side (a 
mirror image). If there are any reflected currents at 
points 5 and 6, these currents will arrive back at points 
2 and 4 in phase and will be absorbed by the matched- 
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load terminating arm 3. For convenience of description, 
the board has been divided into eight sections, I to 
VIII. Consider section IV. Power entering this section 
at point 6 divides equally at points 7 and 8, the currents 
at these points being 180° out of phase; these currents 
then continue on to points 10 and 11. If line length /’ 
is equal to line length /’’, no power will enter arm 12, but 
all will arrive at point 13. If /’/=/ +0.5\X, all the power 
will enter arm 12, none arriving at point 13. If /=1” 
+0.25X, the power divides equally between 12 and 13. 
Therefore, any desired power distribution between 
points 12 and 13 can be obtained by proper choice of 
|i’ —1''|. The power distribution for values of |v —1"| 
ranging from 0 to 0.5 insteps of 0.05 \ was measured 
at 1430 Mc: 


W=Oot, 2.0 210: 


The normalized curve of the experimental results is 
shown in Fig. 5 along with the theoretical curve. The 
maximum error is 5 per cent. The experimental values 
were used in the design of the feedboard. With this in- 
formation, any desired fraction of the power could be 
tapped off at point 12 (Fig. 4), the remaining power 
continuing on from point 13 into section III. Here, 
again by proper choice of line lengths, the desired frac- 


Fig. 3—Microstrip ring network balun. 


Fig. 5—Relative power transmission (theoretical and experimental) 
vs path difference |J’—1'’| for scheme shown. 
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tion of the power could be tapped off at point 17, the 
remainder continuing on from point 18 into section II, 
etc. In this manner, the proper feeding coefficients for 
the array could be obtained. 

Again, consider section IV of the feedboard. Power 
entering the center rat-race at point 12 divides equally 
between arms 14 and 15. Coaxial cables were attached 
at these points in order to feed the spiral doublet. Any 
reflections at points 14 and 15 (or beyond these points, 
down the transmission line to the spiral antennas them- 
selves) will be absorbed by the matched-load terminat- 
ing arm 16. A similar argument holds for the other seven 
doublets. The output arms of the feedboard used to feed 
doublets 1 to 8 are indicated in Fig. 4 by S; to Ss respec- 
tively. A similar extended analysis will show that this 
feedboard provides isolation between the doublets and 
serves to absorb most of the circulating currents by 
series of dumpings into matched loads. 


EXPERIMENTAL RESULTS 


The spiral array was phased by using it as the receiver 
while a pyramidal horn oriented for horizontal polariza- 
tion was used as the transmitter. The transmitter and 
receiver were separated by a distance of 72 feet (ap- 
proximately 3X2d?/d). A bolometer detector was con- 
nected to the output arm of the feedboard (arm 1 in 
Fig. 4) and then by coaxial cable to an amplifier. With 
the receiver directly facing the transmitter, all doublets 
in the array were first set for horizontal polarization. 
The next step was to rotate all doublets by the proper 
amount to obtain an on-axis in-phase condition. When 
the array had been phased, the pattern in the horizontal 
plane was recorded (Fig. 6). The pattern covers the 
range from +108° to —108°. Sidelobe levels obtained 
(over a +90° range) are those predicted. The positions 
of the lobes are also in close agreement with the pre- 
dicted positions, while the 3-db main beamwidth (14°) 
is two degrees narrower than the predicted width (16°). 

It had previously been determined that in order to 
direct the beam to some angle 0, the phase difference 
between adjacent doublets had to be equal to approxi- 
mately 30 (more accurately, 180 sin 6); this 36 relation- 
ship holds for scan angles up to approximately 45°. Thus 
Fig. 7, in addition to the broadside pattern, shows the 
+10° beam displacements obtained by introducing a 
phase difference of +30° between adjacent doublets. In 
one case, the upper spirals of doublets 5—8 were rotated 
clockwise 15°, 45°, 75°, and 105°, respectively, while 
the upper spirals of doublets 4, 3, 2, and 1 were rotated 
counter-clockwise 15°, 45°, 75°, and 105°, respectively; 
in the other case, the rotation directions were reversed. 
Rotation of the upper spiral of a doublet in one direc- 
tion through an angle a causes a corresponding rotation 
of the lower spiral in the opposite direction through an 
equal angle a. Thus constant polarization was main- 
tained while only the phase between doublets was 
changed. In this manner, the phase center was balanced 
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Fig. 6—Experimental broadside pattern 
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Fig. 7—Patterns produced by phase differences between doublets of 
0° and +30° (main beam directed to 0° and +10°, respectively). 
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Fig. 8—Patterns produced by phase differences between doublets of 
0°, +30°, +60°, +90°, and +120° (main beam directed to 0°, 
+10°, +20°, +30°, and +40°, respectively). 
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about the center of the array, and scanning was achieved 
by simple proper rotations of the spiral doublets. In a 
similar manner, beam displacements of +20°, +30°, 
and +40° were obtained as indicated in Fig. 8. 
Table-I shows the variation of maximum sidelobe 
level (over a +90° range relative to the broadside posi- 
tion) and 3-db main beamwidth with angle of scan 
(0° to +40°). In general, the sidelobes are increasing 
with increasing scan; the average sidelobe levels, how- 
ever, are considerably below the figures given in the 
table (e.g., all sidelobes other than the maximum at 
scan angles + 40° are greater than 19 db down from the 
main beam). The 3-db main beamwidths are, in general, 
2° narrower than the theoretically predicted beam- 
widths, with both increasing with increasing scan angles. 


TABLE I 


VARIATION OF MAXIMUM SIDELOBE LEVEL AND 
3-DB MAIN BEAMWIDTH WITH SCAN ANGLE 


Scan An le : Maximum . 3-db : 
(De me ) Sidelobe Level Main Beamwidth 
ates (db) (Degrees) 

0 25.0 14.0 
+10 PIED) 14.5 
+20 20.0 14.5 
+30 20.0 16.5 
+40 Aes 18.0 


At scan angles greater than +40°, sidelobes were ex- 
cessive (e.g., the maximum sidelobe level at +45° was 
16 db down from the main beam). The cross-polarized 
broadside pattern showed lobes 20 db down from the 
main beam. 

An additional feature of this array can be seen in 
Fig. 9. Shown here again is the broadside pattern as ob- 
tained with the detector at arm 1 (see Fig. 4) of the feed- 
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Fig. 9—A broadside pattern and a split-beam pattern resulting 
from a 180° phase shift delivered to one side of the array. 


board, as described previously. If the detector and the 
load-terminating arm 3 are switched, the sharp null 
shown in Fig. 9 is obtained. This null occurs because a 
180° phase shift is delivered to one side of the feedboard 
when the detector and load are switched, so that the 
four doublets on one side of the feedboard are 180° out 
of phase with the four doublets on the other side. The 
pattern thus obtained would be useful for tracking pur- 
poses. 
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Leaky Wave Antennas II: Circular Waveguides* 


L. 0. GOLDSTONE}, SENIOR MEMBER, IRE, AND A. A. OLINER}, FELLOW, IRE 


Summary—The propagation characteristics of leaky waves on 
slotted circular cylinders are derived by the methods described in an 
earlier paper on leaky rectangular waveguides. The cross section of a 
typical leaky circular waveguide is represented by a radial transmis- 
sion line network incorporating a lumped terminating admittance 
which characterizes the slot discontinuity and the external region. 
The expression for the lumped admittance is obtained by combining 
an integral equation solution with a variational procedure. The reso- 
nances of this transverse network, which yield the complex propaga- 
tion constants for the leaky waves, are solved by perturbation tech- 
niques to produce results in simple and practical form. Solutions are 
obtained for leaky waves corresponding to two different excitations 
of a slotted cylinder and these results are shown to compare favor- 
ably with measured values. 


I. INTRODUCTION 


LUSH-MOUNTED antennas of the type classed 
is leaky waveguides radiate all along the length of 
- the waveguide and can be characterized by a com- 
plex propagation constant along the waveguide. The 
characteristics of leaky wave antennas and their analy- 
sis by microwave network techniques have been dis- 
cussed in an earlier paper, where applications were made 
to leaky waveguides of rectangular cross section.* It is 
the purpose of this paper to present the application of 
such network methods to the analysis of slotted circular 
waveguides of the type shown in Fig. 1. 


Fig. 1—Slotted circular waveguide. 


The representation of a slotted circular waveguide by 
a transverse equivalent network involves a nonuniform 
transmission line of the so-called radial type. The modes 
required for such a radial transmission line representa- 
tion are E-type and H-type modes which propagate in 
the radial direction but are distinguished by vanishing 
longitudinal (axial) components of the magnetic field 
and electric field, respectively. These modes are similar 
to the modes employed in connection with rectangular 
leaky waveguides! in that no coupling between the two 
types can occur in longitudinally uniform structures. A 
discussion of these modes is presented in Section II. 


* Received by the PGAP, July 7, 1960. The research described in 
this paper was performed at the Microwave Res. Inst. of the Poly- 
technic Inst. of Brooklyn under Contract No. AF 19(604)-2031 with 
the AF Cambridge Res. Ctr. 

% ce nave Res. Inst., Polytechnic Inst. of Brooklyn, Brooklyn, 

1L. O. Goldstone and A. A. Oliner, “Leaky wave antennas I: 
rectangular waveguides,” IRE TRANS. ON ANTENNAS AND PROPAGA- 
TION, vol. AP-7, pp. 307-319; October, 1959. 


The parameters of the lumped circuit portion of the 
transverse equivalent network which characterizes the 
slot discontinuity and the external region are obtained 
by a combination of small-aperture results, which are 
derived in the Appendix, and a variational procedure. 
The Appendix treats the diffraction by the slot of a 
radial wave arising from within the slotted cylinder, for 
two different exciting fields. The diffraction problems 
are phrased in terms of integral equations for the slot 
electric field; solutions for the slot fields and for the slot 
magnetic polarizabilities are obtained in the static small- 
aperture limit. These polarizabilities are then used in 
Section III in the evaluation of the slot discontinuity 
susceptances, while the slot fields are employed in varia- 
tional expressions for the radiation conductances. 

The calculation of the leaky wave propagation con- 
stants is accomplished by employing a perturbation 
technique for the solution of the transverse resonance 
equation which is similar to that used for rectangular 
leaky waveguides. 

Two types of leaky waves in slotted circular wave- 
guides are considered in this paper. These are the E- 
type and H-type leaky waves which can be regarded as 
perturbations of the Eo and Hu modes, respectively, in 
a closed round waveguide. The theoretical results agree 
very well with measured values for both types of leaky 
waves when the slot is narrow. The agreement is not as 
good for very wide slots because a perturbation proce- 
dure was used and because the lumped susceptances oc- 
curring in the transverse networks are evaluated in the 
small-aperture limit. The problem of calculating the 
propagation constants of these two leaky waves has 
been treated previously by Harrington.? Although there 
exist similarities between his method and the network 
procedure described here, there are also significant 
differences between them. The advantages of the net- 
work procedure have been discussed in the previous 
paper.! It might, however, be mentioned that the net- 
work treatment results in expressions for the propaga- 
tion constant which are much simpler to compute from 
than those of Harrington and which, furthermore, are 
explicit in form. 

As mentioned above, the diffraction problems treated 
in the Appendix also yield the magnetic polarizabilities 
of an infinite slit in an infinite cylinder for both the 
longitudinal and transverse components of the magnetic 
field. It is interesting to compare these results with those 
for an infinite slit in an infinite conducting plane. The 


2R. F. Harrington, “Propagation along a slotted cylinder,” 
J. Appl. Phys., vol. 24, pp. 1366-1371; November, 1953. 
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results in the literature for the plane are available only 
for the case for which there is no field variation along the 
slit. Under these conditions, the transverse magnetic 
polarizabilities, which correspond to weak excitation of 
the slits, are identical when the width of the slit in the 
plane is replaced by the arclength of the slit in the cylin- 
der. For the longitudinal magnetic polarizability, which 
corresponds to strong excitation of the slits, the results 
for the plane and the cylinder are slightly different from 
each other. Further details are included in Section III. 


II. RaprIAL TRANSMISSION LINE MopEs 


Although the cross sections transverse to the radial 
direction in circular cylindrical coordinates are non- 
uniform, it is nevertheless possible to describe the fields 
transverse to the radial direction in terms of mode func- 
tions which possess orthogonality properties over the 
transverse cylindrical surfaces. In terms of such modes, 
circular cylindrical waveguide structures may be repre- 
sented by transverse, radial transmission lines. 

The field expansion employed here differs from that 
presented in the “Waveguide Handbook.”* The mode 
function concept, as discussed there, requires that the 
mode functions be independent of the coordinate along 
the transmission direction, with the result that a vector 
modal representation is not possible for radial geometry. 
However, if the mode function concept is modified to 
permit a variation of the mode functions with 1, it be- 
comes possible to derive a vector modal representation 
with orthogonality properties over the cylindrical $z 

cross section. Such a representation is discussed below 
for the case where there is an arbitrary exponential 
variation in the axial (zg) direction. 

Using such radial mode functions, a representation of 
the transverse (to r) vector fields may be written as* 


E,(r, 6) = 2) Vn'(r)en'(r, 6) + Vn’ (ren (7, ¢), (la) 
Hi(r, 6) = 22 In'(r) Bn’ (r, 6) + In! (1) Bn'(r, 6). (1b) 


The primes in the above representation denote the E- 
type modes for which H,=0, and the double primes de- 
note the H-type modes for which E,=0. 

A set of mode functions appropriate for radial trans- 
mission line representations has been discussed else- 
where by Altschuler and Goldstone.’ The radial trans- 
mission line equations, characteristic impedances and 
radial propagation wavenumbers given in this reference 
are the same as those to be employed here. However, it 
should be noted that the mode functions which are to 


_ 3N. Marcuvitz, “Waveguide Handbook,” M.1I.T. Rad. Lab. Ser., 
McGraw-Hill Book Co., Inc., New York, N. Y., vol. 10, Section 1.7; 
19517 ° 

4 The exponential z dependence has been suppressed. 

5 H. M. Altschuler and L. O. Goldstone, “On network representa- 
tions of certain obstacles in waveguide regions,” IRE TRANS. ON 
MicrowaVE THEORY AND TECHNIQUES, vol. MTT-7, pp, 216-217; 
April, 1959. 
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be used in this paper differ from those given in the refer- 
ence cited in that the ¢ components eng and ng of the 
mode functions used here contain a factor of 1/r which 
does not appear in the corresponding functions given 
in the above reference. The mode functions as defined 
here possess orthogonality properties over the cylindri- 
cal surfaces transverse to 7 and are normalized as fol- 
lows: 


Qa Qa 
f Apo Cre en rae = f hy!” Xr07 en TEP = Onn, 2) 
0 0 


Qr 24 
i hy!* X r0:em! 1rdd = f hi’® Xoro-en rag = 0," GC) 
0 ® 0 


where ro is a unit vector in the radial direction and the 
asterisk denotes the complex conjugate. 

Since both the cases which are treated in this paper 
involve longitudinal field components which are even 
functions of the angular coordinate ¢, the required mode 
functions become: 


E-type modes: 


en 
Ca (/~ cos nd, (4a) 
T 
ia 0: (4b) 
fen cos upd 
i= = — ) (4c) 
T r 
ikz nm sin np 
Cant J = ) (4d) 
hk? — ko? ¢ Va 
H-type modes: 
ae 
Ente a /= cos ng, (Sa) 
Tv 
or = 0, (5b) 
pie En cos nd (5c) 
vg r 
Re nm sin n 
Vase a 2 so pI (5d) 
Ro — kh. 4 Var 
where 
nm = 0,1, 2, 
and 
e, = Ly 2. n=0 
ibe n> 0. 


The orthogonality of the mode functions given in (4) 
and (5) is easily verified by substitution into (3). 
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Il]. TRANSVERSE EQuIVALENT NETWORKS 


The two “lowest” leaky waves in slotted round wave- 
guide can be regarded as perturbations of the Hy, and 
Ey. modes in a closed guide. Each of these cases poses a 
different transverse discontinuity problem for the evalu- 
ation of the lumped circuit parameters which character- 
ize the slot and the external region. 

The susceptances which characterize the slot are in 
each case obtained by a “small-aperture” procedure us- 
ing results which are derived in the Appendix. The con- 
ductances are obtained by the application of a varia- 
tional procedure employing small-aperture quasi-static 
trial fields. 


A. Perturbed Hy, Mode 


The coordinate system and the electric field distribu- 
tion for this case are indicated in Fig. 2. The appropriate 


Fig. 2—Slotted circular waveguide carrying 
the perturbed Hy, mode. 


discontinuity problem is that of the »=1 radial H-type 
mode [in (5) | incident on the slot from within the wave- 
guide. From the solution of this discontinuity problem, 
given in the Appendix, one obtains 


$0 M, 

f To x E,ad@¢ an reel Dy ay a 0), (6) 
—¢%9 Jwe 

where 


E, is the tangential electric field in the slot, 
I,” is the discontinuity in the modal current of 
the n=1 radial H-type mode at the slot, 
hy''(a, 0) is the value of the n=1 H-type mode func- 
tion at the center of the slot, and 


Pa. i k? 
Dh bP) tn (2 ag) 


izability of the slot for the axial component 
of the magnetic field.® 


M is the magnetic polar- 


2 


6 The axial magnetic polarizability of an infinite slit in an infinite 
conducting plane is known for the case k,=0 to be 


ees 
M.=— ————- 
2 In (8/kyd) 


(See N. Marcuvitz, “Waveguide Circuit Theory: Coupling of Wave- 
guides by Small Apertures,” Polytechnic Inst. of Brooklyn, Brooklyn, 
N. Y., Rept. R-157-47, PIB-106, p. 42, eq. (26c); 1947. It is seen that 
this polarizability is similar in form to that for the cylinder, but, in 
contrast to the cylinder result, it is dependent on frequency. 
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The modal voltage of the »=1 mode at the slot is 
given by 


go 
i J hy'*(a, ¢) to X Ex(o)add 


—0 


0 
ro x E,ad¢. (7) 
—$o 


~ hy'’*(a, 0): 


The substitution of (6) into (7) then yields 


Vi" = —I)"hy'"(a, 0) -hy''*(a, 0). (8) 


Jue 


When the polarizability value and the appropriate mode 
functions from (5) are substituted into (8), the un- 
normalized slot susceptance is obtained directly as 


ut De 2 
Zz 


= jB = jue (— : Jaina), 


k? 


Ty 


Vee 


while the normalized susceptance becomes 


B a 
B’ = —— = 2xaln (2/0), 


& (10) 
0 
where 
K 
ys" pel BA, 
wpa 
and 
k= Vk? — k,. 


The susceptance is normalized to the characteristic ad- 
mittance of the angularly symmetric (~=0) H-type 
mode rather than the »=1 mode. This is a matter of 
convenience, since any normalization is satisfactory pro- 
vided all admittances are normalized in the same way. 

The conductance portion of the equivalent circuit 
which terminates the radial transmission line in this case 
can be calculated from 


R 
G= ) 
| Vy" |2 


(11) 


where P is the radial real power flow per unit length and 
Vi" is the mode voltage of the »=1 radial mode at the 
slot. The power P is given by 


$0 0 

P= Re {i E X H*-roadd = Re f E,H,*add. (12) 
—oo —$o 

The magnetic field H, can be represented by 


ie} 
aA tt 
Dd, Ty: Hien ) 


n=0 


H. (13) 


with 
—S 
Yate Vat 


a (14a) 


aoe 
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and 


0 0 
| aos sg f Ae x Eado a hin’ Egadd, (14b) 
—$9 —¢$9 


sire 
where Y,’’ is the admittance seen by the mth mode look- 
ing in the direction of increasing r from the slot. When 
(13) and (14) are substituted into (12), the expression 
(11) for the conductance may be written as 


?0 $0 
[a6 [a6 F9(a, 6)6"(a; 6, 6) Fela, #) 
Ge —o —$o 15) 
0 2 ‘ ( 
hz Eg(a, o)d 
—¢o 


where 


ao = 
G"' (a; 4, ¢’) =! Re ye Vn! (a) hen! (6) Ren’ ($’). 


n=0 


This expression for the conductance can readily be 
shown to be stationary with respect to the variation of 
the aperture field Ey(a, @) about its correct value. 


> 
The expressions for the admittances Y,(a) are given 
in the Appendix as’ 


=> xk Hy, (ka) 
Vota) = c (16) 

jopa H,,'(xa) 

Therefore, one finds 
= 1 —- — 
Re Y,.’(@) = i [Y,/"(a) + Y,/*(a)] 
Z 1 

(17) 


4 Twa” | A,’ (xa) |? 
When the “static” slot electric field, 
1 
~ V1 = (6/60)? 


derived from the integral equation solution presented in 
the Appendix as (53), is employed as a trial field in (15) 
and the integration is performed, the unnormalized con- 
ductance becomes 


Es 


big En J o(ngo) 1]? 
C= read 1a) al roe esa 
where 
rea We n= 0 
= 2, n> 0, 


and J» is the Bessel function of the first kind of order 
zero. When the conductance is normalized in the same 


7 1,’ represents the first derivative of the Hankel function with 
respect to its argument. 
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manner as was the susceptance in (10), one obtains 
G dae? En Jo(ndo) |? 
Gee St E 2). (19) 
Yo xa no | Hn" (xa) |2L So(oo) 


The complete transverse network for this case is 
shown in Fig. 3. 


r=0 be— aq ———_w| =a 


Fig. 3—Transverse equivalent network of slotted circular 
waveguide for the n=1 radial H-type mode. 


Fig. 4—Slotted circular waveguide carrying 
the perturbed Eo, mode. 


B. Perturbed Ey, Mode 


The geometry and magnetic field for this case are 
shown in Fig. 4. The relevant discontinuity problem is 
that of the angularly symmetric (n=0) E-type radial 
mode incident on the slot from within the waveguide. 
The determination of the appropriate circuit parameters 
parallels that in the preceding case. From the results 
given in the Appendix, one obtains 


0 
if Fe EV(a; Wadd — HM he a (Oe 
—oo 


(20) 


where 


Ty’ is the discontinuity in the modal current of 
the n=0 radial E-type mode at the slot, 
hy’ (a, 0) is the value of the »=0 E-type mode func- 
tion at the center of the slot, 


Rifkin 
Me==(=) (2ago)? is the magnetic polariz- 


ability of the slot for the @ component of 
magnetic field,’ and where the symbols ro 
and EF, have the same meaning as in (6). 


8 The polarizability in this case becomes identical with that for an 
infinite slit in an infinite plane, when the width d of the slit in the 
plane is replaced by the arclength 2a¢p of the slit in the cylinder. (See 
C. G. Montgomery, R. H. Dicke, and E. M. Purcell, “Principles of 
Microwave Circuits,” M.I.T. Rad. Lab. Ser., McGraw-Hill Book Co., 
Inc., New York, N. Y., vol. 8, p. 178, Table 6.1; 1948.) 
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As before, the mode voitage at the slot can be ex- 
pressed approximately as 


$0 
Vo = hy'* (a, 0) Z ro X E,(a, )add. (21) 
—%$o 
The unnormalized susceptance is then given by 
To! ooh 
Fogel re 
Vo wuM sho (a, 0) “ho (a, 0) 
=e 
Ejeet iS, (22) 
K 2 
arn (=) po 
and the normalized susceptance by 
B 8 
pi =—= (23) 
Yo Kado” 


A variational expression for the conductance in this 
case may be formulated in a manner similar to that for 
the H-type mode. The resulting expression is 


$0 0 
i dp | dd'E.(a, $)G'(a; $, ¢’) E.(a, ¢’) 
—$o —$0 
G. = ie 3 J (24) 
hoo E.(a, ¢) do 


—$0 


where 
oe 
G'(a; ¢, ¢’) = Re yt Vn’ (a)hon (a, &)hon'(a, 6’). 
n=0 


— 
The admittance Y,’/(a) is the dual of the impedance 
— 


Z,/'(a) given previously and is 


— Jed H,, (ka) | 
Y,,/(a) = ? (25) 
K H,,° (xa) 
so that 
= Ieee ay 
Re ¥q/(a) = = [¥n'(a) + Ya'*(a)| 
2we 1 
(26) 


4 1K? | H, (ka) |? 


When the mode functions from (4) and the “static” 
trial electric field 


E, = V1 — (¢/¢0)?, 


given by (73) of the Appendix are substituted into (24), 
the integration can be performed to yield the following 
expression for the unnormalized conductance: 


Ge Swe Pe En Bex ). 


TK? n=0 | A, (xa) |2 noo 


(27) 
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where 
a) 


m > 0; 


é, = 1, 
= 2, 


and J; is the Bessel function of the first kind of order 
unity. The normalized conductance is then given by 


G = 8 Ss En (ay. (28) 


GQ = — = 
Yo = TK n=0 | H, (ka) |2 noo 


The complete transverse network for this case is shown 
ity Figs. 


O O 
r=0 = c= 


Fig. 5—Transverse equivalent network of slotted circular 
waveguide for the n=0 E-type radial mode. 


It should be noted that the above expressions for the 
susceptances given by (9) and (22) are strictly valid 
only for small-aperture angles (sin ¢o ~@o). However, as 
shown in the next section, the small-aperture approxima- 
tion can be pushed somewhat beyond the above re- 
striction. The conductance expressions, while not lim- 
ited by small-aperture considerations, are practical only 
when the associated infinite series converge rapidly. The 
rapidity of convergence depends upon the magnitude 
of xa. For small values of xa the series converge rapidly, 
but for larger values the convergence is slower. In the 
two cases treated here the values of xa are of the order 
of 1.841 and 2.405, respectively. In this range the series 
converge fairly rapidly; in the calculations of the next 
section four terms were satisfactory for the H-type mode 
and five terms for the E-type mode. 


IV. PERTURBATION SOLUTIONS AND COMPARISON 
WITH MEASURED VALUES 


The propagation constants of the two lowest leaky 
waves in slotted circular waveguide are obtained from 
the solutions of the resonance equations for the two 
transverse equivalent networks derived in the preceding 
section. The networks are shown in Figs. 3 and 5. The 
resonance equations for these H-type and E-type waves 
are, respectively, 


fa 
Z1""(a) 4 R XY! = 0 (29) 
Zy'(@) eee 

and 
re 
Zit) en eae 
ARR M10 (30) 


Zo (a) 


IS. 5 
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where, as shown in the Appendix, 


oe 
Z1'(a) 


Eo _ Ji (ka) 
Pagan eee Te) 
and 
Zila): © Fol) 
Zo (a) 4 Jo’ (xa) ; 


The normalized values of resistance and reactance R’ 
and X’ are those corresponding to the slot discontinuity 
and are computed from the conductances and suscept- 
ances given in (10), (19), (23) and (28) by means of 
/ G’ / B’ 
Gt Be 


Ka 


It should be recalled that the impedances for the n=1 
mode are normalized with respect to the »=0 mode 
characteristic impedance for convenience, as discussed 
in Section ITI. 

As in the case of rectangular leaky waveguides,! the 

resonance equations are complex transcendental equa- 
tions whose exact solutions require either computing 
machines or arduous numerical computations. A per- 
turbation technique is therefore employed to obtain sim- 
ple essentially-closed-form solutions for the propagation 
constants. A perturbation procedure is also appropriate 
because most cases of practical interest are concerned 
with slow leakage. 
- The unperturbed structure for these cases is taken to 
be a completely closed circular waveguide. The unper- 
turbed transverse wavenumbers are thus the wave- 
numbers of the lowest H and E modes in the round 
guide and are koa=1.841 and koa =2.405, respectively. 
The perturbed wavenumbers can then be written in the 
form 


K = Ko + Ax. 


When the resonance equations (29) and (30) are ex- 
panded about their corresponding unperturbed wave- 
numbers ko, and the terms containing the derivatives of 
the terminating impedances are neglected, the perturba- 
tion Ax for the H-type wave is obtained as 


(x oa)” 


ie “( (XG) + 5Ro), (31) 
a@ \(xea}*=— 1 


with 
koa = 1.841, 


and for the E-type wave as 


An = ~ GR) — X’(ko)), (32) 


with 
Ko@ = 2.405. 
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The leaky wave propagation constants are then given by 


he = B — ja = JSR? — 2 = VSR? = ko? — 2koAk 


et ee KoAK 
ay ras c(t = a“). 
k? = Ko" 


Thus, we may write for the H-type wave 


B r 


k Ng 


R 


eae (=) ( sae ) Ra) (34) 
2r dX \a/ \(koa)? — 1 
with , 
ve sa, 
a/ kh? — Ko? 
and for the E-type wave 
r r Ko@ [Ago \? 
a) 
on w= (AY Re), (36) 
with 
we 2a : 
Jk? — ko? 


In Figs. 6-9 theoretical curves computed from (33)— 
(36) are compared with measurements published by 
Harrington.2 Since the susceptance expressions em- 
ployed in the resonance relations are based on small- 
aperture solutions, it is to be expected that the smaller 
angle results in Figs. 6-9 agree well with the measured 
values, while those for the larger angles show a dis- 
crepancy. 

Harrington’s theoretical results? are not shown here, 
but his analysis is similar to Rumsey’s work on rectan- 
gular waveguides? in which a variational procedure is 
used that is equivalent to the determination of the 
equivalent circuit parameters of the slot. The trial 
fields used by Harrington differ from those used in this 
paper for the evaluation of the radiation conductances 
in that the fields he uses do not satisfy appropriate edge 
conditions. In calculating the propagation constants, 
Harrington has also employed a perturbation technique, 
but he retains terms of the order of (Ax)?. Consequently, 
his results for the attenuation constant in the case of 
wide slots are superior to those computed from (34) 
and (36), but at the cost of a substantial increase in the 
complexity and difficulty of actual computation such 
that no explicit expressions for the attenuation con- 
stants can be obtained. 


°V. H. Rumsey, “Traveling wave slot antennas,” J. Appl. Phys., 
vol. 24, pp. 1358-1365; November, 1953. 
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Fig. 6—Guide wavelength vs guide radius for H-type 
leaky wave in circular guide. 
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APPENDIX 


DIFFRACTION OF A RADIAL WAVE BY A SLOTTED 
CIRCULAR WAVEGUIDE 


The equivalent circuit representations of the slotted 
circular waveguide, which are employed in the trans- 
verse resonance calculation for the leaky wave propaga- 
tion constants, require the solutions of two diffraction 
problems, one for each type of leaky wave. These two 
correspond to the n=1 H-type radial mode and n=0 
E-type radial mode incident on the slot from within the 
circular waveguide. 

Each of the two diffraction problems is stated in terms 
of an integral equation for the electric field in the slot 
aperture. The integral equations are solved in the static 
and small-aperture limits, and expressions are derived 
for the aperture fields and the slot polarizabilities. 


A. H-Type Radial Mode Incident (n=1) 


The radial H-type mode incident from within the cir- 
cular waveguide possesses an angular dependence cos ¢ 
and an axial dependence exp(—jk.z). The transverse 
(to r) fields of this incident mode are given by!® 


Ezinc(1, ¢) = Vine’ (r)er’’(r, ?) 


Hiine(r, ) = ling (ry hi (r, ¢), (37) 


where the mode functions e:’’ and Ay’ are those given 


in (5) for n=1. If the region within the guide (r <a) is 
designated as region 1 and that outside (r >a) as region 
2, the transverse (to 7) magnetic fields in the two regions 
can be represented in terms of the complete set of H- 
type vector mode functions as 


i) 


Hi (r, 6) = Dy In” (nr) bn’ (1, 9) 


n=0 


HO, 6) = Stn", 6). (38) 


n=0 
The imposition of the continuity requirement on the 
tangential magnetic field at the aperture then leads to 
an integral equation for the transverse electric field in 
the aperture as follows: Continuity of the tangential 
magnetic field at the aperture requires that 


H,(a, ¢) e.. A, (a, ¢). 


When these fields are expanded in terms of the radial 
- modes, 


> I," (a) An(a, $) = 3 I," (a)hn(a, $) 


n=0 n=0 


is obtained, where J,‘’’(a) and I,’’(a) are the modal 
currents at r=a~ and r=at, respectively. This expres- 


10 The exponential z dependence is the same for all field com- 
ponents and is suppressed. 
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sion can now be rearranged slightly to yield 
T"(a)hy'"(a, ¢) = OT," (a) A, (a, 4), (39) 
n=0 ‘ 


where 

Iy"(a) = 1"(a) — 11a) 

Le MOL G) 2 tae (a), 
and the prime on the summation indicates the omission 
of the m=1 term. 

For all except the »=1 mode, 
Smeal 
L,"(a) = Vx""(a)Vn""(a) 


Tn"(a) = — Yn'"(a)V (a). (40) 
— a 

The admittances Y,’’ and Y,’’ are the admittances 

“seen” by the nth H-type mode at r=a, looking in the 

direction of increasing ry and decreasing r, respectively. 

These two admittances can be shown to be 


> kK Ay (xa) 

Y,(a) == a 
jopa H,,)’ (xa) 

— fol, ENCE 

Y, (a) ve . Sy 


jopa J,’ (Ka) : 


where k= /2?—k2. It should be noted that these are 
not the radial line characteristic admittances since on a 
radial transmission line, unlike a uniform line, the actual 
current-voltage ratio is not equal to the characteristic 
admittance even under matched conditions. 

The modal voltages V,,’’(a) are given by the trans- 
forms of the transverse (to r) electric field as 


$0 
Vn’ (a) = i ha eee x ro E,(a, ¢’)adg’ 
—$o 


%0 
= Iten'' (o') Eg(a, $')adg’. 
—$0 


(41) 
By combining (39), (40), and (41), one obtains the 
integral equation 


I,'’(a) ha"’(a, ¢) 


$0 


<—> 
= Do! Vn! (ayhen’ (a, b) hin’ (a, ¢') Eg(a, ¢’)adg’, (42) 
—oo n=0 
where 


5 ed 


— — 
Vee = Vert + Vee 
The integral equation (42) may be rewritten as 
= Tate rp 
[T1/"(a) + Y1'"(a) Vi" (a) |hea!" (a, 6) — Vo" (a) Vo" (a) ho’ (a, b) 


Ch), APES 
= DY Vn!" (a)hen' (a, b) hen’ (a, 6’) Es(a, o’)adg’. (43) 


—o n=1 
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In order to recast the above integral equation into a 
form permitting diagonalization of the kernel a set of 
“static” admittances are introduced as follows: 


o 2K? 
Yi (a4) =| — = , 
JON 
so that 
fast —> 
V,/’(a) — Vns''(a) > 0 for n> ka. 


The integral equation (43) can also be rewritten, em- 
ploying a compact notation patterned after that used 
in the “Waveguide Handbook,”" in the form 


foo) 


Do Lnhen''(a, ¢) 


m=0 
go 2 
= De ae WG), b) hen’ (a, g \Ey(a, oy) ‘\ad¢’, (44) 
—$o n=1 
where 
a <> 
ie =, Fe Vo Ve; m= 0, 
e = 
Ths — Tene + Vie Vas m= 1; 
a cod = 
Le = (ia ac Pay Vee he wale 


In view of the linearity of (44) the electric field in the 
aperture can be represented in terms of partial fields as 


Spon ICE 


m=0 


When the slot aperture is narrow (sin ¢y9~qo), the 
set of integral equations (44) can be simplified by ap- 
proximating the effective exciting currents 7, in the 
following manner: 


thal jin = al 


> 


yO; m # 1. 


The justification for this approximation is the following. 
—_. 


It is to be noted that the true Y, values are complex; in 
fact, the few lowest ones, which contribute significantly 
to the far-field radiation pattern, are predominantly real 
in character. Alternatively stated, the radiated power is 
not carried by the »=1 mode alone. As a result, the 
admittance derived here for the slot actually should 
contain a conductive contribution. The susceptance 
value of a narrow slit, on the other hand, consists largely 
of contributions from the high modes, which are essen- 
tially purely reactive and for which Y,= Yp, is an excel- 
lent approximation. Furthermore, while it is true that 
for the lowest modes the difference between Y, and VY. 


ns 


is largest, these modes do not have much reactive con- - 


1 Op. cit., Marcuvitz, “Waveguide Handbook,” Section 3.5. 
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tent and will not contribute significantly to the sus- 
ceptance of a narrow slit. For these reasons, the pro- 
cedure below will result in a purely susceptive slot ad- 
mittance which can be expected to be a reasonably accu- 
rate approximation to the slot susceptance. The radia- — 
tion conductance value obtained on this basis, however, 
would be in serious error. A variational procedure is 
therefore employed, in Section III, to evaluate the radi- 
ation conductance values in a separate and accurate 
fashion. ~ 

With the approximations indicated above, tte infinite 
set of integral equations (44) reduces to the single 
integral equation 


h(a, ?) 
go 2 


= ba Yo! (a) en!’ (a, &) Men" (a, $')Eg1(')adg’ . 


—$o n=1 


(45) 


When the mode functions and static admittances are in- 
serted into (45), the integrand is seen to be an even 
function of ¢, and the equation assumes the form 


cos p 

Ve 
In order to diagonalize the kernel of (46), the following 
identity!” is invoked: 


(46) 


4a 7% % cos nd cos n¢’ 
f > ————— £1 (¢') d¢’. 


Jor 


n=1 


2 cos nd cos nq’ 


-2>5 = In 2|cos ¢ — cos ¢’|, (47) 
n=l 
and a change of variable is made according to 
cos ¢ = acosé+ 8B. 
Since 0<@<do, we have that 0<6<z7, and 
a = sin? ¢/2, B-="cos* day. (48) 
Eq. (46) may then be written as 
acos@-+ 6 
2x2 a a A 
BS a wall In 2a | cos@ — cos 6’ | 851(6)d0’, (49) 
where 


/ 


Since the aperture is assumed to be small, 
a~ 07/4, B aoaal 


With the above approximations and the use of the 
identity, (48) and (49) may finally be written as 


2 64 cos 28 cos n6’7) , 
1 of [tna 2 EE | Bcoyaet. (50) 
JOM VT n=% 


2 Cf., ibid., p. 147. 
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If 641 is now represented as a Fourier cosine series, 


841(6’) = >> a, cos 6’, 


(51) 
n=0 
it follows immediately from (50) that 
WANS 1 
one i(=) Gees thet poe (52) 
K 4wear/ In (2/0) 


and 
a, = 0, i> OF 


It follows from (51) and (52) that for a narrow slot 
(sin do~q@o) the electric field in the aperture is given 
approximately by 

di 1 1 


ea a 
rll ) Qav/r In (2/60) W¢0? — ¢? 


E3(¢) ~ (53) 


It is also of interest to note that the solution for the field 
satisfies the edge condition. 

The magnetic polarizability of the slot for this case 
can be defined as follows: 


$0 $0 
f ro X E.(¢)add = Zo E4(¢)add 


—¢o —$o 


dee 


ts 
ey ti 
jwe\ kK 2 In (2/¢0) 
Eq. (54) can then be written in the following form: 
$0 1 
[tx Bio)adg = =— MaN'hy"(a,0), (55) 
—$0 JWE 
where 
hy'"(a, 0) = 20 1/7; 
upon comparison with (54) we find that 
RN 1 
M. = — ( ) (56) 
Js K In (2/0) 


is the axial magnetic polarizability of the slot in the 
cylinder. 


B. E-Type Radial Mode Incident (n=0) 


In this case the transverse (to 7) fields of the incident 
mode are given by 


E, ACE ¢) = Vine (r)eo (7, ¢), 


Hi; ine(7, ¢) = Tine’ (r) ho (r, ¢), (57) 


where the mode functions eo’ and Ay’ are given in (4) 
with n=0. 

When the external and internal tangential magnetic 
fields at the slit are equated as in the preceding case the 
following integral equation is obtained: 
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To’ hgo' (a, ¢) 
PO! 10% ey 
oom D2 Yn! (a) hgn! (a, hon’ (a, ') E-(o')adg’, (58) 
—o n=1 


where 
Ty = 1,00 =O 
is the discontinuity in the current of the »=0 mode at 
¥=@, and 
oO — = 
Y,/(@)-= V, (@) +, @ 
is the dual of that in the first section of this Appendix. 


— — 
The admittances Y,’(a) and Y,'(a) are therefore given 


by 


Fy Sete Hale 
« H,)(Ka) 
— iwead J,/ (ka 
Y,/(a) = —2 eer (59) 
« Ja(ka) 
The “static” admittances for this case are 
<= JwEn 
Ya (G) = Vrs (oe) = — pct (60) 
K 
Integral equation (58) can then be written as 
, Lnltom' (a, ¢) 
m=0 
oo © gy 
=— DS Vns’ (a) hon’ (a, &) hon’ (a, $’) Ez(¢’)adg’, (61) 
—$9 n=1 
where 
f,=I, . m= 0, 


P «> Oo 
Lm = (Yn (2) - Va Ves m> 0. 

As before, the aperture field can be represented in terms 
of a set of partial fields, 


EA) ee yD; Im&em()- (62) 


m=0 


If higher exciting currents /,, are again neglected in com- 
parison with /)=Jo', the set of integral equations re- 
duces to 

oo 2 


S WE 
tess inf — ={ >= n cos nd cos ng’&:0(6')do’. (63) 
Tv K? —0 n=0' 


The integral equation (63) may now be solved in the 
following manner: First, we recognize that 


go 2 
>> 2 cos no cos n¢'&.0(¢') do’ 
—oo n=1 
go © 
Sy > n cos ng cos n'&.0(¢')d¢’. 
0 n=1 
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If the above integral is now integrated by parts, we have 


go © 
2 Se n cos nd cos nd’&0(¢') dd’ 


0 n=1 


oo $0 
= 2 DS cos n@ sin na’) 


mr=1 0 


COMME (hee 
+2f sie 
0 


dq’ n=1 n 


aed 
£08 00 008 OTe ah anne 
dq’ 
Since &20(¢) vanishes at d=», 


go © 


Ds n cos nd cos nd’&0(¢’)d¢’ 


—oo n=1 
Gh ss 
s 2f >» 
0 


dq’ n=1 


cos nd cos nd’ d 


— &20(¢’) dg’. 
ig! o(b')dd 


Employing the identity (47), we see that 


d > cos n@ cos ng’ = 1 


dq’ n=1 n 


sin ¢’ 


2 cos¢ — cos ¢’ 
Eq. (63) can now be written as” 
sin ¢’ 


BOF $0 d 
eye i &20($’)de’. (64) 
2J, cos¢—cos¢’ dd’ 


Ts 3K 


The following change of variable is now made: 
COs @ —==a cos 04-"8, 
where, for 0<¢<do, we have, as before, 0<@<m with 


(6 =sC0s" $o/ 2, a = sin? do/2, 
so that 
T Q’ 
te iH SOT: (65) 
9 cos@— cosé’ 
where 
(0') i, ThE neem 
=— —- — sin &’ — &, : 
§ J mae de’ 0 
The function g(6’) may be represented as 
g(0’) = >, a, cos v6". (67) 
v=0) 
Then, since! 
* cos v6’ sin v0 
f ——— df’ = — r ——» (68) 
0 cos@é— cos@’ sin 6 
(65) becomes 
2 sin v0 
Lo = a, i (69) 


a sin 6 


13 This integral exists in the meaning of Cauchy’s principal value. 

14 W. Magnus and F. Oberhettinger, “Special Functions of Mathe- 
matical Physics,” Chelsea Publishing Co., New York, N. Y., p. 141; 
1949. 
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Since, from (66), 
g(0) = 0, 


we note from (67) that 


> a, = a+ >) @ = 0. (70) 
v=0 v=0 ~ 
From (69) and (70), therefore, 
ao = — a = 1/1, a, = Uw 
and as a result, we find 
"8 we d&, 1 
g(0) = —jA/— a eee Aap —(1—cos@). (71) 


do T 


Since the aperture is assumed to be narrow, sin $0 ~¢o, 
it follows from (66) and (71) that 


hoe tee Ged ¢ 
do oo we V/8r Vbo? — ¢? 


When the differential equation (72) is integrated, the 
result is 


(72) 


2 


K 1 Fen Be ch 
ig = —-—-j-—= 2— ¢?4+ K; 
o(¢) i) nee V0 ¢? + 


however, since &20(@0) =0, K=0. In view of (62), the 
aperture electric field is therefore given by 


2 , 


PO tke! Ti SOkee 
E(¢) =f = Vb ee 


we VW 8r So) 


The magnetic polarizability of the slot for this case 
is then obtained by writing 


%0 0 
fi rut El@plads =) he" | suealage 


—¢0 —o 


jou [ K\? Tado” 
=_- — — {| — Tire 
» (=) ( 2 ) 


= jwuM 4Io'ho (a, 0), 


where 
hy’ (a, 0) = ol/ar/ 2; 


as a result, we have that 


x fe NA 
Me = 35 ({) Gaon 


is the magnetic polarizability for magnetic field excita- 
tion transverse to the slot. 


(74) 
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A Spiral-Grating Array* 


J. R. DONNELLAN} anv R. T. CLOSEt 


Summary—A problem associated with a two-dimensional spiral 
doublet array is concerned with the fact that use of the spiral doublet 
as an element to obtain linear polarization leads to spacing difficul- 
ties, since the distance between elements is always twice as much in 
one dimension as in the other. Thus, the allowable scan is seriously 
limited in one principal plane. A method is proposed for obtaining 
linear polarization from a single spiral antenna by using it in com- 
_ bination with a polarization grating and a ground plane. Experimental 
results were obtained by combining eight of these elements in a one- 
dimensional array. For the broadside position of the main beam, the 
field was found to be linear, with the cross-polarized component of 
the main beam being 40 db down from the principal polarization. As 
the beam was scanned out to +45°, the linearity decreased quite 
rapidly (to approximately 10 db at 45°) and the plane of the linear 
polarization was rotated; sidelobes were higher than predicted. Use 
' of this array rather than an array of doublets would involve a com- 
promise concerning specifications, dependent upon the particular ap- 
plication of the array. 


INTRODUCTION 
\ SUCCESSFUL scanning array employing flat 


two-wire Archimedean spirals has already been 

developed at the U. S. Naval Research Labora- 
tory.' In order to obtain a linearly-polarized beam, this 
system relied upon the properties of the spiral doublet: 
two equally excited spirals of opposite sense placed side 
by side. However, use of the doublet as an element 
leads to a spacing problem in a two-dimensional array: 
the spacing between elements is always twice as much 
in one dimension as in the other. Thus, the maximum 
attainable angle of scan (i.e., the maximum angle to 
which the beam can be scanned without the introduc- 
tion of an extraneous main lobe) is seriously limited in 
one principal plane. 

Several proposals that appear to remedy this situation 
have been considered. For example, the use of the 
doublet configuration can be retained if the size of the 
- individual spiral elements can be reduced, or if, by some 
arrangement of elements, the spacing does not become 
excessive. The system described in this report is the 
result of an attempt to eliminate the need for the doub- 
let in order to produce linear polarization from an array 
of spirals. It makes use of a polarization grating of paral- 
lel wires? to produce the requisite linear polarization. 


* Received by the PGAP, July 15, 1960; revised manuscript re- 
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THEORETICAL CONSIDERATIONS 


The following discussion demonstrates a method of 
obtaining linear polarization from a single spiral an- 
tenna by using it in combination with a polarization 
grating and a ground plane. The grating consists of 
closely spaced wires stretched on a frame, the spacing 
between the wires being much less than a wavelength 
and much greater than the wire diameter. Consider now 
a vertical grating situated a quarter-wavelength, at the 
frequency of operation, above a ground plane and a 
spiral situated at an arbitrary distance d in front of the 
grating (Fig. 1). When the spiral is properly excited at 
its terminals, it will radiate equal amounts of circularly- 
polarized radiation (of opposite senses) in both the for- 
ward and backward directions. The sense of circular 
polarization is that observed when the radiation is 
viewed in the direction of propagation. 
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Fig. 1—Spiral-grating configuration. 


Assume, for example, that the spiral radiates left- 
hand circular polarization in the +Z direction and 
right-hand circular polarization in the —Z direction, 
Consider the radiation (right-hand) which travels to- 
ward the grating. The circular polarization can be re- 
garded as the vector sum of horizontal and vertical 
components in phase quadrature. The vertical (parallel 
to the X axis) component travels a distance d, is re- 
flected at the grating with a 180° phase change, and re- 
turns to the spiral, again traveling a distance d, so that 
the “effective path length” traveled is 2d+)/2. The 
horizontal (parallel to the Y axis) component passes 
through the grating after traveling a distance d, then 
travels a distance \/4, is reflected from the ground 
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plane with a 180° phase change, again travels a distance 
d/4, and finally a distance d, returning to the spiral; 
the “effective path length” traveled is 2d+A. Thus, the 
phase angle between the horizontal and vertical com- 
ponents has changed from +90° to —90°; z.e., the com- 
bination of grating and ground plane serves to retain 
the rotational sense of the radiation incident upon it. 
This reflected radiation of one sense will combine with 
the direct radiation, which is of the other sense, to pro- 
duce linear polarization. The orientation of the plane 
of linear polarization is a function of d for a fixed grating 
position. For a vertical grating, if d=h/4, vertical polar- 
ization will be obtained (if d=\/8, 45° polarization is 
obtained, etc.); for a horizontal grating (with d=X/4), 
horizontal polarization is obtained. It can also be shown 
that, with the spiral fixed at a distance \/2 above the 
ground plane, the spacing between grating and ground 
plane can be varied to produce linear polarization on- 
axis and at one other angle off-axis (this was actually 
done in the experimental work). 

Thus, with the grating in a fixed position and the 


spiral a fixed distance above the grating, the spiral may 


be rotated in order to change the phase of its far field 
while maintaining linear polarization. That is, rotation 
of the spiral through an angle @ advances the phase of 
both the right-hand and left-hand circular polarization 
by an amount 6; therefore, the plane of the linear polar- 
ization is maintained, while the phase of the far field is 
changed by an amount 6. A scanning array could thus be 
composed of many of these elements, and equal spacing 
between elements would be maintained in both dimen- 
sions. 


EXPERIMENTAL CONSIDERATIONS 


Preliminary experiments were carried out with a sin- 
gle spiral and grating in which the spacings between 
spiral and grid, and between grid and ground plane, 
could be varied. These experiments verified that this 
system would indeed provide linear polarization, while 
maintaining the phasing properties of the spiral. It was 
found that, off axis, the plane of linear polarization ro- 
tated from that of the broadside position. This was due 
to the fact that the effective spacing between grid and 
ground plane for an off-axis position is greater than the 
spacing for the broadside position; thus an additional 
relative phase change between the right-hand and left- 
hand circular polarizations is produced and the plane of 
linear polarization is rotated. However, it was decided 
that an array would nevertheless be built since rotation 
of the plane of polarization might not be a serious dis- 
advantage, depending on the particular application. It 
was experimentally determined that, in order to obtain 
suitable linearity for both the broadside and off-axis 
positions (out to approximately +45°), the spacing be- 
tween spiral and grating should be approximately \/8 
and the spacing between grating and ground plane, ap- 
proximately 3/8. 
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A photograph of the eight-spiral array, operated at a 
frequency of 1420 Mc, is shown in Fig. 2. The array 
consisted of one row of spirals (all of the same sense) 
placed in front of a horizontal grid structure. The feed- 


board was a microstrip ring network type, in which the © 


desired aperture illumination was obtained by varying 
the lengths of line between the ring networks. A diagram 
of this feedboard is shown in Fig. 3. It was designed to 
give a 25-db Tchebycheff illumination working into 
matched loads. The grating was constructed of no. 24 
(1/32-inch diameter) tinned copper wire with a $-inch 
spacing between wires. The lateral spacing between 
spirals was \/2. 


Fig. 3—Microstrip ring network feedboard. 


An additional feature of this array was the use of 
“plug-in” spirals; 7.e., two small male conductors 
(UG-260-U) were attached to the input terminals of 
each spiral, while two small female connectors (UG-260- 
U) were attached to each two-wire line. The spirals 
could then be easily removed from, or inserted in, the 
array by simply pulling them out or plugging them back 
in. In the eight-doublet array of Donnellan,! the spirals 
were soldered to the two-wire lines, which made removal 
or replacement of the spirals a rather lengthy task. 

The array was phased for the broadside position by 
using it as the receiver, while a pyramidal horn was used 
as the transmitter. The distance between transmitter 
and receiver was approximately 120 feet (about 5(2D?/)) 
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at 1420 Mc) and the horn could be rotated to provide 
any arbitrary linear polarization. Two elements at a 
time were connected into the feed system (matched 
loads terminating all other connections) and one spiral 
was rotated with respect to the other until both were in 
phase (the technique is described more completely in 
Donnellan.' This process was continued for all elements 
taken two at a time until all were in phase. A second row 
of elements was actually inserted several times and 
measurements were taken. This series of measurements 
indicated that two-dimensional arrays without a spacing 
problem were feasible (with certain limitations, as will 
be seen from the following antenna patterns). Patterns 
of single and double rows were found to be essentially 
the same, so that, in the interest of ease in experimenta- 
tion, pattern measurements were completed with a sin- 
gle row. 

Scanning the beam was accomplished by using again 
the phasing property of the spiral, whereby a mechan- 
ical rotation of the spiral of 1° is equivalent to a change 
in phase of the far field of the spiral of 1°. By setting 
successive angular displacements in each element, the 
beam was scanned in the horizontal plane. The predicted 
angle of scan (@) is related to the phase difference be- 


tween elements with half-wavelength spacing by the 


following relation: 


_ Phase difference (in degrees) = 180 sin @. 


These differences were calculated and the spirals were 
then rotated the proper amount in order to obtain the 
scans shown in Figs. 4 through 9 (pp. 294-295). 


RESULTS 


The experimental program verified the theoretical 
considerations in general, and also pointed out some 
interesting characteristics of the system. For the broad- 
side (0° scan) position of the main beam, the field on 
- axis was found to be linear (approximately 40 db); the 
linearity is taken as the difference, in dbs, between 
the principal polarization of the main beam and the 
cross-polarized component. Fig. 4 shows this broadside 
position of the main beam for the principal polarization 
and for cross-polarization. The plane of polarization was 
found to be rotated 12° clockwise from the predicted 
horizontal position. The sidelobes were higher than the 
—25-db feedboard design which may be attributed to 
the fact that the feedboard was designed for a 25-db 
amplitude distribution looking into matched loads, 
while the introduction of a polarization grating between 
the spirals and the ground plane can conceivably alter 
this amplitude distribution. No attempt was made to 
further improve the sidelobe level, since the principal 
concern was the linearity of the polarization. As the 
beam was scanned out to +45°, it was found that the 
linearity decreased with increasing scan angle (Figs. 
5-9). Concomitant with this, the plane of the linear 
polarization rotated with increasing scan angle. This 
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information is presented with each antenna pattern. 
Another noticeable characteristic of this scanning 
array is the apparent increase in gain as the beam is 
scanned out to +45°. This is due partially to the fact 
that the spiral placed \/8 above the grid and grid placed 
3d/8 above the ground plane had a dip, rather than a 
maximum, in the center of the pattern. Also, no Aqua- 
dag was placed on the outer turns of the spirals (as had 
been done in the eight-doublet array) to improve the 
axial ratio of the individual spirals. (Aquadag serves to 
attenuate any currents which might reach the outside 
edges of the spiral; if these currents were reflected from 
the ends of the spiral, they would radiate in a sense op- 
posite to that of the currents traveling from the input 
terminals and increase the axial ratio. A high axial ratio 
serves to modulate both the amplitude and the phase of 
the radiation, the change in phase no longer being linear 
with spiral rotation and the amplitude no longer remain- 
ing invariant as the spiral rotates.) It might also be 
pointed out that the patterns of the individual spirals | 
predicted better linearity for the large scan angles (not 
less than 14 db) than was actually obtained when the 
spirals were combined in the array. 


CONCLUSIONS 


The principal advantage of the spiral-grating system 
is that it eliminates the need for a spiral-doublet con- 
figuration to obtain linear polarization and thus, in a 
two-dimensional array, allows for equal spacing of 
elements in both principal dimensions. Such a system 
could easily be scanned over a +45° range without the 
introduction of extraneous main lobes (and even further, 
if the spiral pattern permitted). On the basis of the 
experimental work, this can be achieved if one is willing 
to accept linearity of approximately 10 db at the larger 
scan angles. The sidelobes obtained with the spiral- 
grating system were higher than those obtained with the 
eight-doublet array, but no further attempt was made 
to reduce the sidelobe level. Some other considerations 
are that the linearity drops off rather rapidly with in- 
creasing angle of scan and that the plane of polarization 
rotates as the beam is scanned. These factors may or 
may not be disadvantages, depending upon the particu- 
lar application of the system. 

The spiral-doublet array, in contrast with the spiral- 
grating system, produces a field that is linearly polar- 
ized and whose plane of polarization does not rotate with 
the scan angle. If doublets are used as the elements in 
an array, they can conceivably be spaced approximately 
/m (at one frequency) in one dimension, so that a 
+45° scan can easily be obtained. However, the mini- 
mum spacing in the other dimension will be 
2/m(~0.7A), so that the beam can be scanned in this 
plane only + 30° before an extraneous main lobe is intro- 
duced. It seems, therefore, that the choice of either a 
spiral-grating system or a spiral doublet will be dictated 
by the particular use to which the system is to be put. 
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Coupled Surface Waves and Broadside Arrays of End-Fire Antennas* 


J. L. YEN, 


Summary—The concept of coupled surface waves is introduced 
to describe the action of element interaction in a broadside array of 
traveling-wave end-fire antennas. These coupled waves are formu- 
lated in terms of an interaction matrix, whose elements either can be 
determined experimentally, or, for some antenna structures, can be 
computed. Two methods of pattern evaluation based on two well- 
known approaches to end-fire antennas are then presented. The roles 
of interaction in the two methods turn out to be quite different, and 
the method which assumes radiation at the end of the antenna ap- 
pears to be-more favorable. A criterion on the over-all interaction 
effect on array pattern is introduced which indicates under what 
condition coupling can be neglected. Finally, the nature of patterns 
obtainable with different element spacings is examined. 


INTRODUCTION 


‘ 


ROADSIDE arrays of traveling-wave end-fire an- 
B tennas in one form or another such as helices, 

dielectric rods and Yagis have been in use for 
many years. In some instances, the antennas are sepa- 
rated by fairly large distances, while in others the spac- 
ings are small. Despite the variety of existing systems 
one factor common to most is that their proper function- 
ing depends on negligible interaction between adjacent 
elements. The degree of interaction is usually deter- 
mined experimentally from perturbations on pattern 
and impedance produced by coupling, and some qualita- 
tive measure is adopted to indicate its smallness. Only 
in one notable exception! is coupling actually utilized to 
achieve the desired pattern through proper adjustments. 
In view of the many advantages of broadside arrays of 
end-fire antennas, a better understanding of the role 
of element interaction is desirable. The purpose of this 
investigation is to explain the effect of interaction on 
array pattern from the coupled surface waves carried 
by the antenna system and to examine the types of pat- 
terns obtainable with such arrays. The results should be 
of assistance in exploiting element interaction in an- 
tenna systems and to predict when interaction can be 
neglected. 

There are two well-known approximate approaches to 
the radiation of a traveling-wave end-fire antenna; 
one assumes the radiation of a surface wave along the 
antenna, while the other assumes the radiation of a sur- 
face wave at the radiating end of the antenna. To de- 
scribe a broadside array of end-fire elements from either 
viewpoint, the concept of coupled surface waves on the 
antenna structures is introduced. These coupled waves 
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are formulated in terms of an interaction matrix, whose 


elements either can be determined experimentally or, 
for some antenna structures, can be computed. Two 


methods of pattern evaluation based on_the two ap- | 
proaches to end-fire antennas are then presented. The ~ 


roles of interaction are different in the two methods and 
the one which assumes radiation at the end of the an- 
tenna appears more favorable. Finally the nature of 
patterns obtainable with different element spacing is 
examined. 


RADIATION OF A BROADSIDE ARRAY OF 
ENpD-FIRE ANTENNAS 


Before we consider a broadside array of traveling- 
wave end-fire antennas, it is necessary to have a clear 
understanding of the individual elements which make 
up the array. The radiation of a traveling-wave end-fire 
antenna involves the diffraction by the finite antenna 
structure of an incident wave produced by the feed. 
Unfortunately, even the approximate solution of such a 
problem for a very simple geometric configuration in- 
volves considerable mathematical difficulties. A rigorous 
analysis of the general problem is out of the question. 
This leads to the use of crude approximate methods 
which can be more readily applied. The simplest and 
the most commonly employed approximation assumes 
the radiation of a surface wave along the antenna struc- 
ture. Sometimes a direct feed radiation is added for 
higher accuracy. A general discussion of end-fire an- 
tennas of various structural types based on this method 
was given by Reynolds.? A second method capable of 
greater accuracy, but slightly more complicated, as- 
sumes the surface wave to travel without radiation along 
the antenna until it reaches the far end where it is 
diffracted, thereby producing radiation. The diffraction 
can be approximated by the radiation of an aperture 
field at the radiating end of the antenna whose distribu- 
tion is simply that of the surface wave as shown by 
Brown and Spector.* The interference of this radiation 
and the direct feed radiation yields the radiation field of 
the end-fire antenna. Although the two methods assume 
different radiation mechanisms, computation of near 
fields of ideal discrete sources with traveling-wave ex- 
citation by Hyde‘ indicates that in certain situations the 
first method can also be interpreted as an aperture field 
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radiation interfering with the direct radiation from the 
feed. The computed near field distributions greatly re- 
semble the experimental plots obtained by Ehrenspeck 
and Kearns.'*> Zucker and Kay® pointed out that in 


spite of neglecting turbulence in the surface field near 


_ the feed, the first method yields results that are useful 


= 


in engineering design. In particular, it can be readily 


generalized to tapered or modulated antennas.’ Thus, it 
appears that each method possesses certain advantages 


over the other. 


Since both approaches provide useful information on 
traveling-wave end-fire antennas, it is worthwhile to 


_ examine broadside arrays of such elements from both 
_ viewpoints. Fortunately, the surface wave on the an- 


tenna structure is the common basis of the two methods. 
We can, therefore, begin with an investigation of surface 
waves on an antenna system. If the antenna structures 
are well separated, there can be no interaction between 


pe the surface waves on the different structures. The array 
pattern is then given by the product of the element 


- pattern and the space factor determined by the complex 


feed amplitudes and the array configuration, regardless 


; of which approach we use. When they are brought close 
_ together in parallel, the individual surface waves com- 


bine and redistribute into a set of coupled surface wave 


- modes, each with a distinct amplitude distribution 


among the structures and each with a distinct phase 


velocity. If the spacing between antennas is not very 
close, the individual structures can be fed independ- 
ently. The combined action of the feed system, however, 


is to excite in general, not one, but all the coupled 


modes. The coupled modes propagating along the sys- 
tem with different phase velocities result in a continuous 
redistribution of fields about the individual structures in 
the form of spatial beating which can be considered as 
modulation of the waves about the individual struc- 
tures. The manner in which two guided waves interact 
has been discussed at length by Miller.* To calculate the 
radiation pattern using the first approach, we can either 
compute the radiation of the modulated waves on each 
structure first and then combine them to obtain the 
array pattern or, conversely, evaluate the radiation due 


_ to each coupled mode on the antenna system and then 


combine them. When the second approach is to be used, 
the distribution of the surface waves about the indi- 
vidual structures at the radiating end computed from 
the coupled modes can be used. Thus, knowing the 


5 H. W. Ehrenspeck and H. Poehler, “A new method for obtaining 
maximum gain from Yagi antennas,” IRE TRANS. ON ANTENNAS AND 


b, PROPAGATION, vol. AP-7, pp. 379-386; October, 1959. 


6F. J. Zucker and A. L. Kay, “A Surface Wave Antenna Para- 
dox,” presented at the URSI Fall Meeting, Pennsylvania State Uni- 
versity, University Park; October 22-29, 1958. Abstract appeared in 
IRE Trans. ON ANTENNAS AND PROPAGATION, vol. AP-7, p. 111; 
January, 1959. act 
- 7A. S, Thomas and F. J. Zucker, “Radiation from modulated 
surface wave structures—I,” 1957 IRE NATIONAL CONVENTION 
REcorD, pt. 1, pp. 153-160. : : 

8S. E. Miller, “Coupled wave theory and waveguide applica- 
tions,” Bell Sys. Tech. J., vol. 33, pp. 661-719; May, 1954. 
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coupled surface wave modes on the antenna system, the 
array pattern can be derived from the feed amplitudes 
of the individual elements using either method. It turns 
out that the two methods yield different array patterns. 
On the basis of the preliminary experimental data avail- 
able, it appears that the second approach is more accu- 
rate. 

For long antenna elements, because of the length 
available for interaction effect to accumulate, even 
very weak coupling at fairly wide spacing can perturb 
the array pattern considerably. However, knowing the 
coupled surface wave modes we can determine the feed 
amplitudes required for a given array pattern using the 
methods outlined above. At such spacings, if we further 
assume the surface waves are not reflected at the radiat- 
ing end of the antenna system, the impedance of the 
feed on the individual elements will not be affected by 
the interaction. Thus, a feed system can be readily de- 
signed to give the desired feed amplitudes. 

If the spacing between antennas becomes sufficiently 
close, we find that the individual elements can no longer 
be excited independently, 7.e. each element is also ex- 
cited by the feeds attached to the others. When this 
happens, the excitation of the coupled surface waves by 
the feed system becomes complicated. Moreover, this 
may lead to modifications in the feed impedances. To 
avoid these difficulties, such close spacings will be ex- 
cluded from the following discussions. 

Before presenting the two methods of pattern evalua- 
tion and their comparison, we shall begin with a discus- 
sion of the coupled surface waves on an antenna system. 
Without loss in generality, all subsequent discussions 
will be confined to one dimensional arrays with uniform 
spacing for simplicity in analysis. 


COUPLED SURFACE WAVE MODES 


Consider N identical structures, each of infinite length 
and each capable of supporting a surface wave without 
radiation in the absence of the others. In bringing them 
close together in parallel, the N individual modes split 
into N coupled modes. If the spacing between structures 
is large compared with the transverse dimensions of the 
structures, we expect the field distribution of a coupled 
mode in a transverse plane to be approximately given 
by the linear combination of fields centered about the 
individual structures, each with a distribution identical 
to that of the uncoupled surface wave. Let vector? r de- 
note position in a transverse plane and let the reference 
positions of the V congruent structures be located at tn, 
with spacing D as shown in Fig. 1. The transverse elec- 
tric field of the uncoupled mode on the mth structure 
propagating in the +2 direction with the phase constant 


9 We have adopted the convention that all German scripts are 
vectors, all lower case bold face Italics and German characters are 
column matrices with scalar and vector elements respectively. All 
capital bold face characters represent NN matrices while v’ de- 
notes the transpose of matrix. 
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Fig. 1—A one-dimensional uniformly spaced array. 


G can be expressed as 


je — Tne". 


(1) 


Here, as well as in subsequent discussions, the time de- 
pendence e** is omitted. Let (1+A)6 be the phase con- 
stant of a coupled mode. The transverse electric field 
of the mode propagating in the +2 direction can there- 
fore be written in the following matrix form: 


Cs) five Taree 


(2) 


Here, f’ is a row matrix whose elements are the vector 
transverse electric fields in the transverse plane of the 
uncoupled modes on the NV structures as given by 


il = lf ce T1), f(r og T2), or f(v a tw) ], (3) 


and where v is the following column matrix whose ele- 
ments are the complex amplitudes on the JN structures: 


(4) 


UN 


By going all the way to the eigen value problem which 
governs the system, a set of secular equations can be de- 
rived from which the perturbation in phase constant, A, 
and the amplitudes matrix on the structures, v, for a 
coupled mode can be determined. Unfortunately, the 
equations are likely to be integral differential equations 
of a complicated form which cannot be solved using 
present techniques. However, under the above assump- 
tion of weak interaction, we expect the secular equations 
to be of the form 


where 1 is the unit diagonal matrix and K is an NX N 
interaction matrix that describes the coupling effect. 
The off diagonal terms of K should depend on the field 
at one structure due to the wave on another, while the 
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diagonal terms should vanish. Although we cannot de- 
rive the secular equations, e.g., (5) and the interaction 
matrix K for all structural types even at reasonably 


large spacing, by using a multiple scattering approach,’ 


we have been able to show, for an arbitrary structural 
type of small transverse dimensions supporting a dipole- 
type mode, that (5) is valid and that K is of the form, 
0 K( Dye RL) 
K(D) 0 . K(N>ZD) 
K= ; : 
RN. =A DAK GN = Dee 


where K(D), the interaction coefficient, depends only 
on the spacing D in terms of wavelength and the trans- 


verse propagation constant g of the uncoupled mode de-. 


fined as 
2 
q= es —1 (7) 
according to 
K(D) = Say ler) (8) 
: 2 gkD 


The transverse field distribution at a large distance D 
from a structure of arbitrary shape supporting a surface 
wave of phase constant 8 must behave as the asymptotic 
form of a Hankel function with imaginary argument, 
1.€. as 


2 
Se eee nh 
V rgkD : : ©) 
Therefore, as expected, the interaction coefficient (8) is 
proportional to the field at one structure produced by 
the wave on another. For two special structural types, 
namely the dielectric rod!! and the coated cylinder,™ 
more accurate values of the interaction coefficient are 
available. In spite of the fact that the coated cylinder 
supports a circularly symmetric mode, both results agree 
favourably with (8). Therefore, unless more accurate 
values are available, (8) can be safely used. If the struc- 
ture has a large traverse extension, K(D) can be deter- 


May \ 


| 


mined from a simple experiment to be described later. - 


Returning to (5) we see that A is an eigen value and 
v the corresponding eigen function of the characteristic 
equation of K. From (8), observe that K is real, hence 
there are N real eigen values with orthogonal eigen 
functions that can be labelled as A™ and v, respec- 
tively. Real A” implies all the coupled modes are propa- 


10 J. L. Yen, “Multiple scattering, propagation in periodic struc- 
tures and coupled waves.” To be published. 
uM. F. Bracey, A. L. Cullen, E. F. T. Gillespie, and J. A. Stani- 
forth, “Surface wave research in Sheffield,” IRE Trans. oN ANTEN- 
NAS AND PROPAGATION, vol. AP-7, pp. S219-S226; December, 1959. 
_ @ A. A. Meyerhoff, “Interaction between surface wave transmis- 
sion lines,” Proc. IRE, vol. 40, pp. 1061-1068; September, 1952. 
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gating. Furthermore, K has zero diagonal terms, its 
trace vanishes, hence é 


N 


DIN 2) 


r=1 


(10) 


“We shall restrict our attention to situations where all 
the coupled modes are slow waves so that (1+A™)8>k. 
Otherwise, some of the modes become fast waves, and 
the nature of the antenna system will be greatly altered. 
Substituting A and v™ into (2), the transverse electric 
fields of the N coupled modes, e(r, z), are obtained. 

' For later discussion it is convenient to express the 
- transverse electric fields of the coupled modes and their 
relation with the transverse fields of uncoupled modes on 
the individual structures in matrix form. Let the trans- 
verse electric field of the N coupled modes be combined 


in a row matrix, 
2 


v= [eM (x, 2), er, Z), pesrainr eG, zB): (11) 
The N eigen-functions v® can be combined to form an 
_NXN matrix V, according to 
V= [yO, V2 2. ey vy], (12) 
Since v™ are orthogonal, they can be normalized so that 
Vis orthogonal, 1.e., its transpose V’ is also its inverse 
VaVere 17 (13) 
; The propagation of the coupled modes can be expressed 
in a diagonal matrix P(z) defined by 


ei (n)Bz (¢) A @) 
e IA (n)Bz 
i aia ee 67 12? (14) 
@) 0) eA (n) Bz 


Using (2) it is seen that the coupled mode matrix (11) 
can be expressed in terms of the individual mode matrix 
-f according to 


e=TVP): (15) 
The value of the interaction coefficient K(D) as given 
‘in (8) decreases exponentially with increasing distance 
_D between structures. Fig. 2 shows the variation of 
_K(D) vs kD for various values of g. Often it is only 
necessary to consider coupling between immediate 
“neighbors, thus only two off-diagonal terms in K need 
be retained. The characteristic equation, (10), is then 
readily solved, and the resulting A” and V for N=2, 3, 
4, 5 are shown in Fig. 3. Note that the two coupled 
“modes for N=2 correspond to the well-known sym- 
“metric and antisymmetric modes. Also, from these ex- 
amples it can be seen that the larger the number of 
‘structures NV in the system, the larger can the perturba- 
tion in phase constant A” become. As we shall see later, 
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Fig. 3—A® and V for interaction between adjacent elements only. 


for an antenna system of finite length L, the effect of 
coupling on array pattern depends on A™BL. Thus, it 
would seem that the larger the number of structures, JN, 
the more severe will be the overall coupling effect. For- 
tunately, A™ does not spread out over an infinite range 
as N increases, only the density of the NV values of A” 
increases in a finite band. This phenomenon is similar 
to the well-known band theory of conductors where the 
energy levels in a sample of finite number of atoms 
form discrete sets distributed within finite bands. 


q 
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To see the existence of a finite band, consider a one- 
dimensional uniform array with N infinite. The periodic- 
ity in the transverse direction requires the modes to 
satisfy Floquet’s theorem, 1.e., 


ment, the radiation of a surface wave of unit amplitude 
along the first element is simply : 


L p-ikr e 

é : (r) A 
C ela Bedy . 
é 


pea a est Og, 20), (16) r 
Assuming interaction between adjacent structures only, ioe es e-21+A)B(L/2) (7) (9), (21) | 
the characteristic equation, (5), réduces to the single Ri : | 
equation. | 
NO where a form factor C is introduced which may vary 
cos6(”) = ; (17) with 6 slowly in case of structures of large transverse — 
2K(D) dimensions. F\(6) is the well-known sin «/x pattern 


Since all the v,”’s are to remain finite, 6“ must be real. 
Thus we arrive at the finite allowed band 


| AW) 


< 2K(D) (18) 


as illustrated in Fig. 3. A useful criterion in estimating 
the degree of coupling in an antenna system will be de- 
rived from this band. 

In general, in an antenna system all the coupled 
modes are excited. Let the complex amplitudes of these 
coupled modes be expressed by the column matrix a. 
The overall transverse electric field on the system is 
therefore 


C(t, 2) = a = f'VP(2)a. (19) 


The differences in phase velocities will cause spatial 
beating between modes. This results in a continuous re- 
distribution along the system of transverse field and 
energy distributions around the individual structures. 
It is only with very special excitations, that one coupled 
mode is excited. When this happens no redistribution of 
field and energy along the system will occur. 


PATTERN EVALUATION 


Having determined the coupled surface wave modes, 
we now proceed to evaluate the array pattern using the 
two approximate methods of treating end-fire elements. 
For simplicity, only the pattern in the plane of the array 
will be discussed. We shall treat the two methods sepa- 
rately and then compare them. An experimental method 
for the determination of the interaction coefficient K(D) 
will also be given. 

The first approach assumes the coupled surface wave 
modes excited by the feed system to radiate along the 
antenna structure. The radiation of the rth coupled sur- 
face wave mode as given in (2), in the plane containing 
the array and of a particular polarization, is equivalent 
to that of the following traveling-wave current distribu- 
tion on the system 


yMe-74+aM pz. 


(20) 


The radiation of the current on each structure has its 
phase center located at the mid-point of the element 
(g=L/2). Using the coordinates (Ri, @) as in Fig. 1, 
where R; is measured from the center of the first ele- 


sin 
Fy“(0) = CL 


i & 
a x= [1+ A)6 — k cos 6| me (22) . 


To obtain the total radiation of the N travelling waves 
of (20) we need the space factor determined from the — 
amplitudes of the elements and the array spacing D. Let 
q’ be the following row matrix 


q’ = (1, eikD sin 6 gj2kD sin Qos Goa ei (N-1)kD sin ls 


(23) 
It is readily seen that the radiation of the rth coupled — 
mode is 


qv (eH 1 488 L12) F() (6), (24) 
When all the coupled modes are excited with complex 
amplitude matrix a so that the total field on the an- 
tenna system is as given by (19), the radiation field — 
G.(@) can be determined immediately from (24). If we 
introduce a diagonal coupled mode radiation matrix — 
F,(0) whose elements are Fi(@) of (22), the result be- 
comes . 


Gi(6) = q'VP(L/2)Fi(8)a. (25) 
Assume that the elements can be excited independently — 
with the complex feed amplitudes given by the column 
matrix u. From (19), (15) and (14), there follows 

u = Va. (26) 
Using the orthogonal property of V as given in (13), 
(25) reduces to 


Gi(0) = q’/ VP(L/2)Fi(6) V'u. (27) 

If there is considerable direct radiation from the feed 
system we have to take it into account. Let the first 
feed with unit amplitude contribute a direct radiation 
field F;(6) of the same polarization as in the above dis- | 
cussion but measured with respect to the feed end z=0. 
The direct radiation of the entire feed system when 
measured with respect to the mid-point of the first ele- 
ment is then 


Gy(8) = eHELIP) 009 OF (6)! u. (28) 


-f 


1 
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The over-all radiation pattern from the array is therefore 
P,(6) =~ | Gi(0) “+ Gyi(8) | 
= |q/[VP(L/2)F1(6) V!+ e112) 0089 F-(9)1] |. (29) 
Note that if A™B(L/2)<1, then P(L/2) $e~#4/)1, 
F,(0) + F(@)1 where F(@) is given by (22) with AM =0, 
Thus, when A™6L<1, (29) reduces to 
P,(0) = |q’[F(@)e-8 2/21 + ei (L12) cos6F-(9)1] u| . (30) 


The quantity in the bracket is the element pattern while 
q’u is the space factor determined from the feed ampli- 
tudes and the array configuration. Since A” lies within 
the finite band of (18), if coupling is to be neglected the 


following criterion must be satisfied: 
2K(D)BL <K 1. (31) 


The error in array pattern due to neglecting interaction 


is then approximately of the order of 2K(D)BL. 


Next let us consider the second approach which as- 
sumes the coupled surface wave modes excited by the 
feed system propagate to the radiating end of the an- 
tenna system and then radiate. The transverse electric 
field at the radiating end caused by these modes are 


then, from (19), (26) and (13), 


~~ 


G(r, L) =a =f'VP(L)a=f/VP(L)V'u. (32) 


Let the radiation of a particular polarization caused by 
the aperture field e(t—11) on the first antenna element 


measured with respect to the center of the aperture 
(0, L) be 


e ikke 


F (6). (33) 


2 


The reason for using R2 is that e(r—t) is usually sym- 


~ metrical and almost of uniform phase; hence the phase 


center of its radiation is at the center of the aperture. 


- The total radiation field of the N-component aperture 


field distributions e(r—t,) with the complex amplitudes 


~ given by (32) is therefore, 


G.(0) = F2(6)q'VP(L)V'u, (34) 


_ with q’ as given in (23). 


In addition to the radiation of the surface waves at 
the radiating end of the antenna system, the direct 
feed radiation as given by (28) has to be included. 
However, since it is necessary to measure the radiation 


field with respect to the radiating end of the first ele- 


a 
‘ The over-all radiation pattern from the array computed 


4 P2(6) 


ment, but not its mid-point, the contribution of the en- 
tire feed system in the same polarization is now 


Gale) ce ern OF Og! u. (35) 


from the second method is finally 


| Go(@) + Gy2(6) | 
| q’[F2(0) VP(L)V’ + e-## 09 °F,(6)1]u| . 


lI 


I 


(36) 
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Note that in order to neglect interaction, it is necessary 
that AMBL<1, so that 


P,(6) = | q’[F2(6)e-41 + F,()e-*% 08 Jul, (37) 


where the quantity in the bracket is again the element 
pattern but determined from the second method. Intro- 
ducing the allowed band of A”, the criterion of negligi- 
ble coupling (31) is again applicable. 

Comparing the over-all array pattern expressions 
P,(0) and P,(@) given in (28) and (36), the most pro- 
nounced difference lies in the appearance of matrix 
P(L/2)F,(@) in the first, and of matrix P(L) F2(@) in the 
second. Although matrix F,(9) depends on AL, its 
variation is of a different form and is less rapid com- 
pared with ei 8(L/2) Thus, the effect on the over-all 
pattern caused by element interaction predicted from 
the first method is less severe than that predicted on the 
basis of the second approach. To see this difference in 
more detail a simple example will be used, which will 
also provide an experimental determination of the inter- 
action coefficient K(D). 

Consider a two-element array with the first element 
fed while the second unfed so that 


Note that we do not use equal excitations because only 
one coupled mode will then be excited. Using the 
coupled modes for NV = 2 as given in Fig. 3 the radiation 
pattern can be computed from (28) and (36). Assume 
that near the main beam we can neglect the direct feed 
radiation. The array pattern will have its maximum 
pushed to one side as in Fig. 4. The logarithmic deriva- 
tives of P(6) and P.(@) at 0=0 are 


ae | 
ain 
dé 6=0 


2K(D 
[1+ OE az cote - ® 
B—k 


L\-kD 


= kDsin K(D)B8L. 


6=0 


d 


The two approaches yield the factors kD/z and kD, 
respectively, because the traveling-wave radiation along 
an element depends on the coupled surface waves at the 
mid-points of the element, while the radiation of an 
aperture field depends on the coupled surface wave at 
the radiating end of the element. The last two terms in 
the bracket of (38) account for the difference in patterns 
of the two coupled modes caused by different phase 
velocities. They are usually of the order of K(D)BL so 
that the quantity in the bracket of (38) is approximately 
unity. The predicted value of (38) based on the first 
method is, therefore, less than that of (39) using the 
second approach. Note that from the measured pattern 
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Fig. 4—Pushing of main beam caused by coupling 
in a two-element array. 


_ of such a two-element array the interaction coefficient 
K(D) can be determined, using either method. 

Preliminary experiment on long Yagi arrays seem 
to favour the second method. The values of interaction 
coefficient K(D) obtained from measurement based on 
(39) are closer to the values computed from (8). Fur- 
thermore, the mechanism of radiation assumed by the 
second approach is more reasonable. Thus, it appears 
that in spite of the usefulness of the first approach in 
predicting the pattern of an end-fire antenna, it is less 
appropriate in describing the behaviour of a broadside 
array of such elements. However, this approach should 
be readily generalizable to broadside arrays of tapered 
or modulated elements because it describes the indi- 
vidual radiation of such elements quite adequately. A 
more definite assessment of the two methods, however, 
must await further experimental evidences. 


EFFECT OF SPACING ON PATTERN 


We shall now investigate the nature of patterns ob- 
tainable from a one dimensional uniformly spaced array 
in the plane containing the array. On the basis of the 
second method, the array pattern can be considered as 
the radiation of a virtual aperture field at the radiating 
end of the antenna system whose distribution includes 
both the coupled surface wave modes and the feed sys- 


% D. L, Sengupta, D. D. Tang, and J. L. Yen, “Broadside Arrays 
of Long Yagis,” presented at the Canadian IRE Convention, 
Toronto, Ont., October 7-9, 1959. 
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tem direct radiation. The former dominates near the : 


elements and decreases exponentially away from the ele- 


ments, while the latter becomes important only far from | 


the elements. 


With very wide spacing, the array pattern is simply 
the product of the feed space factor q’u and the element — 
pattern as given by (37) or (30). It is well known that — 
the space factor of a discrete source system with spacing 


D> has more than one main beam in the form of a 
grating interferometer pattern as illustrated by the 
equal ripple pattern of Fig. 5. The angular intervals be- 
tween the interferometer beams are small when the spac- 


ing is large. With sufficiently large D; a few of these _ 


beams may lie within the main beam of the element 
pattern, resulting in a multiple beam over-all pattern as 
in Fig. 5. Such an array is said to be of wide-spaced 
type whose virtual aperture field distribution exhibits 
wide gaps as shown in Fig. 5. Obviously, wide-spaced 
arrays are applicable only for special purposes. 

When the spacing is decreased, the angular intervals 
between the interferometer beams in the feed space fac- 
tor gradually increase. Soon, a spacing Do is reached at 
which the pair of first-order interferometer beams fall 
on the first minima of the element pattern at +0. The 
over-all pattern will therefore have all interferometer 
beams suppressed below the element side lobe leve. For 
an element pattern with sufficiently low side lobes, this 
results in a fan beam (pencil beam for a two-dimensional 
array) with unobjectionable side lobes using a minimum 
element density. An array of this type is, therefore, said 
to be optimum-spaced with optimum spacing Dp» given 
by 


h/ Do =="Sin Bo. (40) 


The main beam shape again depends on the illumination 
or the feed amplitudes in the usual manner. Fig. 6 illus- 
trates the space factor, the element pattern and the 
over-all pattern of the array of Fig. 5, with spacing re- 
duced to the optimum value given by (40), together with 
the virtual aperture field distribution which includes the 
direct feed radiation. Observe that fields of adjacent 
elements now overlap slightly resulting in an aperture 
field with gaps partially filled so that the interferometer 
beams are suppressed. In fact Dy is simply the effective 
element aperture width defined by assuming a uniform 
distribution with the same first minima +6) as the 
actual element. With optimum spacing Do the aperture 
of the array is, therefore, completely illuminated over a 
total width NDp. 

Optimum-spaced arrays are economical fan or pencil 
beam antenna systems because a low element density is 
used. This not only reduces the total cost of the antenna 
elements but also simplifies the feed arrangement. 
Furthermore, if the effective element aperture width D, 
does not vary rapidly with frequency, the optimum- 
spaced condition can be satisfied for a reasonably wide 


’ 
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Fig. 5—A wide-spaced array. 


frequency band. Although the fields of adjacent ele- 
ments overlap slightly as in Fig. 6, the field at one ele- 
ment because of its neighbours must still remain small 
because of the exponential behaviour. Therefore, con- 
dition (31) for negligible interaction is expected to re- 
main valid for optimum-spaced arrays. This type of 
array has found wide applications. 

When the spacing D is reduced below optimum Ds, 
the fields and effective apertures of adjacent elements 
~ overlap so much that the system aperture width is only 
a fraction of ND». Such an array belongs to the close- 
spaced type with a beamwidth greater than that of an 
optimum-spaced array using the same number of ele- 
ments unless a supergain technique is employed. In a 
nonsupergain close-spaced array the greater than opti- 
mum element density can be used to provide a finer con- 
trol of the system aperture field so that no interferom- 
eter beam is evident as indicated in Fig. 7. A practical 
example is the dielectric rod array of Mickey and Chad- 
wick!4 which consists of 17 434 rods spaced at about 
0.6\. In this array the coupling effect is present but 
seems to be neglected in the design of the array. In 


“4 L, W. Mickey and G. G. Chadwick, “Closely spaced high di- 
electric constant polyrod arrays,” 1958 IRE NATIONAL CONVENTION 


RECORD, pt. 1, pp. 213-224. 
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Fig. 6—An optimum-spaced array. 


ARRAY PATTERN 


Fa we Se ae 


VIRTUAL APERTURE FIELD 


Fig. 7—A nonsupergain close-spaced array. 
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contrast to this, the Yagi arrays of Ehrenspecks and 
Kearns’ are also close-spaced, but the coupling is util- 
ized for their operation: One of these Yagi arrays has a 
fed 6 centre element and two 3.2\ parasites spaced at 
0.66. The parasites are excited exclusively by coupled 
surface waves. The amount of coupling is controlled by 
the length and spacing of the parasites. When properly 
adjusted, a clean pattern of very low side lobes is ob- 
tained with beam width not appreciably narrower than 
that of the centre element. Antennas of this type are 
easily analyzed using the present theory. In fact, other 
antennas of this type should be derivable using the 
present analysis. 

If it is desired to reduce the beamwidth of close- 
spaced array below that can be normally obtained with 
the total aperture width, it is necessary to adjust the 
feed amplitudes to produce a narrow-beam space factor 
as so aptly demonstrated by the French investigators” 
using nondirective elements. In an array of end-fire 
elements with narrow element beam the level of the 
space factor outside the element beam need not be kept 
very low. This results in a greater freedom in choosing 
the space factor, and hence, it should be easier to achieve 
supergain effect. It should be possible to design the feed 


% P, Aigrain, “Les antennes super-directives,” L’Onde Elec- 
rique, vol, 33, pp. 51-54; February, 1952. 
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system using the coupled surface waves. The spacing, of © 
course, cannot be so close that the elements can no 
longer be independently excited. 


CONCLUSION 


We have presented two approximate methods of 
evaluating the pattern of a broadside array of end-fire 
antennas on the basis of the coupled surface wave modes 
on the antenna structures. The method assuming the 
radiation to occur at the radiating end ofthe array ap- 
pears to be the more accurate of the two. However, the 
method assuming the radiation to occur along the an- 
tennas may be more readily generalized to tapered or 
modulated elements. Both approaches are approximate, 
applicable only for reasonably wide spacings so that the 
elements can be independently excited. Further ex- 
perimental evidences are necessary to determine the 
accuracy of the predicted patterns. In spite of these 
deficiencies, the use of coupled surface waves should 
lead to a better understanding of the role of interaction 
between traveling-wave end-fire elements in a broadside 
array. 
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A Theory of Antenna Performance in Scatter-Type Reception* 


S. STEIN{, SENIOR MEMBER, IRE, AND D. E. JOHANSENT, MEMBER, IRE 


Summary—Many of the diversity techniques currently proposed 
or in use for scatter reception may be viewed as utilizing antenna 
systems with several adjustable parameters—typically two to four— 
which are continually varied in response to fluctuations in the inci- 
dent waveform, so as to maximize received power. An investigation 
is described which defines an ideal upper bound on performance of a 
linear variable-parameter antenna system in a scatter-type reception 
situation. The important results are 1) a standard for evaluation of 
antenna performance, and 2) deeper physical insight into the prob- 
lem of antenna operation in scatter-type reception. 

A reception model is used in which the receiving antenna views a 
set of radio sources on a hypothetical celestial sphere. For each such 
set of sources possible, a “matched” antenna excitation exists which 
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maximizes received power. Assuming random fluctuations in the 
radio source consistent with those observed in the transhorizon 
phenomenon, statistical properties of the received power are derived 
for the optimized antenna. Similar properties are then investigated 
for standard diversity systems. A significant “gap” is found to exist 
between the performance of existing or contemplated systems and 
the theoretical optimum. 


I. INTRODUCTION 


N recent years, a significant effort has been directed 
if towards developing and improving the performance 

of transhorizon communications at frequencies of 
VHF and above, i.e., so-called tropospheric scatter. In 
order to circumvent the deleterious effects of character- 
istically weak, fading received fields, increasing use has 
been made of such techniques as large radiated power, 
large antennas, and special redundancy techniques usu- 
ally categorized as diversity 


1961 


The present study focuses attention on possibilities 

for further improvements in antenna performance with- 
out additional increases in antenna size. The aim of this 
investigation has been to achieve an increased under- 
“standing of antenna operation under conditions typified 
_by scatter propagation, and of the nature of possible 
optimal antenna systems for such applications.! The 
results presented in this paper provide a useful frame- 
work for evaluating the performance and relative effi- 
ciencies of antenna systems for tropospheric scatter, 
- including, e.g., space and angle diversity [2|-|5]. 


II]. ANTENNA RECEPTION EQUATIONS? 


_ We assume a linearly polarized antenna system char- 
acterized by a planar aperture distribution, an assump- 
- tion consistent with most UHF and microwave large 
_ antenna design. From the well-known Fourier transform 
_ relation between aperture excitation and radiation pat- 
‘tern, and by applying the rigorous electromagnetic 
reciprocity theorem, we can obtain the following basic 
~ statements: 
. In transmitting, let the complex-valued distribution 
_ of electric field in the aperture plane be Eo(x, y), as- 
sumed to vanish outside the antenna aperture S. Assume 
that Eo(x, y) gives rise to a radiated field whose com- 
plex-valued angular pattern is R(a, 6B), where a, 6 are 
_ the direction cosines of a far-field point with respect to 
the coordinate axes x, y of the aperture plane. Then if 
the antenna is a relatively narrow-beam (high-gain) 
- antenna, and is nonsuperdirective, the following trans- 
- form relationships hold (with conveniently chosen pro- 
- portionality constants): 


1 J 127 
—{f Eo(x, y) exp | - me (xa + v8) | dxdy 


Ef os sly +3) Jady (1) 
sab s,9) exp | cee ady 


R(a, 8) 


I 


. 1? idm 

@ E(x, y) = a ip fi R(q, 8) exp E (xa + y@) | dadB, 

‘id —.o 

: where \=wavelength and S is the domain covered by 


the physical antenna aperture. Further, the total radi- 
ated power is 


Py = ff | Ras) |2daas = ff | Eo(x, y) |2dxdy 


= f | Eo(2, y) |2dxdy, 
a s 


d 


(2) 


1 The only related previous work appears to be that of J. Gran- 

lund [1]. Granlund, however, sought an essentially fixed parameter 
system for maximizing mean received power, whereas our concentra- 
tion is on minimizing the major deleterious effects due to the (sta- 

istically rare) deep fades. 

oe adele tis paper in length, it has been necessary to delete 
any amplified discussion of the basic antenna equations. Such a 

detailed discussion, including both the “non-superdirective” assump- 
tions referred to in this section, and a coniplete discussion of the 
reciprocity relationships, is available in [6]. 
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where the bars indicate magnitudes of complex numbers. 

Next, suppose that we are receiving signals from a dis- 
tribution of sources in the far field. This distribution will 
be described by a complex “voltage” distribution f(a, 8) 
defined over a distant “celestial” sphere, with our an- 
tenna located at the center of this sphere, and being 
used for reception. (The scale on f(a, 8) is always inter- 
preted indirectly through antenna measurements and 
need be of no concern here.) The distribution of sources 
corresponds qualitatively to the various multipath rays 
which appear to “originate” within the so-called “com- 
mon volume” of the transhorizon propagation phenom- 
enon. In addition, the sources are all assumed to be 
coherent; 7.e., even though the RF amplitudes and 
phases of these signals are randomized, a basic wave- 
form which is identical for all is associated with each 
source signal. Thus, f(a, 8) may be thought of as de- 
scribing RF carrier amplitude and phase for signals 
emanating from the point (a, 8B). The power detected 
at the receiving antenna terminals, if we use the same 
antenna characterized earlier by its transmitting proper- 
ties, is then readily shown to be 


2 


J fi cc B)f(a, B)dadg 


1 
Pr =— 
Po 


2 


(3) 


[fz y) Fo(x, y)dxdy 


1 
ni 


2 


i} Eo(«, 9) Fo(a, y)dady 
Ss 


1 

an 

where we have additionally defined Fo(x, y), the field on 
the aperture plane due to the distant sources, 


4 | 2a 
Fy(x, y) = — f fit B) exp| — (ne + y8) |daas. (4) 


—3o 


Eq. (3) bears a direct analogy to the circuit theory 
problem of the signal-noise ratio at the output of a 
filter. This suggests the existence of an antenna analog 
to the “matched filter.” It is, indeed, readily shown 
that for a fixed source distribution f(a, 8), the maximum 
available received power is obtained by adjusting the 
antenna feed so that, over the aperture, the feed excita- 
tion factor E(x, y) is a conjugate match to the incident 


field, 


E\(x, y) <3 jG a 23 y)s (5) 


where K is an arbitrary constant; and that this maxi- 


3 We will here treat f(a, 8) as a surface density (in a sense which 
will be apparent later). However, all results are fully general and can 
be applied to a collection of discrete sources by, for example, treating 
f(a, 8B) as a collection of appropriate delta functions with complex 
coefficients. 
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A. Choice of Source Distribution 


Much of the early phenomenological interpretation of ; 
transhorizon propagation assumed that the signal ar-_ 
rives over essentially a continuum of directions of © 


mum available received power is 


(Px) max = if Ea y) |2dxdy, ) (6) 


where S is the antenna aperture.* 
It is relevant, as a side remark here, to note an im- 
portant connection of (6) with the diversity concept. 


arrival, with randomly varying associated strengths. 
More recently, it has been argued [8] that the scattered 
signal may be comprised of a few discrete “rays” (such 


as might arrive by specular reflection from a few moving | 
tilted layers in the atmosphere). Unfortunately, the — 
available experimental data is not sufficient to resolve | 
the question satisfactorily [9]. In view of this uncer- 
tainty, it seems desirable to evaluate the performance | 
for both discrete and continuous source models. In this 
paper, we will limit the discussion to the more general | 


The derivation holds even if.the aperture is comprised 
of a sum of disconnected subapertures, in which case 
the maximum received power can be written 


Pr = x | F(x, y) |2dxdy = 2p Pp,. (7) 


But this is just the output from a maximal ratio diver- 


sity combiner. [7] operating on the outputs of the sev- 
eral optimized subapertures if they are associated with 
equal noise levels. Carrying this argument a little fur- 
ther, it becomes apparent that the optimum antenna 
defined by (6) can be regarded as a diversity system 
having essentially an zmfinite number of degrees of free- 
dom. 

The result in (6) gives the maximum power that can 
be received with a non-superdirective antenna and an 
arbitrary source distribution, and has been derived on 
the basis that the intensity of the source distribution is 
nonfluctuating. In the sections to follow, we will deal 
with source distributions which fluctuate in time (as 
does the apparent source in scatter reception), with the 
details of the fluctuations described statistically. How- 
ever, we now argue that it is valid to assume as our ideal 
antenna system one with variable parameters that are 
appropriately controlled to follow the source fluctua- 
tions with a negligibly small time lag. This is physically 
reasonable since fading rates in UHF transhorizon prop- 
agation are in the order of 1-10 cps, while, for example, 
diversity combiners are in operation with time constants 
of the order of a couple of milliseconds. 

Thus, we can use (6) as a description of the instan- 
taneous received power for our conceptual optimum 
antenna, and the problem is to relate, through this 
equation, the statistics of Pr to those assumed for the 
fluctuations of the sources. Similarly, we will utilize (3) 
to determine the instantaneous Pr for appropriate “non- 
optimum” antennas, and subsequently again the related 
statistics. 


III]. PERFORMANCE OF THE OpTIMUM ANTENNA 


In this section, we apply the results of Section II, 
and evaluate the appropriate statistics for the fluctuat- 
ing power output of the optimum receiving antenna, for 
a specific “extended” source configuration fluctuating 
statistically in time in all its elements. 


“It is well known that only half the incident power flux can 
actually be abstracted from the field. But since our assumed source 
function is to be interpreted only in terms of antenna measurements, 
one should not become unduly disturbed about the “missing” nu- 
merical factor in (6). 


continuous distribution case. Limited computations 
which have also been already carried out for discrete- 
source models have yielded the same qualitative conclu- 
sions as will be discussed here [6]. 


B. Formulation as an Eigenvalue Problem 


We first introduce convenient dimensionless variables 
to denote positions on the aperture plane, £=x/A, 
n=y/d. In terms of these new variables, the incident 
electric field is written F(&, 7) = F(x/A, y/A) = Fo(x, y). — 

The distribution of apparent sources on the celestial 
sphere may be described by an ensemble of complex 
source functions f(a, 8) with a corresponding ensemble 
of images on the aperture plane given, from (4), by 


1 me) 
Fa) =~ f if fle, B) exp [—i2m(Ea + n6)|dadg. (8) 


We assume that the fluctuations in signal carrier 
emitted from any point (a, 8) are uncorrelated with the 
fluctuations in signal carrier emitted from every other — 
point,® and we describe the spatial statistics of f(a, B) 
by the correlation functions of its real and imaginary — 
components. Thus, we first write | 


f(a, B) = fla, 8) + if(a, B). 


In addition, we assume symmetry in the statistics of 
the quadrature components,® with expectations given by 


E| f(a, B1) f(a, 82)] = E[f(ar, 81) flax, Bs) | 


2 


nN 
8a Br) (ar — a2)6(B1 — Be) (9) 


I 


E| f(r, B1)f(a2, B2)] = 0. 


* The assumption of zero correlation from point-to-point is a 
valid mathematical approximation for our purposes to the more 
physically likely situation at the ranges and beamwidths involved in 
practice, that the “correlation distance” on the celestial sphere is 
very small compared to the distances required for significant changes 
in the antenna beam pattern. 

* Some insight into the generality of this assumption may be ob- 
tained from [10]. 
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From (8) and (9), it is readily shown that 


el F@, (2) = ff ie piace C10) 


_ Hence g(a, 8) can be interpreted as the contribution- 


To 


per-unit solid angle on the celestial sphere to the average 


power density in the incident wave front. 


The integral form in (8) for F(é, 7) and the assumed 
narrow autocorrelation for f(a, 8) imply via the central 


limit theorem that the real and imaginary values of 
_ F(&,n), taken at arbitrary points on (€, 7), obey a multi- 
- variate Gaussian distribution. The statistics of F(E, 7) 
_ are thus completely specified by its complex autocor- 
_ relation function. 


E Ro(&1, £23 "1, n2) = E|F (é1, m)F* (és, n2) | 


(11) 


& ide dadBg(a, B) exp { —i2r[a(t: — &) + B(m — 02) ]}. 


_ The right-hand side of (11) is a function only of the two 


APA 


differences h=£—& and t2=71—7. Hence it is con- 


; venient to define 


eee, gy) = ATR Ei bes i, 02), (12) 


_where we have also introduced a normalization with 
_ respect to the wavelength Xd such that 


(1/7) R(Er = £5; 1 72) 


_ is the actual autocorrelation of F. We thus obtain the 
~ form 


— R(t, te) 


= Why 
0? ik J aadsete B) exp {— (at: +1) (13) 


It is noted that R(h, t:)=R*(—h, —t), 1.e., Ris Her- 


-mitian. 


We next turn our attention to the probability dis- 


‘tribution for the functional given by (6), 


pe yt ik | FE, 1) |2dédn, (14) 
S 


which describes the power received by a continuously 
optimized antenna having the aperture S on the (E, n) 


plane. The distribution of such a functional has been 
- determined by Kac and Siegert [11], [12] for the one- 


dimensional case where S is a line segment and F(t) is 
a function of the single variable t. The extension to two 
dimensions is straightforward [6]. The result is that the 


‘Laplace transform of the probability density function 


for Pr is given by 
L(z) = Elexp {—2Pe}] = I] (+ 2/0)", (15) 


_ where the g; are the eigenvalues of the integral equation, 


¢,(&, 7) Ee of RE = &1, Dae ni) or (Es, mi) dé idm. (16) 
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C. Numerical Solution 


The integral equation (16) has been solved analyti- 
cally for a circular aperture and a class of rotationally 
symmetric sources [6], [13]. These rotationally sym- 
metric source functions g(a, 8) are given as the ratio of 
two polynomials, as follows: 


LM = Anat 8?) 


a, 8) = — ——_________, 17 
BONDS se Moan? +B] (17) 


where the polynomials N(x) and M(x) are of degree 
n<m, respectively. Substitution into (13) immediately 
provides the required kernel R(t, t2). Assuming the aper- 
ture to be circular with a physical radius 7, the eigen- 
values o for (16) are obtained [6], [13] as the denumer- 
able infinity of roots of the set of m Xm determinants: 


E®(c) = det {E,;®(c)}; k =0,+1,+2, +--+ +,(18) 


whose (/j)th elements are 


a —|k\+1 x f 
Tega ton (eZ) Kanes) 
r r 
+ yihj lI, @ =) Kyi (x <)I, (19) 


where J;(z) and K;,(z) are the usual modified Bessel 
functions. The constants y;? are the roots of the poly- 
nomial M(x), and the /,? are the roots of the polynomial 


(20) 


By(¢) = 


H,(x) = M(x) — oN (a). 


Treating the h,? as functions of o [via (20) |, the deter- 
minantal equations E“(c) =0 are readily solved for the 
eigenvalues gzn, either graphically or by using standard 
iterative numerical procedures. 

We have computed the eigenvalues ox, for the sim- 
plest physically meaningful source distribution of the 
class considered, 

N(«) = oS 


THs” 


2 


FH Ca ated peal a aren «4 
@=(, ) . 


TK s 


The various coefficients in (21) have been chosen so that 
first [by (10), (17), and (20)". 


F(é, 1) |2] 


QO me dadg = 
< Eine — Ana? + 62) 2 


2 
| Sar,” | 


E| 


(22) 


~ is interpretable as the average power density of the 


electromagnetic wave incident on the aperture plane, 
and second, the intensity of the source at its center is 


[by (17) ] 


Ts 
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It is clear from (22) and (23) that r, may be interpreted 
‘as the “radius” of the equivalent uniform source which 
also has a total flux density Q on the aperture plane.’ 

In the later physical interpretations of our results, it 
will be convenient to introduce the parameter 1./7a 
where rq is the radius of a typical beam “normally” 
radiated by the given aperture. Specifically, 7a is taken 
to correspond to the radius of the equivalent uniform 
beam having the same maximum radiated power den- 
sity, and the same total radiated power as does a uni- 
formly excited aperture of radius r. For the wavelength 
A, it is found that 


Vein HY: (24) 
The eigenvalues computed for the source (21) are shown 
in Fig. 1 (opposite), plotted against 7,/ra.* The eigen- 
value oz, is defined as the nth zero of the kth determi- 
nantal equation E®(c). Actually, many more than 
fifteen eigenvalues are defined by Fig. 1, since for k>1 
the roots zn must be counted twice. 


D. Error Rate Computations 


In the region of low signal levels of primary interest 
in determining the performance of most communications 
systems, the probability density function for Pr can be 
accurately approximated via a saddle point integration 
performed on the integral defining the inverse transform 
of L(z). However, a relation of more direct physical in- 
terest is available. Namely, the probability of error for 
matched filter detection of FSK signals in white noise is 
related to the instantaneous signal power Pr by [14] 


1 Pr 
e(Pr) = —exp| — ; 
2 dy N. 


where y is the information rate (bits per second) and 
N the noise power density in the data channel, referred 
to the same point in the system as Pe. For a fluctuating 
Pr, with probability density function w(Ppe), the aver- 
age error rate is 


(25) 


hye Pr 
E = Ele(Pr)| = mal exp |- aa w(Pr)dPp. (26) 


Since in VHF and UHF communications the significant 
noise source generally lies within the receiver, the noise 
power density in (26) can be treated as a constant. 
Comparing with (15), (26) is recognized immediately as 


7 The mean square width of the source (21) does not exist. Hence 
the radius of the equivalent uniform source is used here as the ap- 
propriate measure of source width. 

8 This choice of parameter coincides with physical intuition which 
demands that except for a trivial renormalization on power, merely 
increasing the radius r of the antenna aperture (and therefore de- 
creasing the achievable beamwidth) should be completely com- 
pensated by proportionally decreasing the source radius /s. 


IRE TRANSACTIONS ON ANTENNAS AND PROPAGATION 


q —1 
1 te Maes : 
E=—I]|][ Qn? (27) 
2 kn Cin os 
re 
where 
Qn? 
i (28) 
2y Nr. 


is a normalized average signal-to-noise ratio, and the 
orn(Qd2/r,”) are just the “dimensionless”-eigenvalues 
plotted in Fig. 1. 

Using (27) and the eigenvalues already determined, 
the error rates have been computed for a range of values 
of ;/ra. These error rate curves are plotted as the solid 
lines (M =1) in Figs. 2-4 for 7,/ra=0.19, 0.63, and 3.14, 
respectively. 


IV. PERFORMANCE OF SOME TYPICAL 
NoNnorptimMuM ANTENNAS 


We next compare the potentialities represented by 
the performance of an optimized antenna with the per- 
formance of typical diversity arrays. We consider both 
a space diversity array and an angle diversity antenna 
receiving from the same source as does the optimum sys- 
tem. We will also generalize the optimum antenna re- 
sults to a spaced array of optimum antennas separated 
sufficiently from one another so that the fading of the 
several signals is statistically independent. The results 
of all three computations are also plotted in Figs. 2-4. 


A. Space Diversity Array 


A common diversity arrangement consists of a spaced 
array of two or more antennas, separated so that fading 
of the signal on each may be considered statistically in- 
dependent of that on any other. 

We assume that each diversity antenna has a Gaus- 
sian beam shape? with equivalent uniform radius 7,. The 
average power received on any particular beam is ob- 
tained by taking the average value of (3) with the aid of 
(9) and with the Gaussian beam shape substituted for 
R(a, 8). Numerical integration of the resulting integrals 
yields. (Pr) = 0.15. .O\?/r,7, 0/019, Ohare ean 
0.0052 Qd?/r,? for the average power received on a non- 
optimized circular aperture when the ratio of source-to- 
beam width is 7;/ra=0.19, 0.63, and 3.14, respectively. 

For the case of Rayleigh fading, knowledge of the re- 
ceived power is sufficient to define the fading distribu- 
tion. The Laplace transform for the fading distribution 
of the received signal is 


L,(z) = [1 + 2E(Pr)]—1. (29) 


® The Gaussian antenna pattern is an approximation to the 
“Bessel function” pattern which results when a circular aperture is 
illuminated by a “flat” distribution of electric field E(t, »)=1. Use 
of the Gaussian shape is justified by its closeness to the Bessel func- 
tion for signals arriving not too far from beam axis, and by the fact 
that lack of sidelobes more accurately approximates the design ideal. 
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Fig. 2—Comparison of optimum antenna with diversity arrays. 
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Fig. 1—Eigenvalues for optimum antenna receiving 
from extended source. 


SOURCE WIDTH fs 
ANTENNA BEAM WIDTH — 


10 
107? 
hor” 
cir No 
= 10- 
5 Ss 
a N 
-5 
uw 10 Ni 
= 
oS 
—-6 
= 10 
w 
: ‘ 
w fo” 
M= NUMBER 
OF DIVERSITY ANTENNAS 
for? ~ 
—— SPACED OPTIMUM ANTENNAS \ 
— —— SPACE DIVERSITY 
lor” 
ro7!® 
-10 10) 10 20 30 40 50 60 70 


q— (DB) 
RECEIVED SIGNAL-TO-NOISE RATIO (NORMALIZED) 


E(q) ERROR RATE BITS/BIT 


a 
S, 


10) 10 20 


eee ANGLE DIVERSITY 


SOURCE WIDTH 


30 40 50 60 


q—> (D8) 


RECEIVED SIGNAL-TO- 


NOISE RATIO (NORMALIZED) 


M=NUMBER OF DIVERSITY ANTENNAS OR BEAMS 
ee SPACED OPTIMUM ANTENNAS 
mom SPACE DIVERSITY 


ANTENNA BEAM WIDTH 1, 


70 


309 


Fig. 3—Comparison of optimum antenna with diversity arrays. 
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Fig. 4—Comparison of optimum antenna with diversity arrays. 
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Assuming independence among the several diversity 
signals, the Laplace transform of the combined output 
(maximal ratio diversity) is just the product 


Lu(2) = Li(z)M = (30) 


[1 + 2E(Pa)|M 


where M is the order of diversity. By analogy to the 
development of (27) the corresponding error rate is then 
obtained as 


in 1 
Eu(q) = — 
OQn2/r,” 


where the parameter g is the same normalized average 
received carrier-to-noise ratio defined previously in (28). 
Curves for M=1, 2, 4, are given in Figs. 2-4. 


(31) 


B. Angle Diversity 


Error rate curves for systems using four-fold angle 
diversity!” and having source to beam-width ratios 
r,/ta=0.63 and 3.14, are plotted in Figs. 3 and 4, re- 
spectively. The derivation of these curves is presented 
elsewhere [6]. In view of the possible lack of generality 
of the model! used there, with respect to real systems 
of-angle diversity, these lengthy derivations are omitted 
here. It suffices to add that in general it is not physically 
possible to place the usual kind of feed horns any closer 
than that which will provide beam overlap at about the 
3-db points. This physical restriction—that the angle 
diversity beam patterns overlap outside their 3-db 
points—implies that the ratio of source width to effec- 
tive antenna beamwidth (2.e., 7,/r2) should not be taken 
very small, since otherwise the angular spread between 
diversity beams will be so large that sufficient coverage 
of the source by all diversity beams will not be possible. 
For this reason curves for the case 7,/ra=90.19 have not 
been included. 


C. Spaced Array of Optimum Antennas 


Assuming that the separate structures of the array 
are sufficiently separated, fading on the individual opti- 
mum antennas within a spaced array will be statistically 
independent. Maximal ratio combining of the outputs 
(equivalent to the entire array being optimized) then 
produces an error rate which is derived by analogy to 
the development of (30), as follows: 


Ey (q) = 3[22.(9)|M, (32) 


10 The particular type of angle diversity discussed is as follows: 
The scatter volume is illuminated by a single transmitting antenna 
so as to provide an effective source as seen by the receiving structure, 
approximated by (17) and (21). The receiving structure consists of a 
single dish, with multiple feeds providing four independent over- 
lapping receiving patterns with a symmetrical projection on the 
celestial sphere. The outputs from each of the four receiving patterns 
are combined in a maximal ratio combiner. The angular spread of the 
resulting pattern is, in every case, selected to provide best perform- 
ance, 7.e., minimum power requirements for a specified error rate. 
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May 
where the functions £:(g) are the error rates computed 
earlier for a single optimum aperture, and M is again 


the number of circular structures in the array. These 
results have also been plotted in Figs. 2-4. 


V. CONCLUSIONS 


Figs. 2-4 present a comparison between the _pre- 
dicted performance of optimum antennas and the per- 
formance of existing or proposed systems of space and 
angle diversity. The proper measure of performance for 
these systems is obtained by comparing power require- 
ments for a specified error rate performance. Several 
conclusions seem warranted: 


1) A single optimum antenna offers a performance 
comparable to that of a standard space-diversity 
array of similar sized apertures. For the smaller of 
the three apertures considered, r,/ra=0.19, the 
average received power requirement for an opti- 
mum single aperture is seen to be the same as that 
for two-fold space diversity, at an error rate of 
10-* bits per bit. A similar statement can be made 
for the r;/7,2=0.63 case. However, in the case of 
the largest of the three antennas, 7./ra=3.14, the 
apparent superiority of the optimum antenna over 
even fourfold space diversity should be regarded 
with care, since a large part of the apparent supe- 
riority is due to the so-called “aperture-to-medium 
coupling loss,” resulting when a fixed parameter 
aperture of large radius is receiving from a rela- 
tively large source. Thus the comparison for 7;/fa 
= 3.14 appears to show not so much the ineffec- 
tiveness of diversity vs the optimum, but rather 
a demonstration that the optimum effectively avoids 
aperture-medium coupling loss. 

2) Further significant improvements should be pos- 
sible in the aperture utilization of existing antenna 
systems: Thus over the useful range of system 
error rates (10-4 to 10-7 bits/bit) the optimum 
antenna outperforms a fourfold angle diversity 
system of the same aperture size by roughly 5-15 
db." Though only a fourfold angle diversity sys- 
tem has been analyzed, it is doubted that con- 
sideration of higher-order angle diversity systems 
would significantly change these conclusions, be- 
cause of difficulties inherent in placing feed horns 
any closer together than that which will cause an 
overlap at the 3-db points. As a further example, 
it is seen that for a ratio of source width to effec- 
tive antenna beamwidth of 7,/r,=0.63, optimiza- 
tion of a fourfold space diversity array will provide 
a power advantage of 6-7 db at the 10~* error rate 
level, over the nonoptimized system. 


We may also note that the performance gap between the 
“standard” diversity arrays and the optimum becomes 
even greater if we discuss error rates in the range 1077 


‘\ These figures are based only on the two larger of the three aper- 
ture sizes considered. 
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to 10-" bits per bit; requirements approaching the 
latter seem to be evident in applications of the near 
future, if not in those already existing. 
A final note may again be made of the fact that all 
computations have been performed for the idealized 
- situation involving a circular aperture and a circularly 
symmetric source distribution having a “radially” sym- 
metric intensity distribution of the form g(x) 
=(cx*+1)—!. Although the absolute performance figures 
predicted in this context obviously will not remain valid 
for less symmetric situations, the relative differences 
existing among system types and the qualitative fea- 
tures of individual systems are expected to be preserved. 
One may conclude that further research is warranted 
for filling the “gap” shown by our comparisons. An in- 
_ dication of why the optimum antenna performs better 
than more standard systems is offered by Figs. 2-4. The 
error rate curves for the diversity systems asymptoti- 
cally approach a straight line, whereas those for the 
optimum bend continually downwards towards the 
vertical. It may be shown that the slopes of the straight 
line asymptotes are proportional to the order of diver- 
sity or number of degrees of freedom, while the con- 
tinual bending of the optimum is related to the earlier 
remark (see the conclusion of Section II) that it is 
equivalent to an infinite number of degrees of freedom. 
Clearly, the direction for filling the gap lies in antennas 
which also equivalently have a large number of degrees 
of freedom. It is hoped to report further research in that 
direction in the near future. 
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CORRECTION 


Klaus Walther, author of “Reflection Factor of 
Gradual-Transition Absorbers for Electromagnetic and 
Acoustic Waves,” which appeared on pages 608-621 of 
the November, 1960, issue of these TRANSACTIONS, has 
called the following to the attention of the Edztor. 

1) Eq. (46) on page 614 should read: 


yo = ol; where: (4)? 


s if to) i ole) 


2) In the sentence following (47) on page 614, the - 
first inequality should read: jq<y.; the second 
inequality should be: 9g 2c. 

3) In the caption of Fig. 18 on page 619, the expres- 
sion g(x) =s(x) —jt(x) =1 should be corrected to: 


g(x) = s(x) —jt(x); fx) =1. 
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communications 


A Scattering Measurement 
Technique* 


King and Scharfman have described a 
technique for measuring the radar cross 
section of objects by modulating the position 
of the object.! The inherent simplicity and 
accuracy of this technique even within the 
confines of a small laboratory very much 
favor its use for many scattering problems, 
although it has been restricted to the meas- 
urement of cross sections of rotationally sym- 
metrical obstacles. 

A relatively simple mechanical impro- 
visation on the experimental apparatus pre- 


ree 8h Se 


Fig. 1—Apparatus used for making scattering meas- 
urements of rotationally unsymmetrical targets. 


* Received by the PGAP, December 12, 1960. 

1 For details of the measurement theory, see D. D. 
King and H. Scharfman, “Antenna-scattering meas- 
urements by modulation of scatterer,” Proc. IRE, 
vol, 42, pp. 854-858; May, 1954. 


viously used with this technique has allowed 
us to measure successfully the scattering 
characteristics of three-dimensional, arbi- 
trarily shaped objects. In response to King’s 
suggestion, this improvisation is reported 
here in the hope that it may significantly 
extend the value of what is already an excel- 
lent experimental arrangement. 

Generally the method consists of meas- 
uring the amplitude of the audio signal that 
corresponds in frequency to the Doppler 
shift caused by a relative motion between 


_an antenna and target. A component of prop- 


er motion necessary for the measurement 
is realized by revolving the scatterer in a 
large circle. Since this motion, however, re- 
sults in a continuous change of target as- 


pect, the technique has been considered to 
be valid only for targets exhibiting rota- 
tional symmetry. The design pictured in 
Fig. 1 overcomes this limitation. Gear 1 is 
fixed with respect to the observer. As the 
horizontal arm is rotated about the drive 
shaft axis, the chain drives gear 2 so that 
the styrofoam column B rotates about its 
own vertical axis with respect to the hori- 
zontal arm as it simultaneously revolves: 
about the drive-shaft axis. The rotation of 
the target effected by the chain drive exactly 
counters the rotation that would otherwise 
occur with respect to an observer during 
the course of target revolution. The as- 
pect of the target thus remains fixed at all 
times with respect to the observer.? Since 


Fig. 2— Oscilloscope display of the return signals from a reference sphere and target. The pulse marks the 
point at which the amplitude of the signal should be read, The first signal is from a 23-inch steel sphere 
and the second from a 24-inch lucite cube. 


2 This neglects the change in aspect with respect to 
an observer due to the lateral displacement associated 
with the circular path of the target. The extent to 
which this may be neglected depends on the distance 
between observer and target. 
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Fig. 3—Backscattering cross section (¢) vs angular position of the cube. Aspect of the cube and 
s 5 styrofoam support with respect to the measured pattern is indicated. 


g 
- column A supports a rotationally symmet- 
‘rical standard scatterer, it may be fastened 
_ directly to the horizontal arm. Oscilloscope 
_ display of the return signal together with a 
_ positioning pulse (see Fig. 2) from a micro- 
- switch located on the rotating member will 
_ provide that both backscattering and off- 
angle scattering can be measured. The pri- 
mary advantage of an oscilloscope display 
is the ease with which spurious noise com- 
_ ponents entering or being generated within 
- the system may be detected. For example, 
- modulation of the amplifier power supply 
_ by the ac motor, which was found to be a 
serious source of error, may have otherwise 
_ been unnoticed. The major disadvantage of 
the oscilloscope is of course the difficulty in 
reading accurately and the time required to 
obtain a number of readings in order to 
minimize the error. Although it has not been 
tried, a peak reading digital voltmeter com- 
manded by the positioning pulse might over- 
come these objections. 
In order to determine the linearity and 
general reliability of the system, the back- 
scattered signals from a set of six steel 
_ spheres ranging in size from one to four 
inches in diameter were measured. A com- 
_ parison of the experimental results with a 
theoretical curve (the only reference readily 
available at present) has indicated that the 
_ experimental error is at least as small as the 
- error in reading the curve. All measurements 
_ were reproducible within +5 per cent, and 
_ by averaging for many readings the error was 
reduced to +one per cent. This indicates 
that the primary source of the error is in 
reading the oscilloscope. Generally it may be 
said that the inherent accuracy of the modu- 
- lation technique is not impaired by the new 
_ arrangement. 
In order to measure the backscattering 
cross sections of rotationally unsymmetrical 
targets as a function of their angular posi- 


A Belk 


tion, a moveable boom was fastened to gear 
1 (see Fig. 1). Manual movement of this 
boom results in a repositioning of the chain 
and, consequently, in a corresponding angu- 
lar change in column and target. An angular 
scale associated with the boom facilitates 
the making of an accurate and reproducible 
change in target aspect without stopping its 
rotation. In this manner the angular posi- 
tion of a 23-inch lucite cube was changed in 
5° increments up to 45° to either side of a 
plane of symmetry and the return signals 
recorded. It was found that the physically 
symmetrical properties of the cube would 
provide a measure of error resulting from in- 
accurate positioning of the target with re- 
spect to the incident electric vector and due 
to the influence of the styrofoam supports. 
This was discovered when the first set of 
data taken with the cube was found to be 
unsymmetrical and when changes in the 
means of target support and compensation 
for the centrifugally-induced target mis- 
alignment resulted in the recording of almost 
perfectly symmetrical patterns. (See Fig. 3). 

Since the support columns are not sym- 
metrical, any influence they have may be ex- 
pected to show up in the form of an un- 
symmetrical measurement pattern, as was 
first obtained. It was found that a reduction 
in size of the styrofoam support in the vicin- 
ity of the target improved the symmetry of 
the pattern. It may therefore be concluded 
that the support columns should be as small 
as possible; just how small can be deter- 
mined by the above method. In other words, 
if the pattern is perfectly symmetrical and 
the support columns are unsymmetrical, 
then the columns can have had no influence. 
Fig. 4 is a composite of Fig. 3(a)(c), and 
represents the errors resulting from position- 
ing of the target [(a)-(c) of Fig. 3 each in- 
volved a new position of the cube], the 
effect of the styrofoam column, and the 


to 


Ks 
Shee 


Fig. 4—Superposition of the points shown in Fig. 
3(a)-(c), indicating the deviation of the measured 
points from an average symmetrical curve. 


systematic error in reading the oscilloscope. 
It should be pointed out that the value of 
o is proportional to the square of the voltage 
ratios read on the oscilloscope. 

These measurements were preliminary 
and are presented to demonstrate the tech- 
nique. Correspondence is invited by anyone 
working with a similar project who may 
either have suggestions or questions. 

K. STEINBACH 

F. B. VARNUM 

U.S. Army Engr. Res. and Dev. Labs. 
Fort Belvoir, Va. 
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Maximum Gain in Monopulse 
Difference Mode* 


It is well known that maximum theoret- 
ical gain of a simple antenna of the large 
aperture type is obtained if the aperture is 
uniformly illuminated with no spillover. 
For a monopulse antenna, it has been shown 
that maximum theoretical slope gain in the 
difference mode is obtained with odd linear 
illumination having no spillover.? 

It is generally accepted that slope gain is 
the most significant measure of performance 
in the difference mode, because the slope at 
the center of the difference pattern is a factor 
in the ability of a monopulse antenna to ac- 
curately track a target. However, the gain 
at the peaks of the difference pattern is also 
significant, since it may affect the ability of 
a monopulse antenna to achieve angle lock- 
on of a distant target. Furthermore, peak 
gain is a term which is universally employed 
in the antenna field, and is a standard meas- 
urement in the testing of an antenna. In this 
note the aperture illumination which pro- 
vides maximum theoretical peak gain in the 
difference mode will be determined. 

Consider first the field which a distant 
source would create on the aperture of an 
antenna during reception. If the source is 
located off the antenna axis, the distribution 
of this field may be separated into even and 
odd components. Since only the odd com- 
ponent couples to the difference mode of the 
antenna, only this component is of interest 
here. The shape of this odd field distribution 
is a simple sine curve, as indicated in Fig. 1. 


Fig. 1—Odd component of field received from a 
distant off-axis source. 


It is described by the following relationships: 
F(x) = ¢1 sin kx (1) 
<3 Qa re 2a sin 05 


Baar x (2) 


where F,(x) is the amplitude of the received 
field distribution, c, is a constant, x is the 
aperture coordinate, Aq is the wavelength of 
the sine distribution across the aperture, 0, 
is the angle of the source, and 2 is the free- 
space wavelength. 

What shape should the aperture illu- 
mination of this antenna (during transmis- 
sion) have in order to obtain a maximum 


* Received by the PGAP, December 19, 1960. 
oe S. Silver, “Microwave Antenna Theory and De- 
sign,” M.I.T. Rad. Lab. Ser., McGraw-Hill Book Co. 
Inc., New York, N. Y., vol. 12, pp. 177-178; 1949. 
2G, _M. Kirkpatrick, General Electric Co., 
ASTIA Doc. No. AD 18458, pp. 13-15; August 1952. 
Also see E. J. Powers, ASTIA Doc. No, AD 
231217, pp. 8-15, 66-72; July, 1959, 
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signal (during reception) from the above 
field distribution? By analogy with the 
“matched filter criterion”® the illumination 
should have exactly the same shape as the 
field distribution. Thus an aperture illu- 
mination having a sine shape should yield 
maximum gain the difference mode. The 
wavelength of the sine illumination should 
equal the wavelength of the field distribu- 
tion, which, in turn, is directly related to the 
angle of the source. The following relation- 
ship describes this matching aperture illu- 
mination: 


sin kx 
L/2 1/2 
( sin? ) 


-L/2 


F,(x) = ce 


(3) 


where F,(x) is the amplitude of the aperture 
illumination during transmission, C2 is a con- 
stant, and L is the extent of a one-dimen- 
sional aperture. The denominator term con- 
serves power across the aperture. 

When the antenna aperture is finite in 
extent, there is one source angle which per- 
mits the greatest signal to be obtained. This 
angle is the angle of the peak of the differ- 
ence pattern. Correspondingly, there is one 
wavelength of the sine illumination which 
yields the maximum peak gain with a given 
finite aperture. This may be determined in 
the following manner. 

The voltage received by an antenna is 
given by 

Lie 


E=¢ F(x) F,(x)dx (4) 
/2 


where £ is the received voltage and ¢; is a 
constant. This is easily demonstrable by 
means of superposition and reciprocity, or 
by analogy to filter theory. Now combining 
(1), (3), and (4), the received voltage is 

sin? kx 


La 
k= af ——_———_————dx. (5) 
L/2 12 
cae (f sin? kad ) 
-L/2 


Evaluating this integral gives the following 
simple relation for received voltage as a 
function of k and L: 


C4 sin RL \}/2 
E = —({L -—- —— 4 
2 ; (6) 


In order to obtain the relation between 
k and L which gives the maximum received 
voltage, (6) is differentiated with respect to 
k and the derivative is set equal to zero. This 
gives 


tan kL = kL. (7) 
Solving this transcendental equation gives 
woe 0.715 
2 Sw Tv 
L 
== = 128.7°, (8) 


Thus the wavelength of the sine illumina- 
tion for maximum peak gain in the differ- 
ence mode has been determined. 


5 J. H. Van Vleck and D, Middleton, “A theoreti- 
cal comparison of visual, aural, and meter reception of 
pulsed signals in the presence of noise,” J. Appl. 
Phys., vol. 17, pp. 943-944, 964; November, 1946. 


The resulting illumination is shown, 
normalized to its peak, in Fig. 2. This illu- 
mination is a truncated sine curve having an 


edge taper of 2.15 db, and happens to — 


closely resemble that which exists in many 


monopulse antennas of the type in which a_ 


parabolic reflector is illuminated by a mono- 


pulse feed. However, in those cases the gain — ) 
in the difference mode is far below its maxi- — 


mum possible value, because of the large 
amount of spillover. Naturally, there is no 
spillover for this theoretical case. 


The antenna pattern may be calculated | 


for this illumination by means of the usual 
techniques, or may be inferred from the 
truncated-sine illumination. The pattern, 


normalized to its peak, is shown in Fig. 3. _ 


iS 


Fig. 2—Aperture illumination for maximum 
gain at the difference peak. 


-1.0 


Fig. 3—Difference pattern having 
maximum peak gain, 


Its value, normalized to the unit-amplitude 
illumination of Fig. 2, is given by the follow- 
ing formula: 


1 sin (u — Um) 1 sin (2 + 2m) 


E(u) = = 
2, Uu — Um, 2 uw + Um 
where 
aL sin 6 3 
“= —— 
: (9) 
and 
wL sin Om 
Un = — = 0.7154 
ON 
4J. F. Ramsay, “Fourier transforms in aerial 


theory—Part II, Fourier sine transforms,” Marconi 
Rev., vol. 10, p. 17; January—March, 1947. 


: tay | 
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where E(u) is the pattern amplitude, @ is the 
pattern angle, and 6,, is the angle of the pat- 
tern peak. This pattern is simply the differ- 
ence of two ordinary maximum-gain pat- 
terns, squinted to the difference peak angles. 


‘It is interesting that the squint angle is such 


that the first sidelobe peak of one of these 
“lobe” patterns coincides with the main 
peak of the other; this condition would be 
expected to correspond to maximum peak 
gain in the difference pattern if two such 
lobe patterns could be independently 
created with a single antenna structure. 

The value for difference peak gain may 
be expressed by comparing it with ordinary 
peak gain. This ratio is determined as fol- 
lows: 


En, 
coe (ttm) (10) 
Zo 1 ue. Ww 
= sin? (0.715 =) dx 
L J_ie Ie 


where g is the maximum difference peak 


gain, and go is the ordinary maximum peak 


gain obtained with uniform illumination of 
the same aperture, both cases, of course, as- 
suming no spillover. The denominator repre- 
sents the power of the unit-amplitude illu- 
mination of Fig. 2. Evaluating this expression 
yields 2 


eo int 01155 — 0,009. 
£0 


(11) 


Thus the maximum possible gain in the dif- 
ference mode of a monopulse antenna is 2.15 
db below the maximum possible gain in the 
sum mode of the same antenna. It is inter- 
esting that this number exactly equals the 
2.15-db edge taper mentioned in connection 
with Fig. 2. It may also be noted that this 
gain ratio is reasonably close to the 3-db 
ratio which would be expected if the two 
previously-mentioned “lobe” components of 
the difference pattern were separated by a 
squint angle great enough so that only one 
component contributed significantly at the 
difference peak. 

In conclusion, it is appropriate to men- 
tion that a knowledge of the illumination 
which provides maximum theoretical peak 


- gain in the difference mode of a monopulse 


antenna can be helpful in the design of such 
an antenna, even if one chooses to avoid this 
illumination because of the high sidelobes 
it yields. Furthermore, the value of the 
theoretical peak gain can provide a useful 
standard for rating the efficiency of a mono- 
pulse antenna. Of course, the results pre- 
sented here apply exactly only-to a onesdi- 
mensional aperture or to a rectangular aper- 
ture with separable distributions. 

It is perhaps also worth noting that the 


’ “matching” method employed to determine 


the optimum illumination involved a simple 
physical process rather than a complex 


4 mathematical one. The method is applicable 


to other similar problems; for example, the 

illumination for ordinary maximum peak 

gain and for maximum difference slope gain, 

mentioned at the beginning of this com- 

munication, may be solved by inspection 
with this method, 

PETER W. HANNAN 

Wheeler Labs. 

Smithtown, N. Y. 
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Current Distributions on Cylinders 
Excited by Spherical 
Electromagnetic Waves* 


INTRODUCTION 


Relationships are established between 
the current distributions on infinitely-long 
perfectly-conducting cylindrical scatterers 
of arbitrary cross section, illuminated by 
1) plane waves and 2) spherical waves. 
These are of value because, although most 
analytical treatments of scattering by cylin- 
ders assume plane-wave excitation, most 
physical situations involve spherical- or 
cylindrical-wave illumination. 

It is known that the cross-sectional dis- 
tributions of current on circular and elliptic 
cylinders! excited by a spherical wave with 
a far-zone phase center are only locally de- 
pendent on the incident field. An important 
example of this is the following: Suppose 
that the incident spherical wave at some 
cross section of a circular or elliptic cylinder 
is identical locally to a plane wave propagat- 
ing at some angle to the cylinder axis; then, 
the transverse-current distributions at that 
cross section for both illuminations are also 
identical. This result is valid for cylinders of 
more general cross section, subject to certain 
restrictions on cylinder size and source dis- 
tance, and is proven for waves due to electric 
sources polarized in a plane perpendicular to 
any radial line [see p: of Fig. 1(b)] and, spe- 
cifically, for sources polarized parallel and 
perpendicular to the cylinder axis in that 
plane. Vector addition recovers the more 
general case. 

An additional relationship between the 
transverse-current distributions at different 
parallel cross sections of a cylinder when 
excited by a spherical wave, or at any cross 
section when excited by plane waves inci- 
dent at different angles—a relationship that 
is known to hold for circular and elliptic 
cylinders!*—is shown to extend to arbitrary 
cylinders: The components of the transverse 
current distributions for 9>=7/2 are identi- 
cal to those same components at any other 
value of 6) on a larger cylinder—larger by a 
cross-sectional scale factor (sin 00)~!. [See 
Fig. 1(a) and (b).] Therefore, the func- 
tional dependence of these components de- 
pends upon the product of sin 4o and a char- 
acteristic dimension. 

With these relations, much of the data 
on current distribution induced on a cylin- 
der by a spherical or an oblique plane wave 
can be obtained from the normal plane- 
wave case. Experimental simplification is 
also suggested, since current distributions 
for normal plane-wave incidence on a range 
of cylinders of the same cross section but of 
various sizes can be obtained through the 
use of one cylinder of the largest size de- 
sired, and a source of spherical waves. The 
current distributions on all smaller cylinders 


* Received by the PGAP, January 6, 1961. This 
research was made possible by support extended by 
Cruft Lab., Harvard University, Cambridge, Mass., 
under contract No. AF19(604)-786. 

1W. S. Lucke, “Electric dipoles in the presence of 
elliptic and circular cylinders,” J. Appl. Phys., vol. 22, 
pp. 14-19; January, 1951. 

2G, Sinclair, “The patterns of antennas located 
near cylinders of elliptical cross section,” Proc. IRE, 
vol. 39, pp. 660-668; June, 1951. ; 

3 Pp, S. Carter, “Antenna arrays around cylinders,” 
Proc. IRE, vol. 31, pp. 671-693; December, 1943. 
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may then be obtained by varying, axially, 
the cross section at which transverse-current 
measurements are made. 

Proof of these properties utilizes integral 
equations for the electromagnetic field* re- 
cast into an appropriate form by a Fourier 
integral operation, and integrated by the 
method of steepest descents to obtain the 
first terms of series expansions in powers of 
(kRo). Retaining only the first terms pro- 
vides results valid in the far zone of the 
spherical-wave source, kRo>1. 


CROSS-SECTION COORDINATES 


E! (CASE A) 
H! (CASE B) 
POINT OF 
OBSERVATION 


| 
DIRECTION OF 
PROPAGATION 


| 
ner ee 


(Ri, 8), $0) 
PHASE CENTER 
OF SPHERICAL WAVE 


DIRECTION OF 
INCIDENCE 


(b) 


Fig. 1—(a) Configuration for plane-wave scattering; 
(b) configuration for spherical-wave scattering. 


4jJ. A. Stratton, “Electromagnetic Theory,” 
McGraw-Hill Book Company, Inc., New York, N, Y. 
p. 466; 1941. 
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FORMULATION OF THE INTEGRAL EQUATIONS 


The integral equations are: 


Case A 


H’, the incident magnetic field, per- 
pendicular to cylinder axis 2: 


BAS) = Gan i di! 


$.(1 +a <) G(S, S)K.(S')de't ; (la) 
Case B 


E’, the incident electric field, perpendicu- 
lar to cylinder axis, 2: 


H+(S) = K.(S) — § az} if Ku(S’) sin (y’ — ) 


c 


+K,(S’) E Ge) at 
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Case A 


A (p, g, 25 4, 
sin Oo 


Case B 
—A(p, , 2; 80, bo) sin (¢o — y) Cos Bo 


= K,(p,¢, 2) — Gp ae’ jit cos OK v(0', ¢’, z) sin (7y/ — y) Ho (| o- 0’ | Rk sin 80) 


fe) 
+ Kilo, 2) | cos (y — 4) 
Op 


A(p, 9g; 23 4, ¢) sin > tne Ki(o, 9g, z) 
in (g’—y) 9 


" fe) 
~ 2 § Kul, 6,9) [== — — cos (¢’ — 7’) | Ho(|@ — 9’ |ksin@)de’. (3c) 
c! p a 


p ag’ 


G(S’, S) 


02’ 


en 


GES. , 
F = (s',Syac't 


(1b) 


HG) == KS) — $ az if Ks | = ae ee ae, = Jew’ sac » (1c) 


, 
- p 


where coordinates are as defined in Fig. 
1(a), and where either E* or H* is perpendic- 
ular to the cylinder axis; where c is cylinder 
cross-section surface, S is shorthand for po- 
sition coordinate on c, K(S)=%XH(S) 
=surface-current density, 7% is the outward 
normal from c, and 

G(R, R’') = exp ik|R — R'|/4cR — R’|. 


PLANE-WAVE INCIDENCE [FIc. 1(a)] 
For Case A: 


EZ = A(p, $, 2; 90, $0) sin Oo. (2a) 
For Case B: 
Hy = — A(p, ¢, 2; Oo, bo) sin (69 — y) Cos > (2b) 
Hz = A(p, ¢, 2; 90, 60) Sin Oo. (2c) 


A=exp —ik[z cos 00+p sin 09 cos(¢ —¢o)]; 
40, #0 indicate the incident wave’s propagat- 
ing direction. 

Utilizing the Fourier integral form of 
G(R, R’)’—the fact that K(.S) can be repre- 
sented as K(s) exp —zkz cos @>—and the 
identity 


2n5(4) = f exp tx2'dz! 


in the right sides of (1a) through (1c), and 
(2a) through (2c) on the left sides, we have 
the following result (s is the two-dimensional 
position coordinate on c): 


5 These may be obtained by applying the boundary 
conditions for a perfectly conducting surface to the 
equations given by Stratton,’ then specializing them 
to the appropriate coordinate system and incident- 
field polarization. Rationalized mks units and time 
dependence e~*”? are assumed. 


6 Zt C3 erve§ 
G(R, R) =f 10(\g - | VE =a 
Ae 


exp ta(z — z’)da, 


In each case, p and p’ are on ¢ and c’. Note 
that the cross-sectional dependence of Kz, 
case A, and K;, case B, are identical to 
those on a cylinder of identical shape illumi- 
nated by a normally-incident plane wave 
but with transverse dimension scaled-down 
by sin 00. Note also that for normal inci- 
dence K,, case B vanishes. 


SPHERICAL INCIDENT WavEs [Fie. 1(8)] 


The sources are polarized in the z direc- 
tion, case A, or perpendicular to the pi—z 
plane, case B, at Ri, #1, ¢@o. The incident 
fields are representable as follows: 


Case A 


Ef =B f Eye, a, 6) 


Ho™+(| 9 — 01 | Vk? — ae'(2 da; (4a) 
Case B 
Hi = weuC sin (bs — 7) fault, « 42) 


Ho™+(| o— 01 | Jk? — a®)eia(z21)q (4b) 


. } . = d 
H,;= iweuC f Fy(k, a, 62) -— 
Ea Op. 


H+ (| Ace o1| Je — a®)eia(a)da: (4c) 


where the F’s are field factors, assumed to 
be well-behaved in a and ¢2 with possible 
singularities in the a-plane only at points 
a= +k, the branch points of Hy, 

Substituting (4a) through (4c) for the 
left sides of (1a) through (1c), and utilizing 
the Fourier integral of K, 


K(p1, po, 2, 21, 5’) 


= if K(p1, do, 2, 21, Sy B)e®'dB 


$e) _ ras Ho(| 9 — 9'| # sin 00)Kz(o', #, 2)dc’. 


sin (9! 


(the p1, ¢0, $2 and % factors are introduced to 


7 


aD Soe melbsinanf | 


| 


exhibit the functional dependence on source 
position)—and the aforementioned 6-func- 
tion identity in the right sides, equations are 
obtained in which both sides are of the form 
of Fourier integrals. We shall refer to them 
as (5a)-(5c) despite the fact that space limi- 
tations prevent their explicit presentation. 
Since the vector components of the inci- 
dent fields are, at worst, sectionally con- 
tinuous on the surface, the integrands (or 
transforms) are identical (in the sense of — 
Fourier convergence). 

In the limit 2| 9:—e@| >>1, pi>>p, with p on 
c, and the cylinder in the far zone of the 
source, (4a) through (4c) [=the left sides of 
(5a) through (5c)] are integrable by the 
method of steepest descents.’ First-order 
integrations valid in the above limit yield 
terms of order (kRo)~1.8 These represent the 
first terms in series of inverse powers of kRo 
(in which higher terms are negligible), and 
hence are equal to like terms on the other 
sides [(5a) through (5c)]. 

Equality of integrands (in the Fourier 
sense) on each side of (5a) through (5c) 
specifies also equality of form and hence an 
identical far-zone relationship on the right 
of (5a) through (5c) between each integral 
and its integrand.’ Further examination of 
(5) shows that # may be decomposed in the 
limit pi>p, R| e:—| >>1, since 


K (01, $0, $2; 21; 5, a) 
= A(p1, do, $2, 5, a) exp i[orVk? Sr Sie a8) | 


where h is a slowly-varying function of a. 
Utilizing the above relationship and the 

Fourier integral, the first term in the series 

of powers of (kRo)~ for K is, in the limit, 


K(p1, po, 2, 21, S) 


2a 


K(p1, bo, p2, 2, 5, & COS Oo)R sin O 
kRo 


of & 
“exp — i(= = ks c05 6). 


7 The integrands are of the form oot Sn 
D(p1, bo, $2, p, >, k, x) exp ila(z — z1) + pi VR? — a], 
where D is a slowly-varying function of a compared to 
the exponential. The saddle points are at a =k cos >. 


_ § Each term is k ¥2a/kRo sin 0) exp —ix/4 times 
its yoke grand: evaluated at the saddle point, 
@= COS Go. 


we 
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Inserting the indicated saddle-point 
evaluations on both sides of (5) and then 
substituting the above representation for K, 
the following are obtained, in the limit, for 
terms of order (kRo)!: (It is important to 
note that the (kRo)~! dependences on the 
right-hand sides of these equations are car- 
ried implicitly by the K’s.) 
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INTRODUCTION 


In the ever-growing study of plasmas, 
perhaps the most-sought property is the 
electron density. The interaction between 
an electromagnetic wave and the plasma 
has provided a useful tool for measuring this 
density. It is well-known! that in the ab- 


exp ik[Ro — p sin 0) cos (6 — ¢o)] = {right side of (3a) } c (6a) 


B2ikt 
aR C sin (¢o—+) cos O.F's(k, k cos 0, $2) exp ik[Ro—p sin 0) cos (6 — $0) | = {right side of (3b)} (6b) 


Case A 
2kBFa(k, k cos 0, $2) 
RRo sin? Oo 
Case B 
wkRo 
2ik4 
wkRo 


each evaluated in the limit for p on ¢, p’ on c’ 
‘and with K(pi, ¢0, ¢2, 2, 5s, 2) replacing 
K(p’, ¢, 2). 
These are the desired forms of the (kRo) 
terms of the integral equation for spherical- 
wave illumination. If the ratio of the largest 
dimension of the cylinder cross section to the 
source distance is large, and if the far-zone 
- relationship is valid (the limiting conditions 
assumed above), F, and F; are constants in 
¢2 on the cylinder surface. Then, for equal 
values of 0) (polar angle of incidence), the 
cross-sectional dependences of the left sides 
of (6) and (3), respectively, are equal, and 
hence the respective cross-sectional distribu- 
tions of current are identical to within a 
multiplicative constant (dependent upon 
source strength, distance, and axial posi- 
_tion). As before, the cross-sectional de- 
pendences of the components of current 
existing on a cylinder scaled-down under 
normal illumination by sin 0 are equal to 
those on the unreduced cylinder at the cross 
section where the angle of incidence is 4. 
D. B. Brick 
Applied Research Lab. 
Sylvania Electric Products, Inc. 
Waltham, Mass. 


* The Potential Utility of Scanning 
Microwave Beams in Plasma 
- Diagnostics* 
SUMMARY 


A technique is proposed for improving 
the data rate in certain plasma experiments 
wherein the electron density is measured. A 
scanning microwave beam is used in con- 
junction with a de magnetic field to create a 
_ plasma resonance dependent on scan angle. 
Electron density in the plasma can then be 
- measured by noting the beam angle where 

resonance occurs. 


* Received by the PGAP, January 9, 1961. This 
study was supported by Lockheed Missiles and Space 
Division. . ; 


CF,(R, & cos Oo, 2) sin A exp tk [Ro — psin 6 cos (6 — ¢o)| = {right side of (3c)} san (OG) 


sence of a dc magnetic field, the plasma can 
exhibit an equivalent dielectric constant 


given by 
Wp” 
emai 


in which eo is the free-space permittivity, 
w is the angular frequency of the electro- 
magnetic wave, and 


O50 (2) 


with 7 the electron density, e the magnitude 
of the charge on the electron, and m its mass. 
wp is called the plasma frequency. A plane 
wave propagating through the plasma will 
experience a resonance, when wp=w, which 
is clearly discernible in either the transmit- 
ted wave or the reflected wave. Hence, for 
a plasma with fixed properties, varying w 
till resonance is observed gives a measure of 
w,and thus x. Similarly, for a plasma of vary- 
ing properties, use of a judiciously chosen 
fixed w will give knowledge of w, and ” at the 
instant resonance is observed. Both experi- 
ments form a part of plasma diagnostics, 
but the second experiment is far more com- 
mon. 

A major difficulty in the technique just 
described is that, unless the boundaries are 
well defined, only one piece of data at reso- 
nance is obtained for each performance of the 
experiment. In many situations, it is desired 
to know w, as a function of time. For a re- 
peatable experiment, w,(t) can be synthe- 
sized by performing the measurement over 
and over again for a sequence of w values. 
However, this can be costly and introduces 
additional sources of error. 

By using a dc magnetic field perpendicu- 
lar to both the electric and Poynting vectors 
of the electromagnetic wave, a triple reso- 
nance phenomenon can be created? offering 
the possibility of an increased data rate. 
However, the increase is modest, and clearly 
what is preferred is a technique for measur- 
ing the entire curve w,(#) in one experiment. 


1 See, e.g., S. K. Mitra, “The Upper Atmosphere,” 
The Asiatic Society, Calcutta, India, 2nd ed., ch. 6; 
1952. 

2 Ibid., Fig. 9. 
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Accordingly, a search has been under- 
taken for such a technique, with some small 
measure of success, as will be shown in what 
follows. 


ANALYSIS 


In a plasma region pervaded by a uni- 
form dc magnetic field, which for conven- 
ience is taken to be z-oriented, it can be 
shown? that the equivalent permittivity is a 
tensor, given by expression 

er «2 O 
e=e|—ex en 0 (3) 
0 0 €33 


in which e is the permittivity of free space 
and 


@p” vy + jw 


an = ome Gt ja)? best rary ere (4) 
ee rg ee es ae 
ie ae GRO one ee 
ane 1 
oy Ce Pea ees 
SSR rag. ome a “ 


In the foregoing, v is the collision frequency, 
and 


eBy 
m 


(7) 


op = 


is the cyclotron frequency, with Bo the 
strength of the dc magnetic field. 

Assume a plane wave is traversing the 
plasma region with its Poynting vector mak- 
ing an angle @ with respect to the dc mag- 
netic field. No loss in generality results from 
placing the Poynting vector in the x-z plane. 
Then each component of the plane wave is 
expressible as a constant times the expo- 
nential factor 


ei (@t-kgX—heZ) (8) 
with 

kr = k sin 0 (9) 

kz = k cos0 (10) 


in which k is the complex propagation con- 
stant of the plane wave. 
Using (8), the tensor wave equation 


VXV X E= wcE (11) 
yields the three component relations 


(Ro%e11 — Re?) Ex + hoei2Ey + kekzEz = 0 


Rorer2 Le + (Rk? = Ro*e11) Ey = 0 (12) 
kakzEz + (Ro2€33 aa hy”) Ez =5()), 


To obtain a nontrivial solution of (12), one 
sets the determinant formed from the co- 
efficients of the components of E equal to 
zero. When this is done and use is made of 
(9) and (10), we obtain the dimensionless ex- 
pression 


BN\4 
[en sin? 6 + €33 cos? 6| (=) 
0 


Rk 2 
—~ [ (ex1”) +12?) sin? 6+ e11€33(1 + cos? @) ] (~) 
0 
+ (en? + e12”)es3 = 0. (13) 


3 W. P. Allis, “Motion of ions and electrons,” in 
“Handbuch der Physik,” vol. 21, pp. 383-444; 1956. 
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It is instructive to examine (13) for the 
special case 9=7/2. The solutions are then 
easily shown to be 


( ) 
al 
en? + ee? + enress + [eu? + €12? — e11€29| 


= - (14) 


2e11 


We shall refer to the minus sign solution as 
case 1 and the plus sign solution as case 2. 
inserting (14) back into (12) and noting 
that k,=0 for both cases (since @=7/2) we 
find, neglecting loss: 


Case 1 
k\? Wp? 
eee tl ea 
( =) a w 


E, = E,=0 E, #0. 


(15) 


(16) 


This is the conventional case referred to 
in the introduction and is independent of the 
presence of the magnetic field. One reso- 
nance occurs when wp=w. 

Case 2 


( k y = e117 + €12” 
Ro €12 


E,=—-“n,#0 E,=0. 


€11 


(18) 


This is the second case referred to in the 
introduction and provides three resonances, 
when 


wp? = w(w — wr), (19) 
wp? = (w + wr) (w — wr), (20) 
Wp? = w(w + ws). (21) 


In principle, these three resonances can 
be spread apart at will merely by choosing 
wp to be a big enough fraction of w. In prac- 
tice, it is difficult to get large fractions if 
» is high up in the microwave range because 
of the large values of By required. However, 
even as high as X band, w can be made large 
enough to separate the three resonances 
sufficiently to argue that a triple data rate 
is possible. 

We therefore draw the important con- 
clusion that the two branch solutions of (13) 
have, respectively, one and three resonances 
_as 0-7/2. It seems reasonable to expect 
that these resonances will still be evident 
at angles of @ other than 7/2. Since the co- 
efficients in (13) are functions of 0, wy at res- 
onance should be a function of 6 as well as a 
function of # and w. This opens up the possi- 
bility that a scanning microwave beam (6 
variable) will yield a resonant relation wp, 
(@) for fixed values of w and wy. Thus, in ex- 
periments in which w, is varying, recording 
the variable angle @ at which resonance is 
observed might give data which could be 
translated into w,(t). 

The branch of case 1 is initially more 
attractive, since a triple resonance is not 
needed if sufficient dynamic range is avail- 
able in the resonance relation.w,(@) for the 
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single resonance case. It is convenient to 
make the substitutions 


a= en, (22) 
b? = 11633, (23) 
co? = en? + e102”, (24) 
r= (4), (25) 
ko 
which convert (13) to 
[a? sin? @ + 6? cos? 6c? 
— [ac? sin? @ + ab?(1 + cos? 6) |¢ 
+be2= A+ BE+C=0, (26) 


the general solution of which is, of course, 


B fe ( B y G 
=—=—+ —)— —- (27 
eee yew ary le ee 
From (15), we note that the resonance in 
case 1 is a zero. Likewise, from (17), it is 
clear that the resonances defined by (19) 
and (21) are zeros. The motions of all three 
of these zeros with @ can be deduced by set- 
ting ¢=0 in (27). But this results in the con- 
dition 
€33(€117 + €12”) 


€11 sin? 6 + €33 cos? 6 


Cc 
(= 28 
0 (28) 


We assume that By +0 and that w, is finite. 
Then e163, and in the range 0<@<7/2, 


the denominator is never zero nor infinite. ~ 


Either 
(29) 


€33 = 0 
or 
11” + €12” = 0. (30) 


But (29) is equivalent to wp,=w, the zero 
resonance of case 1, and (30) is equivalent 
to (19) and (21), the zero resonances of case 
2. Thus none of the three zero resonances 
move with @. case 1 is ruled out as a possi- 
bility, and in case 2, the only hope is the 
pole described by (20). 
To study the behavior of the pole, we set 
1/£=0, which leads to the condition 
A 0 €11 sin? @ + €33 cos? 0 
é €33(€112 + 12?) 


which has as its principal result 


(31) 


€11 Sin? 9 ++ €33 cos? = 0 O oes (32) 


Using (4) and (6) and neglecting loss, (32) 
can be put in the form 


_ oe 
@). w 
oS 4 2 : (33) 
(23) 
1 — —cos?@ 


and the plasma frequency at pole resonance 
is seen to depend on the angle @. A plot of 
(33) is shown in Fig. 1. 

Inspection of this figure reveals that if 
the cyclotron frequency w» is a small fraction 
of the plane wave frequency w, only a small 
variation of plasma density required for 
resonance with scan angle occurs. It takes a 
precisely controlled small excess of w over 
wp, to give a big dynamic range in w» at reso- 
nance vs @. It is highly improbable that this 
precise a control can ever be achieved. 
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Fig. 1—Plasma density at resonance vs scan angle 
with dc magnetic field as parameter. 


The situation is much more interesting 
if ws exceeds w. It is evident that several 
orders of magnitude of wp can be covered in 
a reasonable range of scan angles if w» is of 
the order of twice w, with the exact ratio not 
too critical. 

This situation is not without its difficul- 
ties. By should not be too large, as it might 
seriously affect the decay process in the 
plasma. This places an upper limit on w, 
which raises the problem of achieving a 
suitable beamwidth with the available aper- 
ture space. Thus, compromises must be 
made, but it appears that the technique 
should be practical in certain experiments 
and it seems well worth the try. The present 
state of the art insures reading the scan 
angle with sufficient accuracy to make use 
of the steep slope of the upper curves of 
Fig. 1. An extremely precise reading of scan 
angle might permit the measurement of 
plasma densities which are beyond present 
techniques. 


APPLICATION 


The type of experiment in which this 
technique might be useful is suggested by 
Fig. 2. The transmitting and receiving an- 
tennas are electronically phased arrays. 
The phasing is synchronized so that the 
transmitting and receiving patterns always 
point in a common direction. If the phasing 
is cycled at a high rate compared to the 
plasma decay, the received signal can record 
the progress of the resonance with scan 
angle as the plasma density diminishes. Al- 
ternately, individual detectors can be placed 
on each element of the receiving array. 
Angular sensitivity can be increased by the 
use of conventional sum-and-difference tech- 
niques. 
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Fig. 2—Possible experimental setup. 


CONCLUSIONS 


Four resonances are possible in the inter- 
action of an electromagnetic wave and a 
plasma in the presence of a dc magnetic 
field. Three of these are zeros and are unaf- 
fected by the angle between the Poynting 
vector and the dc magnetic field. The fourth 
is a pole and is angle sensitive. For values 


_of the cyclotron frequency above the plane- 


‘wave frequency, big dynamic ranges occur 


in the resonant value of plasma density as 


* 


the angle is changed. Thus, a scanning beam 
offers the possibility of becoming a useful 
diagnostic tool. 
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_ The Zone Plate as a Radio- 


Frequency Focusing Element* 
INTRODUCTION 


More than a hundred years ago, Fresnel 


invented the optical zone plate, an interfer- 
ence device which has focusing properties 


similar to those of a lens. Application of the 
device seems to have been limited, and it has 
remained a laboratory curiosity. In the 
course of design studies on large parabolic 
reflectors for use in satellite communication, 
investigation of the zone plate was prompted 
by the apparent simplification it would of- 
fer in construction, in that it can be con- 


* Received by the PGAP, January 20, 1961. 
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structed on a plane surface. This note de- 
scribes some of the properties of zone plates 
as potential substitutes for large paraboloids. 


Basic THEORY 


The zone plate consists of a series of con- 
centric opaque and transparent rings. 

According to Huygen’s principle, each 
point of a transparent zone acts as a second- 
ary source which emits a spherical wave. 
Let the point P be the desired focal point. 
Then, if the zones are of such a diameter 
that radiation from all the transparent zones 
arrive in phase, there will be constructive 
interference. This requires that the distance 
from P to the edges of the successive zones 
be integral multiples of a half-wavelength. 
That is, 


Beet toy (1) 


where py is the distance from the focal point 
to the Nth zone edge, F is the focal length, 
and \ is the wavelength. Therefore, the 
radius of the Nth edge is 


ty = Vey? — F?, (2) 
or 
Nx if AF B) 
. <= ae 
aS 2 NX 


Eq. (3) may be used for calculating the radii 
of the zones. However, we have found it 
usually simpler and quicker to use a graph- 
ical layout. 

It is interesting to note that the focusing 
action is indifferent to whether the first 
zone is opaque or transparent. That is, the 
actions of complementary zone plates are 
identical. = 

This means that a given zone plate has 
a focus for reflection as well as transmission, 
since the opaque rings can be considered as 
reflection sources if they are conducting. If 
two feed antennas are located at these foci, 
and are properly phased, it should be pos- 
sible to utilize both foci. 


THEORETICAL GAIN 


We now attempt to compute the theo- 
retical gain of the zone plate. The amplitude 
Am, due to the mth zone, is proportional to 
the area of the zone and inversely propor- 
tional to the distance from the zone to the 
focal point. Also, we must include an “ob- 
liquity factor” in Huygen’s principle. The 
resulting equation is 


Am = K * (1 + cos6), (4) 


m 


where dm is the area’of the mth zone, dm is 
the average distance from the zone to the 
focal point, (1+ cos @) is the obliquity factor 
and K is a proportionality factor to be de- 
termined. 

It requires only a bit of algebraic manip- 
ulation to derive the theoretical gain formu- 
las. 

Expressed in decibels, the gain for the 
aperture center zone plate is 
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Guiad = 20 log ey Saal 
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and for the obstacle center zone plate, the 
gain is 


Ga.,obs = 20 log 


where M is the number of zones in the plate. 

Strictly speaking, (5) and (6) are correct 
only for zone plates imbedded in the center 
of an infinite screen. However, for a reason- 
ably large number of rings, the obliquity 
factor comes into play, so that the outer 
zones become less and less important. 

It is interesting to compare these theo- 
retical gains with those of a paraboloid. As a 
specific example, the gain of a zone plate 
with F/A=10 has been calculated as a func- 
tion of the number of zones. The correspond- 
ing diameters were then computed. The re- 
sults are plotted in Fig. 1. 
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Fig. 1—Comparison of theoretical gains for 
zone plates and paraboloids. 


Use oF A GROUND PLANE 


As mentioned earlier, any zone plate has 
two principal foci, one for transmission and 
one for reflection. It therefore seems natural 
to attempt to fold the system with a ground 
plane and bring the two foci into coincidence. 

The radiation that is reflected from the 
opaque zones will have undergone a 180° 
phase change. The radiation that passes 
through the transparent zones travels a 
distance d to the ground plane, undergoes a 
phase reversal upon reflection, and returns 
through the distance d. In order for the radi- 
ation from both portions to arrive at the 
focal point in phase, the two portions should 
be 180° out of phase at the plane of the zone 
plate. Therefore, the spacing d should be a 
quarter wavelength. 
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EXPERIMENTAL PROCEDURE 


Experimental zone plates were built for 
the X band by cementing aluminum foil to 
plexiglass sheets, then laying out the rings 
with a beam compass, cutting along the 
lines with a razor blade, and finally peeling 
off the unwanted foil. Fig. 2 shows such a 
plate mounted on the antenna range turn- 
table. 


Fig. 2—Zone plate with ground plane mounted 
on antenna range turntable. 


Two types of feed were tried—a straight 
dipole above a ground plane, and a turnstile 
in which the arms were drooped toward the 
ground plane in an effort to somewhat widen 
the pattern of the feed. 

Fig. 3 is an example of the pattern ob- 
served with a transmission plate of three 
zones, open aperture center. Both transmis- 
sion and reflection patterns are shown. Note 
' that the beam width is about 4° in each case. 
The diameter of this plate was about 16.5 
wavelengths. The beamwidth is almost iden- 
tical with that of a paraboloid having the 
same diameter. And this is the salient feature 
of zone plates. Even though the gain is much 
lower than that of a paraboloid of the same 
diameter (as shown in Fig. 1), the beam- 
width is the same. Therefore, in applications 
where a narrow beamwidth is required and 
the gain can be achieved by amplification, 
it becomes feasible to use the zone plate. 

The reason for the partial independence 
of gain and beamwidth is in the very nature 
of the zone plate as a diffraction device. 
Instead of scattering all of the incident 
radiation into the feed like a paraboloid, the 
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Fig. 3—Zone-plate patterns showing transmission 
and reflection maxima. 
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Fig. 4—Pattern of zone plate with ground 
plane (folded system), 


zone plate scatters only a portion of it into 
the feed and hence, has a much lower gain. 
Beamwidth, on the other hand, is deter- 
mined by the physical diameter. 

Referring again to Fig. 3, we note that 
the sidelobe level on transmission is about 
15 db, whereas on reflection there is a pair of 
near sidelobes about 8 db below the peak. 
These sidelobes have not been explained 
theoretically, although they may be due to 
obscuration by the feed antenna. 

When a ground plane is added and the 
transmission and reflection foci are made to 
coincide, the pattern is as shown in Fig. 4. 
Here again, near sidelobes appear about 7 
db down. 

The theoretical gain equations were com- 
pared with experimental data by sequen- 
tially covering the open zones of a three- 
zone plate. Table I compares experimental 
vs theoretical values. 


TABLE I 


CoMPARISON OF EXPERIMENTAL AND THEORETICAL 
GAINS IN RELATION TO NUMBER OF ZONES 
(APERTURE CENTER ZONE PLATE) 


Number Theoretical Gain Experiment 
of Zones Over Feed (db) Gain (db) 
1 6.0 4.5 
2 11.8 9.0 
Kj 1552 11.5 


CONCLUSION 


Certain aspects of zone-plate perform- — 
ance have been investigated: Because of the 
constructional simplicity, it appears that in 
certain applications where large-aperture 
focusing elements are required, the zone 
plate may prove useful. Such applications 
include low-cost radio telescopes and space 
communications antennas. 


Fig. 5—Fixed zone-plate radio telescope. 


Fig. 5 shows a possible configuration 
based on the folded system. Here the zone 
plate is fixed in position above the ground 
plane, and only the feed antenna is moved 
about. Owing to the deterioration of per- 
formance at large scan angles, such a system 
is limited to rather small angles away from 
the zenith, perhaps + 20°. However, the cost 
of constructing such a device would be quite 
low, compared to even a comparable para- 
boloid. It is felt that the system might be 
useful, for example, to a small college desir- 
ing to do some work in radio astronomy. 

L. F. VAN BUSKIRK 
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Michelson Lab. 

U. S. Naval Ordnance Test Station 
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